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MODELS  AND  COMPUTATIONAL  METHODS 


FOR  DYNAMIC  FRICTION  PHENOMENA 


J.T.  ODEN  and  J.A.C.  MARTINS 

ABSTRACT 


This  work  addresses  the  general  problem  of  formulating  continuum  models 
of  a  large  class  of  dynamic  frictional  phenomena  and  of  developing  computation 
methods  for  analyzing  these  phenomena.  Of  particular  interest  are  theories 
which  can  adequately  predict  stick-slip  motion 3  frictional  damping  in  structural 
dynamics 3  and  sliding  resistance.  This  work  is  divided  into  three  principal 
parts .  In  Part  Is  a  large  body  of  experimental  and  theoretical  literature  on 
friction  is  critically  reviewed  and  interpreted  as  a  basis  for  models  of  dynamic 
friction  phenomena.  In  Part  II 3  continuum  models  of  interfaces  are  developed 
which  simulate  key  interface  properties  identified  in  Part  I.  Variational  prin¬ 
ciples  for  a  class  of  dynamic  friction  problems  are  also  established.  In  Part 
IIIj  finite  element  models  and  numerical  algorithms  for  analyzing  dynamic  fric¬ 
tion  are  presented.  Also3  a  dynamic  stability  analysis  is  presented  in  which 
it  is  established  that  stick-slip  motion  can  be  associated  with  dynamic  insta¬ 
bility  of  the  governing  nonlinear  system  for  certain  ranges  of  slip  velocity  and 
coefficient  of  friction.  Numerical  results  suggest  that  the  new  models  derived 
here  can  satisfactorily  depict  a  large  and  important  class  of  dynamic  friction 
effects . 


I .  INTRODUCTION 


The  contact  and  rubbing  of  solid  bodies:  to  many  it  must  seem  remarkable 
that  an  event  so  common,  so  intrinsic  to  the  mechanics  of  everyday  life,  so 
important  in  a  multitude  of  applications  of  mechanics  to  engineering  problems,  and 
so  often  the  subject  of  experimental  research,  has  not  been  satisfactorily 
depicted  by  a  sound  continuum  model  to  date.  However,  to  those  who  have  taken 
more  than  a  superficial  look  at  the  subject,  the  absence  of  a  universal  continuum 
model  of  friction  may  not  be  surprising.  The  nature  of  dynamic  friction  forces 
developed  between  bodies  in  contact  is  extremely  complex  and  is  affected  by  a  long 
list  of  factors:  the  constitution  of  the  interface,  the  time  scales  and  frequency 
of  the  contact,  the  response  of  the  interface  to  normal  forces,  inertia  and 
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thermal  effects,  roughness  of  the  contacting  surfaces,  history  of  loadings,  wear 
and  general  failure  of  the  interface  materials,  the  presence  or  absence  of 
lubricants,  and  so  on.  Thus,  dynamic  friction  is  not  a  single  phenomenon  but  is  a 
collection  of  many  complex  mechanical  and  chemical  phenomena  entwined  in  a  mosiac 
whose  features  cannot  be  grasped  through  isolated  simple  experiments. 

Successful  computational  methods  in  nonlinear  mechanics  are  generally  those 
based  on  a  thorough  familiarity  with  the  natural  phenomena  being  simulated  and  on 
a  good  understanding  of  the  mathematical  models  available  to  describe  it.  A  major 
issue  in  the  computer  simulation  of  many  natural  events  is  the  absence  of  a 
physically  and  mathematically  sound  model  of  the  phenomena.  This  may  not  be  a 

critical  issue  in  large  areas  of  fluid  and  gas  dynamics  where  the  models 
characterized  by  the  classical  Navier-Stokes  equations  or  the  Euler  equations  are 
regarded  as  acceptable,  but  it  is  an  especially  important  issue  in  many  areas  of 
nonlinear  solid  mechanics,  where  debates  still  rage  over  appropriate  models  of 
finite  elastoplasticity ,  thermomechanical  phenomena,  etc.  The  study  of  friction 
effects  in  solid  mechanics  is  particularly  hampered  by  the  absence  of  sound, 
generally  accepted,  models. 

The  recognition  that  it  may  not  be  reasonable  to  attempt  to  describe  all 
dynamic  friction  effects  by  a  single  theory  seems  to  be  a  major  step  forward. 
Just  as  a  theory  of  linear  elasticity  may  provide  an  appropriate  model  of  the 
behavior  of  mild  steel  only  for  a  limited  range  of  deformations,  so  also  can  any 
single  model  of  friction  be  expected  to  be  valid  only  under  a  limited  set  of 
conditions.  Thus,  we  shall  focus  here  on  a  limited  but  very  important  class  of 
friction  problems:  the  dynamic  contact  and  relative  sliding  of  dry, metallic 

surfaces  with  engineering  finishes  (ground,  abraded,  scraped,  polished,  etc.) 

Still  further  restrictions  are  required.  We  choose  to  divide  this  general  class 
of  frictional  effects  into  three  general  categories: 


■*  fan*  wt*  Mrm  \r»  vm*  tej-w  *r*  wv  .*r»  w**  \tSm  w**  ..  *  >  .  ■»  wnr  vw  vs*  V  >  -  >  u  *  u>  v  >  v  ►  V  v  >  v  \  v  f  \  v*  m;\>-  /VKA KM. 


I.  Quasi-static  Dry  Friction 

II.  Dynamic,  Sliding  Friction 

III.  Wear  and  Plowing 


3 


While  these  classifications  are  by  no  means  sharp  and  exclusive,  they  are  adequate 
for  our  present  purposes. 

By  quasi-static  dry  friction  (Type  I-Friction),  we  mean  frictional  mechanisms 
present  when  two  possibly  polished,  metallic  surfaces  are  pressed  slowly  together 
and  ire  in  static  equilibrium  or  are  slowly  displaced  relative  to  one 
another,  and  for  which  normal  loads  are  sufficiently  small  that  no  severe  penetra¬ 
tion  and  gouging  of  the  contact  interface  occurs.  Here  the  most  conspicuous 

mechanisms  contributing  to  friction  are  the  plastic  deformation  of  asperities, 
the  formation  of  elastoplastic  junctions  on  the  contact  surface,  the  strong 
coupling  between  normal  and  tangential  plastic  deformations,  and  the  dissipative 
mechanisms  afforded  by  the  history-dependence  of  junction  deformation  and  possible 
fracture  of  these  junctions.  The  frictional  forces  may  depend  upon  histories  of 
micro-tangential  displacements  of  particles  on  the  contact  surface.  Theories 
describing  such  frictional  effects  are  basically  static  in  nature;  inertia  effects 
are  generally  not  important  and  the  modeling  should  capture  the  elastoplastic 
deformation  of  the  interface  in  some  way.  Models  and  computational  procedures 
for  studying  these  classes  of  problems  have  been  advanced  by  several  investigators, 
and  we  mention  as  examples  the  works  of  FREDRICKSON  [1976],  ODEN  AND  PIRES  [1983], 

PIRES  and  ODEN  [1983],  MICHALOWSKI  and  MROZ  [1978],  and  CAMPOS,  ODEN,  and  KIKUCHI  [1982]. 


By  dynamic  sliding  friction  (Type  II-Friction) ,  we  refer  to  that  large  and 
important  class  of  truly  dynamic  problems  which  include  such  effects  as  frictional 
damping,  dynamic  sliding,  stick-slip  motion,  chattering,  etc.  The  constitution  of 
the  material  interface  is  essentially  stable;  there  is  no  marked  penetration  or 
normal  plastic  deformation  of  the  interface  and,  at  least  from  a  global  point  of 
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view,  the  frictional  forces  developed  on  the  contact  surface  appear  to  depend  on 
the  sliding  velocity  of  one  surface  relative  to  another.  To  obtain  reproducible 
experimental  data  on  these  types  of  frictional  effects,  it  is  generally  necessary 
to  rub  the  contacting  surfaces  together  for  a  period  of  time  until  a  steady-state 
condition  of  the  interface  is  reached.  Such  pre-conditioning  may  fracture  and 
work-harden  the  surface  asperities  and  result  in  a  smoother  contact  surface  for 
which,  unlike  category  I,  plastic  deformations  normal  to  the  interface  are  not 
significant  or  do  not  occur  at  all. 

A  significant  distinguishing  feature  between  Type  I  and  II  Friction  is  that 
of  scales:  In  Type  I,  the  local  frictional  forces  on  the  interface  depend  upon 
surface  micro-displacements  which  may  be  of  an  order  of  magnitude  of  ten  microns; 
in  Type  II,  much  larger  rigid  relative  motions  can  occur  and  frictional  forces 
vary  with  sliding  velocity. 

The  wear  and  plowing  phenomena  (Type  III)  encompasses  cases  in  which 
substantial  damage  and  penetration  of  the  metallic  interfaces  occurs;  there  may 
be  phase  changes  in  the  interface  materials,  and  portions  of  the  material  may  be 
removed  or  relocated  during  contact.  These  types  of  contact  situations  may 
involve  the  transfer  of  large  normal  forces  across  a  small  contact  area.  Wear 
may  be  experienced  in  all  three  of  the  categories  mentioned,  but  in  categories 
I  and  II  it  is  either  ignored  or  assumed  to  be  negligible.  We  shall  not  consider 
these  types  of  frictional  phenomena  in  the  present  work.  We  also  do  not  consider 
here  problems  of  collision  and  impact,  which  differ  from  those  mentioned  earlier 
due  to  the  occurrence  of  significant  jumps  in  particle  velocities  at  the  contact 
surface,  an  instantaneous  transfer  of  momentum  from  one  body  to  another,  and  the 
propagation  and  reflection  of  stress  waves.  Moreover,  while  all  contact  problems 
with  friction  involve  thermal  effects  to  varying  degrees,  we  shall  ignore  such 
effects  in  the  present  investigation. 
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The  present  paper  is  a  study  of  models  and  computational  procedures  for 
analyzing  contact  and  friction  problems  that  fall  into  category  II:  sliding 
friction  and  stick-slip  motions  of  dry,  rough,  metallic  surfaces. 

The  first  and  perhaps  most  difficult  step  in  the  analysis  of  dynamic  fric¬ 
tion  is  the  development  of  an  acceptable  model  of  the  contact  interface.  It 
is  well  known  that  the  interface  through  which  bodies  contact  is  a  mechanically- 
complicated  media,  the  constitution  of  which  is  different  from  that  of  the 
parent  material  of  either  contacting  body.  The  interface  may  be  a  rough  ir¬ 
regular  pair  of  surfaces  composed  of  impurities,  oxides,  work-hardened  or 
crushed  materials,  gases  and  lubricants.  TheAe^OAe,  It  Is  natuAal  In  develop¬ 
ing  phe.nome.ho logical  models  o£  fcnlctlon,  to  assign  to  the  InteA&ace  a  separate 
stmictune  chaA.acteJvlzed  by  constitutive  equations  Independent  ofi  those  chaAac- 
tentzlng  the  paAent  metals. 

Toward  assessing  what  features  these  interface  models  should  exhibit,  we 
devote  the  first  part  of  this  study  to  a  review  and  critique  of  a  substantial 
body  of  experimental  literature  on  this  subject.  Then,  on  the  basis  of  our 
interpretations  of  this  data,  we  propose  forms  of  constitutive  equations  for 
the  contact  interface,  particularly  its  response  to  normal  stress. 

In  part  II,  we  incorporate  these  dynamic  friction  effects  into  a  nonlinear 
continuum  model  of  contact  and  sliding  friction  of  elastic  bodies.  The  con¬ 
tacting  bodies  are  assumed  to  be  linearly  elastic,  but  the  overall  theory  is 
highly  nonlinear,  owing  to  nonlinearities  in  the  contact  constraints,  fric¬ 
tional  behavior,  and  interface  response.  We  derive  new  variational  principles 


governing  these  dynamic  problems. 
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Part  III  of  the  study  is  devoted  to  the  numerical  analysis  of  the  dynamic 
friction  models  developed  in  Part  II#  There  we  describe  several  numerical  al¬ 
gorithms  for  studying  dynamic  friction.  Dynamic  stability  and  other  qualita¬ 
tive  features  of  sliding  friction  phenomena  are  studied  numerically  using  a 
simple  mechanical  model  and  finite  element  discretizations  of  our  continuum 
model.  Several  representative  problems  are  solved  numerically. 

There  are  several  results  of  this  study  which  we  feel  are  particularly 
significant: 

1.  Contrary  to  a  widely  accepted  belief  that  has  prevailed  since  the  time 
of  Coulomb 3  the  coefficient  of  friction  may  not  necessarily  change  with  increas¬ 
ing  relative  sliding  velocity.  Our  results  confirm  the  experimentally -based 
conjectures  of  TOLSTOI  {1967]  and  some  other  researchers.  This  apparent  change 
is  traditionally  interpreted  as  a  decrease  from  a  static  to  a  kinetic  coefficient. 
It  is3  however 3  only  the  average  value  of  the  friction  force  that  may  decrease 
after  the  initiation  of  sliding  and  not  necessarily  instantaneous  ratios  tangential 
to  normal  stress  components  on  the  contact  sir  face.  Of  course 3  if  the  sliding 
body  is  modeled  as  a  single  ( tangential )  degree-of- freedom  system 3  as  is  so 

often  the  case3  then  the  reduction  of  the  coefficient  of  friction  upon  sliding 
is  the  only  possible  devise  for  incorporating  this  experimentally  observed 
phenomena.  Such  crude  models3  experimental  and/or  analytical 3  cannot  take  into 
account  normal  force  oscillations 3  and  thus  omit  a  critically  important  prop¬ 
erty  of  dynamic  friction. 

2.  A  key  feature  of  acceptable  interface  models  for  dynamic  friction  is 
the  characterization  of  the  response  of  the  interface  to  normal  forces.  This 
mechanical  response  for  most  metal-on-metal  interfaces  is  highly  nonlinear. 
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3.  Stick-slip  motion  may  be  a  manifestation  of  dynamic  instabilities  inherent 
in  the  coupling  of  normal  and  tangential  relative  motions  of  contacting  bodies . 

This  phenomena  is  not  necessarily  the  result  of  a  decrease  in  the  coefficient  of 
friction  with  changes  in  sliding  velocity,  and,  can  in  fact,  be  observed  when  the 
coefficient  >f  friction  is  constant  and  equal  to  its  so-called  static  value. 

4.  It  is  possible  to  develop  and  to  incorporate  simple  nonlinear  interface 
models  into  a  continuum  model  of  elastodynamics  which  captures  many  fundamental 
features  of  dynamic  friction.  In  particular,  frictional  damping  can  be  depicted 
in  such  models  in  a  natural  and  straightforward  way . 

5.  Finite  element  approximations  of  the  continuum  models  can  be  developed 
which  feature  consistently-derived  frictional  damping/ stiffness  matrices.  These 
finite  element  methods,  together  with  numerical  schemes  for  solving  associated 
systems  of  nonlinear  ordinary  differential  equations,  are  capable  of  modeling 
stick-slip  motion,  dynamic  sliding,  friction  damping,  and  related  phenomena  in 

a  significant  range  of  practical  problems. 

PART  I.  PHYSICAL  ASPECTS  OF  DYNAMIC  FRICTION 

This  portion  of  our  study  is  devoted  to  a  review  of  a  selected  body  of 
literature,  mostly  experimental  in  nature,  that  deals  with  the  subject  of  dry 
friction  on  metallic  surfaces.  It  is  necessarily  an  imcomplete  review,  since 
this  subject  is  an  old  and  active  one  in  which  a  large  volume  of  literature 
exists;  but  we  have  attempted  to  focus  on  a  selected  sample  of  results  that 
collectively  provide  some  basis  for  friction- interface  models  to  be  proposed 


later 
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2.  SOME  PROPERTIES  OF  METALLIC  SURFACES 
2.1  Preliminaries 

We  begin  with  a  quotation  from  WHITEHOUSE  and  ARCHARD  11970]: 

"All  surfaces  are  rough.  This  is  the  starting  point  from  which 
current  ideas  about  friction,  wear,  and  other  aspects  of  surfaces  in 
contact  have  evolved" . 

It  is  well  known  that  the  most  polished  metallic  surfaces  are  not  perfectly 
flat.  Under  magnification,  one  observes  that  these  surfaces  have  undulations 
that  form  hills  and  valleys  the  dimensions  of  which  are  large  in  comparison 
with  molecular  dimensions.  These  deviations  from  the  plane  are  called  asperities. 

Many  experimental  techniques  have  been  used  to  observe  and  measure  the 
roughness  of  the  surfaces.  Summary  accounts  of  these  techniques  and  references 
for  more  detailed  descriptions  can  be  found,  e.g.,  in  BOWDEN  and  TABOR  [1950], 
MOORE  [1975]  and  THOMAS  [1975].  Recent  developments  can  be  found  in  STOUT  and 
THOMAS  [1982].  A  popular  instrument  to  measure  surface  roughness  is  the  Taylor- 
Hobson  Talysurf,  with  which  the  vertical  movement  of  a  stylus  moving  over  a  re¬ 
presentative  length  of  the  surface  under  examination  is  amplified  electrically  . 
and  recorded  on  a  moving  paper  chart.  Since  thB  length  traversed  by  the  stylus 
is  very  large  compared  with  the  height  of  the  asperities,  a  different  scale  is 
used  for  the  horizontal  and  vertical  magnifications  (see  Fig.  1). 

Depending  on  the  method  of  production  of  the  surface,  the  height  of  the 
peaks  may  vary  between  0.05ym  to  50ym  while  the  spacings  between  them  range 
from  0.5ym  to  5mm  (SARKAR  [1976]).  The  average  slope  of  the  asperities  is 
usually  very  small,  of  the  order  of  5°  to  10°  (TABOR  [1981]). 

Various  quantities  have  been  proposed  as  measures  of  surface 
roughness.  For  example,  let  z(x)  denote  the  actual  height  of  particles 
on  a  material  surface  about  a  fixed  datum  plane  at  a  point  x  on  the 


Figure  1.  Profiles  of  mild  steel  specimen  after  three  surface 
treatments,  (a)  Surface  ground  only,  (b)  Surface 
ground  and  then  lighly  polished,  (c)  Surface  ground 
and  then  lightly  abraded  on  600  paper  (GREENWOOD  and 
WILLIAMSON  [1966]). 
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plane  (Indicated  in  one  dimension  in  Fig.  2).  The  center-line — average 
of  asperity-heights  over  a  test  length  L  is  obviously 

T  /L|*(x)|dx 
0 

Ue  can  relocate  this  datum  on  the  mean  profile  plane,  and  compute  as 
roughness  measures  the  r.m.s.  ("root-mean-square”)  deviation  o  ,  the 
r.m.s.  of  derivatives  a*  ,  r.m.s.  of  second  derivatives  a  ",  the  auto-corre- 
lation  C  or  the  power  spectral  density  P  : 


C(A) 


z(x  +  A)dx 


CO 

/  P(X)exp(iXA)dX 
—00 


P(X)  »  J 


00 


— CO 


C(A)exp(-iXA)dA 
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Here  A  is  the  delay  interval  and  X  is  the  angular  frequency  of  asperities. 

The  parameters  a*  and  a"  are  regarded  as  sensitive  indicators  of  the  slope 
of  the  asperities  and  of  the  degree  of  sharpness  at  asperity  peaks,  respectively. 

The  .auto  correlation  function  of  a  profile  C(A)  describes  the  extent  to 

which  (on  an  average)  the  height  at  any  point  on  the  profile  is  correlated  with 

the  height  at  another  point  downstream.  For  a  random  signal,  C(A)  declines 
2 

from  a  at  A  «  0  to  zero,  as  A  becomes  large.  The  power  spectral  density 
function  P(X)  is  the  Fourier  transform  of  C(A)  and  represents  the  rate  of 
change  of  mean  square  value  with  frequency,  where  the  mean  square  value  is  taken 
in  a  narrow  frequency  band  at  various  center  frequencies.  For  this  reason  P(X) 
is  also  known  as  the  mean  square  spectral  density. 

2.2  Constitution  of  Metallic  Surfaces. 

The  frictional  properties  of  unlubricated  metals  are  greatly  affected 
by  the  presence  of  surface  films.  Figure  3  illustrates  schematically  the 
surface  layers  encountered  on  metallic  surfaces  in  industrial  environments 
(RABINOWICZ  11965]).  According  to  RABINOWICZ,  the  contaminant  layer  may  be 
derived  from  a  variety  of  sources,  among  them  the  oil  drops  found  in  industrial 
environments  and  the  lubricants  applied  during  the  preparation  of  the 

surface.  When  in  contact  with  the  air,  the  main  constituents  of  the 

adsorbed  layer  are  generally  molecules  of  water  vapor  and  oxygen.  The 
oxide  layer  is  produced  by  reaction  of  oxygen  from  the  air  with  all 
metals  (except  with  the  noble  metals,  such  as  gold  and  platinum,  which 
are  generally  free  of  oxide  or  other  surface  films).  The  wo rk-hardened 
layer  results  from  the  heavy  deformation  of  the  metal  due  to  the 

mechanical  preparation  of  the  surface. 

According  to  BOWDEN  and  TABOR  {1964],.  if  the  preparation  of  the 

surface  involved  polishing  at  high  speeds,  the  surface  layer  may  consist 
of  a  smeared  "fudge"  of  metal,  metal  oxide,  and  polishing  powder.  Below 
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Contaminant  layer 
Adsorbed  gas  layer 
•Oxide  layer 
Work- hardened  layer 

Metal  substrate 


Figure  3.  Schematic  illustration  of  films  on  a  metal  surface 
(not  to  scale)  (RABINOWICZ  [1965]), 
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the  severely  deformed  layers,  the  metal  usually  shows  vestigial  traces  of 
deformation  produced  in  earlier  stages  of  abrasion.  A  schematic 
representation  of  this  is  shown  in  Figure  4. 

3.  THE  NORMAL  CONTACT  OF  METALLIC  BODIES 


3.1.  Static  Contact  of  Rough  Surfaces 

In  this  section  we  summarize  several  theoretical  and  experimental  results 
on  the  static  normal  contact  of  rough  surfaces.  For  the  theoretical 
developments  we  follow,  essentially,  the  survey  papers  of  ARCHARD  [l974]  , 
THOMAS  [l975j  and  WHITEHOUSE  [l98o].  For  the  experimental  work  we  follow, 
essentially,  the  survey  papers  of  BACK,  BURDEKIN  and  COWLEY  [l973]and 
WOO  and  THOMAS  [l98o].  We  are  especially  interested  in  answers  to  the 
following  questions:  What  is  the  stiffness  of  two  asperity-covered  sur¬ 
faces  when  pressed  against  each  other?  What  is  the  corresponding  true 
area  of  contact? 

3-1.1  Various  Models  of  Contacting  Rough  Surfaces 

Because  surfaces  are  rough,  the  true  area  of  contact  is  much  smaller 
than  the  apparent  area  of  contact.  As  a  result,  the  true  contact  surfaces 
may  often  support  pressures  so  large  that  they  are  comparable  with  the 
strengths  of  the  materials  of  the  contacting  bodies.  BOWDEN  and  TABOR 
[l95o]  suggested  that  these  pressures  would  be  always  high  enough  to  produce 
plastic  flow.  The  contacting  regions  would  deform  plastically  until  the 
true  area  of  contact  would  be  just  large  enough  to  support  the  load.  The 
true  area  of  contact  A^  would  be  then  proportional  to  the  applied  load 
N  .  With  this  proportionality,  Ar«^N,  it  is  possible  to  provide  a  simple 
and  elegant  explanation  of  Amontons*  laws  of  friction  (see  Sections  4.1,2). 

Objections  to  the  concept  of  plastic  deformation  of  the  asperities 
were  put  forward  by  ARCHARD  [l 957] .  According  to  this  author,  when 
materials  of  comparable  hardness  rub  together > a  protuberance  may  be  plas- 
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Figure  4.  Schematic  diagram  showing  topography  and  structure 
of  a  typical  polished  metal  specimen  (BOWDEN  and 
TABOR  [1964]). 
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tically  deformed  at  its  first  encounter  with  the  other  surface,  but  its 

relaxation  would  be  elastic.  At  its  many  subsequent  encounters  with 

the  other  surface  the  protuberance  would  bear  the  same  load  by  elastic 

deformation.  ARCHARD  shows  then  that  although  the  simple  Hertzian 

theory  does  not  predict  the  proportionality  between  A^  and  N  ,  a 

generalized  model  in  which  each  asperity  is  covered  with  microasperities, 

and  each  microasperity  with  micro-microasperities,  gives  successively 

closer  approximations  to  the  law  A^o-cn  as  more  stages  are  considered 

(Figure  5).  Archard  explained  that  the  essential  part  of  the  argument 

was  not  the  choice  of  asperity  model:  it  was  whether  an  increase  in  load 

creates  new  contact  areas  or  increases  the  size  of  existing  ones  •  For 

physically  plausible  surfaces  any  elastic  model  in  which  the  number  of 
.  2/3 

contacts  remains  constant,  we  find  that  A^-<  N  ;  but  if  the  average  size 
of  contacts  remains  constant  (and  its  number  increases)  then  A^^c  N.  Al¬ 
though,  admittedly  artificial,  Archard' s  models  were  thus  very  important 
in  showing  the  admissibility  of  elastic  deformation  of  the  asperities  and 
the  effect  of  the  superposition  upon  the  surfaces  of  asperities  of  widely 
differing  scales  of  sizes. 

The  next  step  came  with  the  development  of  surface  models  based  on 
knowledge  gained  from  the  examination  of  the  actual  surface  topography. 
GREENWOOD  and  WILLIAMSON  [  1 9 6 6 J ,  and  others,  have  shown  that,  for  many 
surfaces,  the  distribution  of  heights  is  very  close  to  Gaussian.  These 
authors  also  investigated  the  distribution  of  peak  heights  and  concluded 
that  it  was  also  close  to  Gaussian  (see  Fig.  6).  The  model  developed  by 
GREENWOOD  and  WILLIAMSON  thus  assumed  a  Gaussian  distribution  of  peak 
heights  (asperity  heights).  They  also  assumed  that  the  asperities,  at 
least  near  their  summits,  were  spherical  with  a  constant  radius  of  cur¬ 
vature.  The  topography  of  the  surface  was  described  in  terms  of  three 
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Figure  6.  Cumulative  height  distribution  of  bead-blasted  aluminum. 
(x)  Distribution  of  all  heights 
(•)  Distribution  of  peak  heights 

A  perfect  Gaussian  distribution  would  give,  with  the 
vertical  scale  used,  a  straight  line. 

(GREENWOOD  and  WILLIAMSON  [1966]). 
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parameters:  a  ,  the  standard  deviation  of  the  distribution  of  asperity 
heights,  R  ,  the  mean  radius  of  curvature  of  the  asperity  tips  (assumed 
constant  in  the  model),  and  n  »  the  surface  density  of  asperities  (later, 
WHITEHOUSE  and  ARCHARD  [1970]  showed  that  these  parameters  are  not  in¬ 
dependent). 


The  deformation  of  the  asperities  was  assumed  to  be  elastic  (Hertzian). 
The  material  property  required  was  thus  the  equivalent  elasticity  modulus 


of  the  two  surfaces  E'  ,  defined  as 
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1-v 


1 


i-v; 
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(3.1) 


where  Ej  and  E£  are  the  moduli  of  elasticity  of  each  of  the  contacting 

bodies  and  Vj  and  \>2  are  the  corresponding  Poisson's  ratios. 

The  results  obtained  using  this  theory,  for  typical  values  of  a  , 

R  ,  r)  and  E'  are  shown  in  Figures  7  and  8.  In  Figure  7,  the  separation 

* 

h  is  the  distance  between  the  mean  planes  of  the  distributions  of  peak 
heights  of  the  two  surfaces.  It  is  clear  that  the  separation  decreases 
approximately  proportionally  to  the  increase  of  the  logarithm  of  the 
normal  load  and  that  the  real  area  of  contact  is,  approximately  proportional 
to  the  normal  load.  GREENWOOD  and  WILLIAMSON  observed  that  the  above 
proportionalities  hold  exactly  if  the  distribution  of  asperity  heights  is 
assumed  exponential  rather  than  Gaussian.  For  such  a  distribution,  which 
is  a  fair  approximation  to  the  uppermost  25  percent  of  the  asperities  of 
most  surfaces,  it  is  possible  to  obtain  closed  form  expressions  for  the 
relations:  load-separation,  area  of  contact  separation  and,  consequently, 
area  of  contact-load.  These  relations  are  the  following: 


N  (Kg) 


Figure  8.  Relation  between  true  area  of  contact  and  load. 

2 

nominal  area  A  =  10  cm 

2 

-  nominal  area  A  =  1  cm 


(GREENWOOD  AND  WILLIAMSON  ,[1966j). 
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j$-  *  W  (n^a*)  Va*/R  exp(-h*/a*) 

A  ^  ^  ^ 

■  it  (nRo  )  exp(-h  /a  )  (3.2) 

,  YiT  •J'R-/a*  P 
A  *  E* 


where  p  is  the  nominal  contact  pressure,  A  is  the  nominal  area  of 
contact  and  Ar  is  the  real  area  of  contact. 

TSUKIZOE  and  HISAKADO  [l965,  1968]  developed  a  statistical  model 
for  the  contact  of  rough  surfaces  under  the  assumption  of  plastic  defor¬ 
mation  of  the  asperities.  They  further  assumed  a  Gaussian  distribution 
of  surface  heights  and  that  the  material  displaced  due  to  plastic  deforma¬ 
tion  could  be  neglected.  This  approach,  which  works  well  for  light  to 
moderate  loads,  leads  to  the  following  relations  involving  load,  area  of 
contact  and  separation: 


j  erfc  ( 
■  j  erfc 
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(3.3) 


Here  the  separation  h  is  the  distance  between  the  mean  planes  of  the 
distributions  of  surface  heights  (see  Fig.  9),  and  H  is  the  indentation 
hardness.  (H  =  applied  load/permanent  indentation  area  when  a  hard  in- 
denter  is  pressed  into  a  "flat"  surface.) 

Some  of  the  limitations  of  these  early  statistical  models  of  elastic 
or  plastic  contact  have  been  successively  overcome  by  more  sophisticated 
models  proposed  by  other  authors. 

WH1TEH0USE  and  ARCHARD  11970]  used  the  fact  that,  if 


current  mean  line  of  the  rough  surface 

mean  line  of  the  rough  surface  at  the  beginning 
of  contact 

current  separation 

separation  at  the  beginning  of  contact 
penetrating  approach 

Separation  and  penetrating  approach  between  a  rough 
and  a  flat  surface. 
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a  surface  profile  is  of  a  random  type,  then  it  can  be  defined  completely 
(in  a  statistical  sense)  by  two  characteristics:  the  probability  density 
of  the  height  distribution  and  the  autocorrelation  function.  More  recently, 
WHITEHOUSE  and  PHILLIPS  [1978,1982]  have  generalized  the  previous  model  in 

order  to  incorporate  a  general  autocorrelation  function  C(A)  and  in  order 
to  have  a  fully  2D  surface  model  (rather  than  a  ID  profile  model). 

An  alternative  model  was  developed  by  NAYAK  [l97lj.  In  his  model, 
the  height  of  a  rough  surface  z  ■  z(x,y)  is  considered  to  be  a  two  dimen¬ 
sional  random  variable,  the  cartesian  coordinates  x  and  y  being  the 
independent  variables.  The  essential  differences  between  this  model  and 
those  of  WHITEHOUSE  et  al.  are  the  following:  this  model  is  based  on  con¬ 
tinuous  (rather  than  discrete)  properties  of  random  waveforms  and  uses  the 

power  spectral  density  function  (rather  than  the  autocorrelation  function). 
In  this  work  NAYAK  also  emphasized  the  distinction  between  the  statistics 
of  the  surface  and  the  statistics  of  a  profile  of  the  surface.  This  author 
then  shows  how  the  higher  order  surface  statistics  of  interest  can  be 
correctly  obtained  from  the  parameters  a  ,  a'  and  a"  (r.m.s.  values 
of  z  and  its  first  and  second  derivatives,  respectively)  of  a  single 
profile,  for  Gaussian  isotropic  surfaces. 

BUSH,  GIBSON  and  THOMAS  [l975]  developed  a  model  based  on  NAYAK's 
random  surface  theory  with  the  additional  assumption  that  the  cap  of 
each  asperity  is  replaced  by  a  paraboloid  having  the  same  height  and 
principal  curvatures  as  the  summit  of  the  asperity.  The  asperity  deforma¬ 
tion  was  assumed  to  be  elastic  (Hertzian).  The  authors  conclude  that 
for  large  separation  the  contact  area  is  proportional  to  the  load  and 
for  smaller  separation  the  proportionality  is  only  approximate.  The 
normal  stiffness  becomes  vanishingly  small  at  large  separations  and,  at 
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constant  separation ,  the  stiffness  is  proportional  to  the  load.  Also 
the  stiffness  is  inversely  proportional  to  the  square  of  the  r.m.s.  rough¬ 
ness.  The  explicit  relations  obtained,  for  large  Separation,  h  ,  are 
the  following: 


exp 


exp 
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dh 
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exp 


(3. A) 


where  S  is  the  normal  stiffness  per  unit  area  of  contact. 

O'CALLAGHAN  and  CAMERON  [1976],  using  a  model  with  some  features  in 
common  with  the  previous  one,  obtained  results  (also  for  elastic  contact) 
that  are  qualitatively  analogous  to  those  of  previous  authors:  approximate 
proportionality  between  real  area  of  contact  and  normal  load  and  also  ap¬ 
proximate  proportionality  between  separation  and  minus  the  logarithm  of  the 
load. 

Other  authors  have  also  modeled  conditions  on  the  contact  that  are 
more  complex  than  a  simple  elastic  or  plastic  deformation  of  the  asperities. 
MIKIC  [1974]  modeled  elasto-plastic  conditions.  He  assumed  that  the  real 
contact  pressure  is  equal  to  the  lower  hardness  of  the  contacting  materials 
and  that  the  total  deformation  of  each  asperity  is  the  sum  of  elastic  and 
plastic  deformations  produced  by  that  pressure.  Analogous  assumptions  were 
later  used  by  ISHIGAKI  et  al.  [1979]  to  model  loading-unloading  situations. 
Work-hardening  of  the  asperities  was  first  considered  by  HALLING  and  NURI 
[1975].  Work-hardening  and  the  mechanics  of  unloading-reloading  were  incor¬ 
porated  in  the  model  of  FRANCIS  [1977J. 
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3X2  The  Mode  of  Deformation-Plasticity  Indices 

While  many  elastic  and  elasto-plastic  models  of  surface  response  to 
normal  loads  have  been  proposed,  the  question  of  what  circumstances  dictate 
the  use  of  an  elastic  or  a  plastic  representation  remains  central.  In  par¬ 
ticular,  when  can  plastic  deformations  be  expected  to  be  sufficiently  large 
to  be  a  necessary  ingredient  in  a  reasonable  model  of  surface  behavior? 

In  an  attempt  to  resolve  this  question,  several  authors  have  proposed  para¬ 
meters  to  indicate  which  mode  of  surface  deformation  may  be  dominant  in 
various  situations. 

GREENWOOD  and  WILLIAMSON  [1966]  used  their  model  of  the  surface  topo¬ 
graphy  to  study  the  limits  of  the  elastic  deformation  of  rough  surfaces. 
Although  their  model  assumed  elastic  contact,  they  used  it  to  study  the 
early  stages  of  growth  of  plastic  contact  areas  with  increasing  load.  The 
criterion  for  the  initiation  of  plastic  flow  at  the  asperity  level  was  based 
on  the  analysis  of  ball  indentation  tests.  The  criterion  for  the  limit  of 
elastic  behavior  of  the  contacting  surfaces  was  based  on  the  ratio:  plastic 
area  of  contact/total  area  of  contact.  For  values  of  this  ratio  greater 
than  around  0.02,  the  surfaces  can  no  longer  be  considered  to  behave  elasti¬ 
cally. 

These  authors  found  that  a  nondimen sional  parameter  could  be  defined 
which  predicts  the  mode  of  deformation  of  the  surface  over  very  wide  ranges 
of  loads.  This  parameter  is  the  so-called  plasticity  index 


When  ijj  is  less  than  0.6  plastic  contact  could  be 
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caused  only  if  the  surfaces  were  forced  together  under  very  large  nominal 

pressures;  when  is  greater  than  1.0  plastic  flow  will  occur  even  at 

trivial  nominal  pressures.  The  load  has  thus  little  effect  on  the  mode 

of  deformation  of  the  surface.  The  essential  factors  are  the  material 

properties  (--)  and  the  surface  finish  (2|o^.  For  most  engineering 

materials  and  surface  finishes  \^»  >  1  .  We  can  conclude  that  during 

the  initial  contact  of  these  surfaces  the  deformation  will  be  mainly 

plastic  even  at  the  lightest  loads , 

The  repeated  loading-unloading-loading  of  the  metal  surfaces ,  as 

in  normal  sliding  or  in  metallographic  polishing ,  will  produce  changes 

in  the  shape  of  the  asperities  (decrease  of  ~  )  and  the  deformation 

R 

will  become  elastic.  It  turns  out  that  over  a  wide  range  of  experimental 
conditions  the  real  contact  pressure  on  the  asperities  becomes  of  the 
order  of  0. 1H  to  0.3H  (TABOR  [l975]).  However,  GREENWOOD  and  WILLIAMSON 
[l966]  observed  that  the  wear  debris  produced  in  sliding  can  reverse  this 
trend  and,  by  damaging  the  surface,  can  prevent  ip  from  reaching  the 
elastic  range. 

The  effect  on  the  plasticity  index  \{>  of  the  continued  sliding  of 

two  metal  surfaces,  under  conditions  which  prevented  severe  wear  damage 

of  the  surfaces,  is  shown  in  Figure  10,  which  is  due  to  GREENWOOD  and 

WILLIAMSON  [  1 9 6 6 J .  In  this  figure  the  decrease  of  \p  down  to  values 

which  correspond  to  an  elastic  normal  deformation  of  the  asperities  is 

easily  observed.  Several  other  such  indices  have  been  proposed  by  other 
authors,  and  we  mention  as  typical  those  in  the  works  of  WHITEHOUSE  and 

ARCHARD  [1966],  ONIONS  and  ARCHARD  [1973],  and  MIKIC  [1974]. 

A  simplified  analysis  of  TABOR  [1975a]  based  on  the  indices  proposed 
by  GREENWOOD  and  WILLIAMS  [1966]  and  WHITEHOUSE  and  ARCHARD  [1966],  showed 
that  a  critical  geometrical  factor  in  the  transition  from  elastic  to  plas 
tic  deformation  of  a  rough  surface  is  the  average  slope  of  the  asperities. 
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Figure  10.  Effect  of  continued  sliding  on  the  plasticity  index. 

The  behavior  of  three  mild  steel  specimens  slid  against 
a  stainless  steel  woven  grid  at  100  cm/s  under  loads  of 
1.5  kg.  (GREENWOOD  AND  WILLIAMSON  [1966]). 
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3.1.3  Some  Experimental  Results 

Early  experimental  results  on  the  stiffness  of  contacting  surfaces 
and  their  real  area  of  contact  can  be  found  in  KRAGELSKII  [1965].  This 
author  refers  to  the  work  of  SOKOLOVSKII  [1946],  BOBRIK  [1947]  and 
VOTINOV  [1940]  who  proposed  a  relationship  between  the  approach  (a)  and 
the  nominal  pressure  (p)  of  the  form: 

a  =  C  pm  +  K  (3.5) 

where  C,  m  and  K  are  constants  determined  experimentally.  According  to 

SOKOLOVSKII,  these  constants  are  usually  m=0.3-0.5  and  K=0.  Figure  11, 

also  from  the  same  book,  illustrates  the  experimental  relationship  between 

the  real  area  of  contact  (A  )  and  the  nominal  pressure  (p)  for  different 

r 

metals  in  contact  with  a  glass  prism.  It  can  be  seen  that  some  of  the 
curves  obtained  depart  appreciably  from  linearity.  According  to  KRAGELSKII, 
the  relationship  may  only  be  considered  linear  at  low  pressures  or  for 

materials  which  do  no  work  harden  appreciably. 

Important  sources  of  experimental  results  on  the  contact  of  metallic 
surfaces  can  be  found  in  the  literature  devoted  to  the  design  of  machine 
tools.  The  research  of  several  workers  who  studied  the  stiffness,  damping, 
friction  and  wear  characteristics  of  fixed  and  sliding  joints  has  been 
summarized  in  a  paper  by  BACK,  BURDEKIN  and  COWLEY  [ 1973] .  The  essential 
conclusions  related  with  the  stiffness  of  contav.  ing  surfaces  are  the 
following: 

(1)  For  low  nominal  pressures  ,  characteristic  of  sliding  connec¬ 
tions,  no  plastic  deformation  was  observed  and  the  approach-normal 
load  relationship  is  of  the  form  of  (3.5)  with  K=0  and  with  m 
having  values  in  the  range  m=0.32-  0.5,  in  perfect  agreement  with 
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Figure  11.  Real  area  of  eontaet  vs.  pressure  for  different  metals  in 
contact  with,  a  glass  prism. 

(1)  lead;  (2)  cadmium;  (3)  magnesium;  (4)  aluminum; 

(5)  copper;  (6)  steel  10 

Maximum  height  of  surface  asperities:  40  y.  (KRAGELSKII  [1965]) 


the  observations  of  SOKOLOVSKII  mentioned  above. 

While  these  authors’  reference  to  "no  plastic  deformation"  represents 
a  mild  contradiction  to  those  results  described  earlier,  a  study  of  avail¬ 
able  discussions  of  the  experimental  set-up  indicates  that  some  nib¬ 
bing  of  the  surfaces  may  have  been  done  prior  to  the  experiments  in  order 
to  obtain  reproducible  results  (e.g.  DOLBEY  and  BELL  [1971]). 

(2)  For  higher  nominal  pressures  ,  characteristic  of  fixed  or 
bolted  connections,  the  general  behavior  to  be  expected  when 
successive  loadings  and  unloadings  of  the  surfaces  are  performed, 
is  shown  in  Fig.  12.  It  is  clear  that  the  total  approach  of 
the  surface  has  an  elastic  and  a  plastic  component  and  that  (if 
the  joint  orientation  is  not  disturbed)  the  unloading  and  re¬ 
loading  paths,  are  practically  coincident  .  For  these  elastic 
paths,  CONNOLLY  and  THORNLEY  [1965]  proposed  the  following 
relationship  between  the  approach  and  the  normal  load: 

a  =  4"  In  P  +  Q  (3.6) 

D 


where  B  and  Q  are  constants. 

Tables  with  experimental  values  of  the  constants  C  and  m  in 
(3.5)  and  B  in  (3.6)  can  also  be  found  in  the  same  paper  of  BACK  et  al 
[1973]  together  with  the  corresponding  ranges  of  validity. 

The  elastic  stiffness  of  the  surface  per  unit  area  (S)  is  thus 
given  by: 
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da 
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(3.7) 


when  equation  (3.5)  is  used,  and  by: 


SOME  EXPERIMENTAL  POINTS  HAVE 
BEEN  OMITTED  FOR  CLARITY 
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(^UI/UOl) 


(CONNOLLY,  SCHOFIELD  and  THORNLEY  [1967]). 

The  experimental  values  a  contain  the  deformation  of  the  surface  asperities 
(penetrating  approach)  and  some  deformation  of  the  solid  material  in  the 
neighborhood  of  the  interface. 
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when  equation  (3.6)  is  used.  It  was  also  found  experimentally  that  for 
mild  steel  specimens  with  shaped  or  turned  surfaces  the  constant  B  is 
inversely  proportional  to  the  surface  roughness  in  c.l.a..  Consequently, 
the  stiffness  (3.8)  is  inversely  proportional  to  the  surface  roughness. 

In  recent  years,  much  experimental  work  has  been  done  to  test  the 
statistical  models  of  the  surface  geometry  and  contact.  The  essential 
conclusions  of  a  survey  paper  of  WOO  and  THOMAS  [l98o]  are  the  following: 

(i)  The  bearing  area  ratio  (A^/A)  increases  as  the  0.8  power  of 
the  dimensionless  load  (p/H)  up  to  values  of  p/H  of  the  order 
of  0.1  (see  Fig.  13). 

(ii)  The  dimensionless  separation  (h/a)  decreases  proportionally 
to  the  increase  of  the  logarithm  of  the  dimensionless  load  (p/H) 
(see  Fig.  14) . 

31.4  Conclusions 

Since  the  assumptions  on  the  roughness  of  the  surfaces  and  their 
mechanical  behavior  are  different,  the  relationships  involving  area  of 
contact,  load  and  separation  derived  from  the  models  described  earlier 
are  also  different  and,  in  some  cases,  difficult  to  compare.  However, 
comparing  the  theoretical  predictions  of  those  models  with  the  experi¬ 
mental  results  presented  above,  some  broad  conclusions  can  be  drawn. 

We  summarize  them  as  follows: 

On  the  Mode  of  Deformation  of  Metal  Surfaces 
(i)  The  essential  factors  affecting  the  mode  of  deformation  of  a 

rough  surface  are  the  material  properties  (E',  H)  and  the  surface 
finish.  The  normal  load  is  expected  to  have  little 


effect  on  the  mode  of  deformation  of  the  surface. 
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Figure  13.  Bearing  area  ratio  vs.  dimensionless  load. 

Least  squares  straight  line  obtained  on  the  basis 
of  the  experimental  results  of  6  different  papers 
by  various  authors. 

(WOO  and  THOMAS  [1980]). 
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( ii )  For  most  engineering  materials  and  surface  finishes  the  initial 
contact  of  the  surfaces  is  expected  to  be  plastic  even  at  light 
loads. 

( iii )  The  repeated  loading-unloading-reloading  of  the  metal  surfaces , 
as  in  normal  sliding  or  in  metallurgical  polishing ,  produces 
changes  in  the  shape  of  the  asperities,  which  lead  to  a  subse¬ 
quent  elastic  deformation,  provided  that  severe  wear  is  prevented 
during  the  process  of  sliding. 

On  the  Stiffness  of  Compressed  Rough  Surfaces 

( iv )  At  small  penetrating  approaches  (large  separations )  the  stiffness 
of  rough  surfaces  becomes  vanishing ly  small 

(v)  The  stiffness  of  a  surface  is  inversely  proportional  to  its 
roughness 

(vx)  The  normal  load  increases  roughly  as  cm  exponential  function  of 
the  penetrating  approach  (the  separation  decreases  proportion¬ 
ally  to  the  increase  of  the  logarithm  of  the  load). 

( vii )  For  light  loads,  because  of  (iv),  the  normal  load  is  closely 
proportional  to  a  power ,  in  the  range  ~~  to  •jJ—-  ,  of  the 
penetrating  approach. 

(viii)  For  the  same  normal  load  the  stiffness  during  the  first  loading 
is  smaller  than  the  stiffness  during  unloading  or  reloading , 
due  to  the  plastic  deformation  which  occurs  during  the  first 
loading. 

On  the  Real  Area  of  Contact 

(ix)  Most  of  the  statistical  models  developed  predict,  for  both  elastic 
and  plastic  contact,  an  exact  or  an  approximate  proportionality 
between  real  area  of  contact  and  normal  load. 
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(x)  Experimental  results  show  that,  often,  the  real  area  of  contact 

increases  with  a  power  of  the  normal  load  which  is  slightly  smaller 
than  one,  even  for  light  loads. 

A  definitive  precise  conclusion  seems  difficult  due  to  the  absence 
of  a  really  satisfactory  method  for  the  experimental  determination  of  the 
real  area  of  contact  —  a  fundamental  difficulty  pointed  out  by  TABOR 
[1981]. 

3.2  Dynamic  Contact  of  Rough  Surfaces 

In  section  3.1  we  presented  a  review  of  theoretical  models  and  experi¬ 
mental  results  on  the  static  contact  of  rough  surfaces.  Our  purpose  in 
this  section  is  to  determine  whether  or  not  the  conclusions  drawn  for  the 
static  deformation  of  rough  surfaces  can  be  expected  to  hold  for 
dynamic  situations.  We  are  only  interested  in  dynamic  situations  for  which 
no  bulk  plastic  deformation  of  the  metallic  bodies  occurs.  In  the  follow¬ 
ing  we  describe  experimental  observations  of  two  phenomena  for  which  the 
above  criterion  is  satisfied:  the  low  velocity  impact  of  metallic  bodies 
and  the  dynamic  response  (to  a  sinusoidal  excitation)  of  metallic  bodies 
compressed  together  with  some  preload. 

The  problem  of  the  impact  of  two  elastic  bodies  (two  spheres  or  a 
sphere  and  a  massive  body)  with  relative  velocity  before  impact  smaller 
than  the  critical  velocity  required  to  initiate  plastic  deformation, 
has  been  studied  by  several  authors  on  the  basis  of  the  HERTZ  theory 
of  impact  (see  LOVE  [l944]).  It  has  been  observed  (LIFSHITZ  and  KOLSKY 
[l963] )  that  the  predictions  of  this  theory  for  the  time  of  contact  of  a 
steel  sphere  with  a  polished  face  of  a  mild  steel  block  is  given,  with 
good  accuracy,  by  the  HERTZ  solution  (even  up  to  velocities  of  impact  that 
are  eight  times  larger  than  the  critical  velocity  corresponding  to  the 
onset  of  plastic  deformation  in  the  block).  On  the  other  hand,  the  HERTZ 
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theory  does  not  predict  any  energy  dissipation  and,  consequently,  for 
velocities  of  impact  smaller  than  the  critical  velocity  mentioned  above  the 
velocity  of  rebound  would  be  expected  to  equal  the  velocity  of  impact, 
i.e.,  the  coefficient  of  restitution 

e  -  V  /V  (3.9) 

r  o 

would  be  expected  to  equal  unity  (in  (3.9)  Vq  and  Vr  are  the  absolute 
values  of  the  velocities  of  impact  and  rebound  of  the  sphere,  respectively). 

However,  the  above  predictions  are  not  confirmed  experimentally.  Co¬ 
efficients  of  restitution  of  the  order  0.90  -  0.95  (energy  losses  of  the 
order  10  -  20  percent)  are  reported  by  HUNTER  [1957]  for  collisions  between 
ball  bearings  and  steel  blocks  and  coefficients  of  restitution  of  the  order 
0.985  (energy  loss  j  of  the  order  3  percent)  are  reported  by  TILLETT  [1954] 
for  steel  balls  on  glass  blocks.  These  discrepancies  between  the  theory  and 
the  experimental  observations  prompted  several  authors  to  study  the  possible 
sources  of  energy  dissipation.  HUNTER  [l957j  estimated  the  energy  lost 
by  the  impacting  ball  due  to  the  excitation  of  elastic  vibrations  on  the 
massive  body  subsequently  to  the  impact.  JEHCKEL  and  KLEIN  [l952]  , 
and  TILLET  [1954]  correlated  the  energy  lost  by  the  impacting 
ball  with  the  energy  dissipated  in  the  block  due  to  internal  friction. 

It  was  concluded  that  the  effects  considered  above,  together  with  surface 
friction,  air  damping,  etc.,  were  sufficient  to  explain  the  small  losses 
observed  in  the  impact  of  the  steel  balls  on  glass  blocks.  However  the 
same  was  not  true  for  the  case  of  steel  blocks.  LIFSHITZ  and  KOLSKY  [l963j 
then  investigated  the  possibility  of  those  losses  being  the  result  of 
local  plastic  deformation  of  the  steel  block  —  flattening  of  microscopic 
protuberances  (asperities).  To  investigate  this,  measurements  of  coeffi¬ 
cients  of  restitution  of  steel  balls  impacting  on  mild  steel  blocks  with 
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different  surface  finishes  were  carried  out.  Those  authors  found  that 
with  blocks  with  finely  ground  surface  finishes  the  coefficients  of 
restitution  were  very  scattered  and  never  exceeded  0.82;  for  surfaces 
polished  with  very  fine  emery  paper  the  results  were  much  less  scattered 
and  a  maximum  value  of  0.87  was  obtained;  for  a  highly  polished  surface 
(highly  reflecting  mirror  finish)  the  coefficient  of  restitution  was 
found  quite  consistent  at  0.95.  LIFSHITZ  and  KOLSKY  thus  conclude  that 
the  coefficient  of  restitution  is  very  sensitive  to  the  surface  finish 
and  that,  probably,  further  polishing  would  produce  higher  values  of  e  , 
closer  to  the  theoretical  predictions. 

We  note  that  these  observations  on  the  effect  of  the  surface  roughness 
on  the  rebound  of  metallic  bodies  can  be  easily  predicted,  at  least  quali¬ 
tatively,  from  the  models  discussed  in  Section  3.1.  We  recall  that  for 
most  engineering  surface  finishes  and  materials  the  plasticity  index  assumes 
values  which  correspond  to  a  plastic  deformation  of  the  asperities,  upon 
first  loading,  for  a  wide  range  of  normal  loads.  This  obviously  leads 
to  important  energy  dissipation  which  contributes,  in  the  manner  exper¬ 
imentally  observed,  for  the  low  values  of  coefficient  of  restitution. 

In  particular,  for  mild  steel  (the  material  of  the  blocks  used  by  LIFSHITZ 
and  KOLSKY)  the  model  of  GREENWOOD  and  WILLIAMSON  [l966j  predicts  that 
only  for  a  highly  polished  surface  and  small  loads,  can  the  deformation  of 
the  asperities  be  essentially  elastic  (see  Fig.  15  due  to  TABOR  {l975b]). 

We  can  thus  conclude  that  GREENWOOD  and  WILLIAMSON'S  model  allows  for  an 
interpretation  of  the  experimental  results  which  corroborates  the  conclu¬ 
sions  of  LIFSHITZ  and  KOLSKY. 

Experimental  results  of  more  quantitative  nature  on  the  dynamic 
characteristics  of  metallic  surfaces  were  provided  by  ANDREW,  COCKBURN 
and  WARING  [l967-68j.  In  their  experiments  the  annular  interfaces  of 


(Kg/mm2) 
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(<7*/R) 

Figure  15.  Nominal  pressure  at  which  the  transition  from  elastic  to  plastic 

deformation  occurs  for  a  flat  metal  of  specified  roughness  pressed 
on  to  a  flat  ideally  smooth  hard  surface. 

The  model  used  for  the  calculations  is  the  model  of 

GREENWOOD  and  WILLIAMSON  [1966]. 

Material  properties  used: 

_2 

Aluminum  H  =  40  Kg  mm 
Mild  Steel  H  =  120  Kg  mm"2 
Tool  Steel  H  =  400  Kg  mm"2 

-2 

Ball-bearing  Steel  H  =  900  Kg  mm 

For  each  material  below  the  "elastic  boundary"  the  deformation 
is  elastic.  The  onset  of  plastic  deformation  of  the  asperities 
occurs  on  the  "elastic  boundary".  (TABOR  [1975b]). 
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several  mild  steel  discs  compressed  together  with  some  preload  were  sub¬ 
jected  to  a  normal  harmonic  force.  From  their  dynamic  response,  the 
stiffness  and  damping  of  the  interfaces  was  evaluated  for  different  dimen¬ 
sions  of  the  apparent  area  of  contact,  surface  finishes,  and  for  various 
values  of  the  normal  preload.  Although  some  of  their  experimental  results 
(influence  of  apparent  area  of  contact  and  roughness  on  the  stiffness) 
cannot  completely  be  interpreted  from  the  static  contact  results,  other 
results  are  clearly  similar  to  those  observed  in  the  static  case.  The 
linear  dependence  of  the  stiffness  on  the  normal  load  and  the  (small) 
difference  between  the  initial  loading  and  unloading  stiffnesses ,  shown  in 
Fig.  16,  could  be  predicted  from  eq.  (3.8)  and  Fig.  12,  respectively, 
obtained  experimentally  for  static  contact. 

Although  the  experimental  evidence  available  is  not  abundant,  we 
think  that  the  facts  collected  in  this  section  suggest  clearly  that 
the  dynamic  behavior  of  most  rough  surfaces  is  consistent 
with  the  stiffness  characteristics  analyzed  in  detail 
in  the  static  case . 

4.  DRY  SLIDING  FRICTION  OF  METALLIC  BODIES 

4 . 1  Classic  Laws  of  Friction 

When  two  metallic  bodies  in  contact  are  subjected  to  applied  forces 
which  tend  to  produce  relative  sliding  motion,  friction  stresses  develop 
on  the  interface  that  tend  to  oppose  that  motion. 

In  the  following  only  the  resultants  of  the  stresses  on  the  contact 
surface  will  be  considered:  N  ,  the  compressive  normal  force  and  F  , 
the  friction  force.  The  metallic  bodies  are  considered  essentially  as 
rigid  bodies  with  a  well  defined  tangential  relative  velocity  u^,  . 


TOTAL  STIFFNESS  (Ibf./in.  x  I06) 
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According  to  MOORE  [l975],  the  classic  laws  of  friction,  as  they 
evolved  from  early  studies  in  the  past  centuries,  are  the  following: 

(i)  The  friction  force  (at  the  onset  of  sliding  and  during  sliding) 
is  proportional  to  the  normal  contact  force. 


(4.1) 


The  coefficient  of  proportionality,  y  ,  is  known  as  the 

coefficient  of  friction.  Often  two  values  of  y  are  quoted:  the 

coefficient  of  static  friction ,  y  ,  which  applies  to  the  onset  of 

s 

sliding  and  the  coefficient  of  kinetic  friction,  y^  ,  which  applies  during 
sliding  motion. 

(ii)  The  coefficient  of  friction  is  independent  of  the  apparent  area 
of  contact 

(iii)  The  static  coefficient  is  greater  than  the  kinetic  coefficient 

(iv)  The  coefficient  of  kinetic  friction  is  independent  of  the  sliding 
velocity . 

Another  important  characteristic  of  the  friction  force  is  the 
following  (RABINOWICZ  [l965]): 

(v)  When  tangential  motion  occurs,  the  friction  force  acts 

in  the  same  direction  of  the  relative  velocity  but  in  opposite 


sense. 


u„ 


F  =  -  y  N 


(4.2) 


The  first  two  laws,  usually  known  as  the  Amontons'  Laws  of  friction, 
are  generally  observed  to  hold  for  gross  motions  of  effectively  rigid  bodies, 
However,  we  notice  that  deviations  from  the  first  law  have  been  reported  ai. 
various  circumstances:  an  increase  of  the  coefficient  of  friction  for  lighs. 
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loads  (BOWDEN  and  TABOR  [l964]  )  ,  or  an  increase  of  friction  for  loads 
greater  than  a  fairly  well  defined  value  which  corresponds  to  the  breaking 
of  the  oxide  films  on  the  surface  (BOWDEN  and  TABOR  [l 964^ )  or  yet  a 
decrease  of  the  friction  coefficient  for  very  high  loads  when  the  true 
area  of  contact  approaches  the  magnitude  of  the  apparent  area  of  contact 
and  bulk  plastic  deformation  of  the  bodies  in  contact  occurs  (BAY  and 
WANHEIM  [1976]). 

The  third  law  derives  from  classical  experiments  (EULER  [l75oJ): 

if  a  body  rests  in  equilibrium  on  an  inclined  plane  and  the  inclination 

of  the  plane  is  slowly  increased  up  to  the  angle  at  which  the  sliding 

initiates  (this  angle  is  known  as  the  angle  of  static  friction  0  = 

s 

tan  *  y  )  it  is  often  observed  that,  when  the  motion  starts,  it  starts 
s 

abruptly  and  the  body  acquires  a  large  velocity  in  a  short  time.  If 
the  kinetic  friction  force  were  equal  to  the  maximum  static  friction 
force,  the  slippage  down  the  inclined  plane,  at  the  critical  angle  0  , 

would  be  infinitely  slow,  since  no  source  for  tangential  acceleration 
would  exist.  The  source  for  the  observed  acceleration  is  thus  the  differ¬ 
ence  between  the  maximum  static  friction  force  (y  N)  and  the  kinetic 

s 

friction  force  (y.N)  .  The  stick-slip  phenomenon  discussed  in  Section  5 

K. 

is  also  intimately  associated  with  this  law. 

The  fourth  law  is  now  known  to  be  invalid.  A  large  volume  of  exper¬ 
imental  data  and  empirical  formulas  for  the  variation  of  the  friction  co¬ 
efficient  with  sliding  velocity  can  be  found,  e.g.,  in  KRAGELSKII  {1965J 
and  MOORE  [1975].  However,  as  noted  by  RABINOWICZ  [1965],  for  many  pur¬ 
poses  in  which  only  limited  velocity  ranges  are  of  interest,  the  kinetic 
friction  coefficient  may  be  taken  to  be  a  constant  independent  of  the 
sliding  velocity. 
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Also  according  to  RABINOWICZ  [l96s],  the  fifth  law  has  been,  essen¬ 
tially,  confirmed  by  experiment:  for  surfaces  without  pronounced  directional 
properties ,  the  instantaneous  friction  force  may  fluctuate  by  a  degree 
or  so  from  its  assigned  direction,  changing  direction  continuously  and 
in  random  fashion  as  sliding  proceeds. 

4.2  The  Origins  of  Friction  —  A  Brief  Review 

In  this  section  we  present  a  summary  account  of  the  theories  that 
have  been  proposed  to  explain  the  origins  of  the  friction  force.  Detailed 
historical  accounts  and  thorough  discussions  of  these  theories  can  be 
found  in  BOWDEN  and  TABOR  [l964],  KRAGELSKII  [l965  ,  1978]  and  DOWSON  [l978]  . 
4.2*L  The  Early  Theories  of  Friction 

Many  of  the  early  theories  of  friction,  developed  during  the  eighteenth 
and  nineteenth  centuries,  attempted  to  explain  the  frictional  behavior  in 
terms  of  the  surface  roughness  and  the  interlocking  of  the  surface  asper¬ 
ities.  It  was  thought  that  the  sliding  of  two  contacting  bodies  involved 
the  riding  of  rigid  asperities  of  one  surface  over  the  other.  If  the 
average  asperity  angle  is  0  the  coefficient  of  friction  would  be 
equal  to  tan  0  and  consequently  it  would  be  independent  of  the  load 
or  the  size  of  the  contacting  surfaces.  This  provided  an  explanation 
for  the  Amontons'  laws.  The  assumption  that  the  asperities  on  one  surface 
could  traverse  the  gap  between  asperities  on  the  other  provided  an 
explanation  for  the  fact  that  often  the  kinetic  friction  is  smaller 
than  the  static  one.  According  to  TABOR  [1981],  the  main  weakness  of 
this  early  work  on  friction  (AMONTONS  [1699],  EULER  [1750],  COULOMB 
[1785])  is  that  basically  the  models  proposed  are  non-aissipative  and 
friction  is  certainly  a  dissipative  process  —  a  critique  which  was 
first  expressed  by  LESLIE  [1804], 
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4w2.2  The  Adhesion-Plowing  Theory 

A  substantially  different  explanation  for  the  friction  phenomena 
was  provided  by  BOWDEN  and  TABOR.  Their  theory  (reviewed  in  BOWDEN  and 
TABOR  [1950,1964],  TABOR  f  1972, 1975^  198 1] )  has  been  the  most  widely 
accepted  in  recent  decades  among  the  researchers  of  solid  contact 
phenomena.  The  inter  facial  friction  between  metallic  bodies  is  attri¬ 
buted  essentially  to  two  causes:  the  formation  and  shearing  of  metallic 
junctions  between  the  surface  asperities  and  the  plastic  deformation  of 
the  softer  surface  by  hard  asperities.  As  a  consequence,  the  friction 
coefficient  can  be  given  as  the  sum  of  two  components  resulting  from  each 
of  the  above  effects, 

y  =  y  +  (4.3) 

a  p 

where  y  results  from  the  adhesion  (welding)  and  y  results  from  the 
a  p 

plastic  deformation  (plowing). 

The  Adhesion  Component  of  Friction.  As  noted  earlier ,  BOWDEN  and 
TABOR  [l964]  assume  that  when  two  clean  metal  bodies  are  put  in  contact, 
plastic  flow  at  the  tips  of  the  asperities  and  local  welding  between 
opposing  asperities  occur.  The  true  area  of  contact  is  then  proportional 
to  the  normal  load: 

Ar  =  N/H  (4.4) 

Here  H  is  the  hardness  of  the  softer  of  the  contacting  materials. 

Neglecting  the  contribution  of  plowing,  the  friction  force  is  then 
equal  to  the  force  Fa  required  to  shear  the  junctions  formed  in  the 
manner  described  above.  If  t j  is  the  average  shear  strength  of  the 
junctions,  then: 

F  =  At.  =  (t./H)N 
a  r  J  J 
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or,  dividing  by  N, 

y  «  t./H 

^a  j 

Due  to  plastic  yielding  and  work-hardening,  the  interface  between 
welded  asperities  is  as  strong  or  probably  stronger  than  the  undeformed 
material  in  the  hinterland.  For  this  reason,  the  average  shear  strength 
of  the  junctions  (Tj)  can  be  assumed  to  be  equal  to  the  shear  strength 
of  the  softer  of  the  contacting  materials  (t)  .  Consequently: 

U  -  t/H  (4.5) 

a 

The  Amontons’s  laws  of  friction  are  then  verified:  the  friction 
force  is  proportional  to  the  normal  load  and  independent  of  the  apparent 
area  of  contact. 

For  most  materials,  t  is  of  the  order  of  0.2H  so  that,  for  this 
simplified  model: 

Ua  «  0.2 

However,  for  clean  metals  enormous  values  of  \i  may  be  obtained 
and  even  for  metals  in  air  u  may  be  of  the  order  unity. 

According  to  the  same  authors,  this  discrepancy  can  be  overcome  if 
the  plastic  junction  growth  due  to  combined  normal  and  tangential  loading 
of  the  asperities  is  considered.  The  yield  criterion  assumed  for  this 
situation  is  of  the  form 

2  2  „2  2 

p  +  a  s  =  H  =  a  x  (4.61 

*r  r 

where  pr  and  sr  are  the  normal  and  tangential  stresses  at  each 
asperity,  i.e.,  the  average  true  stresses  on  the  true  area  of  contact 
and  a  is  a  constant  of  the  order  10  . 

The  above  criterion  predicts  that,  when  an  asperity,  initially  in 
pure  normal  contact,  has  its  tip  plastically  deformed. 
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p  -  H 
*r 

When  a  small  tangential  stress  sr  is  applied  further  plastic 
flow  will  occur.  However,  although  the  newly  added  stress  is  tangential, 
the  normality  rule  dictates  that  the  initial  plastic  movement  will  be  a 
normal  approach  and,  because  of  the  shape  of  the  asperity,  it  will  lead 
to  an  increase  in  the  area  of  contact.  As  is  increased,  the  plastic 

displacements  will  become  more  nearly  tangential. 

If  Aq  =  N/H  is  the  initial  true  area  of  contact,  due  to  the  normal 

load  alone  it  can  be  shown  from  equation  (4.6)  that  the  equation  governing 
the  increase  of  the  true  contact  area  with  the  tangential  load  is 


1  +  a 


2 


(4.7) 


With  real  surfaces  in  air  this  process  of  junction  growth  is  termi¬ 
nated  due  to  the  presence  of  the  contaminant  layers  which  are  weaker 
than  the  bulk  material.  If  x^  =  K  x  (with  K  <  1)  is  the  critical  shear 
stress  of  these  interface  films,  the  junction  growth  terminates  when  the 
tangential  stress  on  the  interface  reaches  the  value  x^  ;  then  gross 
sliding  occurs.  From  equation  (4.6)  sliding  occurs  when: 

2  2  2 
Pr  +  a  Xi  =  a  x 


The  adhesion  component  of  the  coefficient  of  friction  is  then: 

1 


A  x* 
r  l 


a 


A  p  Y7Z-2  . . 

r*r  a  (K  -  1) 


(4.8) 
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Figure  17  correlates  the  tangential  force  coefficient  $  =  F  /N 

Si  Si 

with  the  growth  of  the  true  contact  area  A^/A^  for  different  values  of 
K.  It  can  be  seen  that: 

(i)  for  K  +  1  ji  +  «  ;  this  corresponds  to  perfectly  clean  surfaces 

of  very  ductile  metals,  where  junction  growth  goes  on  indefinitely 

(ii)  a  very  small  amount  of  weakening  of  the  interface  (say,  K  =  0.95) 
reduces  the  coefficient  of  friction  to  reasonable  values  (of  the 
order  unity) 

CLii)  For  a  very  weak  interface  (say,  K  <  0.2)  the  junction  growth  is 
very  small;  this  situation  corresponds  to  the  presence  of  a  good 
boundary  lubricant  or  a  thin  film  of  a  softer  metal;  for  this 
case  can  be  written  with  sufficient  accuracy  as 

Pa  -  K/VcT 

and  using  equation  (4.6): 

V  =  T./H  (4.9) 

a  i 

which  is  formally  analogous  to  equation  C4.5). 

We  wish  to  point  out  that  the  theory  described  above  not  only  gives 

an  explanation  for  the  Amontons '  laws  of  friction,  but  also  allows  for 
interpretations  of  the  other  classic  laws  (section  4.1).  Following 
arguments  of  RABINOWICZ  [l951,  1965J,  the  static  friction  is  often  greater 
than  the  kinetic  one  because  the  strength  of  the  junctions  would  increase 
with  the  time  of  stationary  contact  (we  will  discuss  this  point  in  sec- 
tion  5.1.);  the  weak  dependence  of  the  friction  force  on  the  sliding 
velocity  would  be  a  consequence  of  the  small  rate  dependence  of  the 
strengths  of  most  solids;  the  opposite  directions  of  friction  and  sliding 
velocity  would  be  a  consequence  of  the  plastic  deformation  of  the  isotropic 
material  on  the  contact. 
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We  also  point  out  that  all  the  theory  is  based  on  the  proportionality 
'N  (cf.  aq.  (4.4)).  We  recall,  from  Section  3.1.,  that  the  assumption 
of  plastic  deformation  of  the  asperities  is  not  essential  to  obtain  this 
proportionality.  Consequently,  arguments  similar  to  those 
presented  earlier  can  also  be  used  to  explain  Amonton's  laws  in  the  case 
of  elastic  deformation  of  the  asperities.  On  the  other  hand,  deviations 
from  the  proportionality  Ar<?<N  ,  would  lead  to  deviations  from  the 
Amontons '  laws . 

Finally,  we  note  that  in  recent  years  an  extensive  effort  has  been 
devoted  to  the  understanding  of  the  physical  nature  of  the  interfacial 

bonds  between  metallic  bodies  and  the  factors  that  may  influence  the 
existence,  strength  and  persistence  upon  unloading  of  these  bonds 
(mutual  solubility  of  the  metals,  ductility,  roughness,  surface  films, 
etc.).  Reviews  of  these  developments  can  be  found  in  TABOR  [l977,  1 98 l] 
and  MAUGIS  [l98l]. 

The  Plowing  Component  of  Friction.  The  plowing  component  of  the 
friction  can  be  estimated  using  the  simplified  model  of  a  hard  conical 
asperity  grooving  on  a  softer  surface,  as  shown  in  Fig.  18  .  Assuming 
that  the  plastic-yielding  metal  is  isotropic,  the  plowing  component  of 
the  coefficient  of  friction,  pp  ,  can  be  shown  to  be 

u  «*  —  tan  0  (4 . 10) 

HP  71 

Due  to  the  small  slopes  of  the  asperities,  this  expression  leads  to 
values  of  the  order  of  0.07  (for  0  *  5°)  to  0.14  (for  0  =  10°).  As  a 
consequence  of  these  small  values,  the  contribution  of  the  plowing  compo¬ 
nent  of  friction  will  be  negligible  unless  the  adhesion  is  small. 


55 


0N 


Figure  18.  Plowing  component  of  friction  for  a  cone  sliding  on 
a  softer  metal. 
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4.2.3  Recent  Theories 

Despite  the  wide  acceptance  of  the  theory  described  above  and  its 
effectiveness  in  explaining  the  basic  laws  of  friction,  several  serious 
criticisms  have  been  offered  and  some  alternative  theories  have  been 
proposed,  especially  in  recent  years. 

Some  of  the  criticisms  (BIKERMAN  [l976])  result  from  the  fact  that, 
in  the  air,  most  of  the  solids  are  covered  with  various  films  which  pre¬ 
vent  atomic  contact  between  the  bulk  materials  of  two  contacting  bodies. 

\ 

As  a  consequence,  the  force  required  to  pull  apart  two  contacting  bodies 
is,  in  daily  life,  zero  (if  gravity  is  compensated).  Furthermore,  accord 
ing  to  the  same  author,  the  formation  of  adsorbed  layers  in  the  air  is  so 
fast  that  all  the  effects  associated  with  sliding  are  comparatively 
slow.  From  these  facts  and  the  critical  review  of  several  experimental 
observations  on  adhesion  and  friction,  BIKERMAN  concludes  that  the 
importance  of  the  adhesion  in  frictional  phenomena  is  negligible,  as 
long  as  the  term  friction  is  used  in  the  sense  of  a  force  which,  for 
a  given  path,  is  reproducible  and  is  a  definite  function  of  the  normal 
load. 

Other  alternative  theories  have  been  developed  by  authors  who  are 
especially  concerned  with  the  evolution  of  the  friction  force  during 
prolonged  sliding  and  with  the  interaction  between  friction  and  the 
wear  damage  of  the  surfaces.  In  the  following  we  summarize  some  of  the 
more  elaborate  theories  that  have  been  proposed.  Our  purpose  is  only 
to  exemplify  the  type  of  ideas  which  have  been  advanced:  a  renewed 
interest  on  the  interlocking  effects  of  the  roughness ,  now  associated 
with  plastic  deformation  —  plowing  by  hard  asperities  and  entrapped 
wear  particles  and  deformation  of  the  asperities  and  of  the  subsurface 


layers. 
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RIGNEY  and  HIRTH  [l979j  developed  a  new  model  to  explain  the  origin 
of  friction  during  the  steady-state  sliding  of  metals.  This  steady-state 
sliding  is  the  advanced  stage  of  sliding  during  which  a  near  surface 
micro-structure  of  the  type  indicated  in  Fig.  19  (developed  during  the 
initial  break-in  stage)  remains  constant.  In  particular,  the  average 
thickness  of  the  highly  deformed  region  with  cell  structure,  remains 
constant  during  the  steady-state  sliding.  With  the  assumptions  that 
the  frictional  force  arises  from  plastic  deformation  in  the  near-surface 
region  and  that  most  of  the  deformation  work  is  confined  to  the  cell 
region,  those  authors  equate  the  work  of  the  friction  force  to  the  work 
of  deforming  plastically  the  subsurface  volume  and  obtain  an  expression 
for  the  friction  coefficient. 

KUHLMANN-WILSDORF  [l98l]  proposes  the  modification  of  the  adhesion 
theory  of  friction  to  include  the  interlocking  effect  due  to  the  micro¬ 
roughness  superimposed  on  the  hills  and  tops  of  the  asperities.  This 
micro-roughness  is  thought  to  be  due  to  inhomogeneities  in  the  near 
surface  deformations  and  its  importance  during  sliding  is  suggested  by 
its  persistence  during  the  sliding  process  and  by  the  observation  of 
micro-grooves  on  the  wear  tracks.  The  same  author  analyzes  the  importance 
of  the  adhesion  during  sliding  and  concludes  that  only  for  very  smooth 
surfaces  (average  slopes  <  0.9°)  the  adhesion  will  be  strong.  A  detailed 
analysis  of  the  plastic  deformations  in  the  subsurface  material,  using 
concepts  of  the  theory  of  dislocations,  is  then  presented  and  quantita¬ 
tive  expressions  for  friction  and  wear  are  proposed. 

SUH  and  SIN  [l98l]  offered  the  following  explanation  of  the 


genesis  of  friction: 
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Figure  19. 


Longitudinal  section  of  a  wear  specimen. 

The  curved  lines  indicate  strain  and  the  arrow  indicates 
the  sliding  direction  (RIGNEY  and  HIRTH  [1979]). 


Very 

Highly  Plastically  Bulk  Material 

Deformed  Deformed  Undeformed 
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" The  coefficient  of  friction  between  sliding  surfaces  is  due  to 
the  various  combined  effects  of  an  asperity  deformation  component  yd  , 
a  component  y^  from  plowing  by  wear  particles  and  hard  surface  asperi¬ 
ties  and  a  component  y^  from  adhesion  between  the  flat  surfaces .  The 
relative  contribution  of  these  components  depends  on  the  condition  of 
the  sliding  interface ,  which  is  affected  by  the  history  of  sliding,  the 
specific  materials  used,  the  surface  topography  and  the  environment. " 

The  asperity  deformation  coefficient  yd  is  expected  to  be  very  im¬ 
portant  at  the  onset  of  sliding  (static  coefficient  of  friction)  and  is  also 
expected  to  contribute  for  the  steady-state  sliding  if  new  asperities 
are  continuously  generated  as  a  consequence  of  wear;  yd  is  not  expected 
to  contribute  significatively  in  early  stages  of  sliding.  High  values 
for  the  plowing  component  y^  are  associated  with  two  identical  metals 
sliding  against  each  other  with  deep  penetration  by  wear  particles  while 
low  values  are  obtained  either  when  wear  particles  are  totally  absent 
from  the  interface  or  when  a  soft  surface  is  slid  against  a  hard  surface 
with  a  mirror  finish.  The  component  y^  is  expected  to  contribute  con¬ 
tinuously  during  the  sliding  process  and  its  contribution  is  expected  to 
increase  when,  due  to  wear,  the  number  of  particles  entrapped  between  the 
two  surfaces  increases.  The  low  values  for  the  adhesion  component  y  are 

di 

associated  with  the  sliding  of  well  lubricated  surfaces  and  high  values 
with  identical  metals  sliding  against  each  other  without  any  contaminants 
or  oxide  layer;  y  is  expected  to  be  zero  (or  negligible)  at  the  onset  of 

cl 

sliding  due  to  the  presence  of  contaminants.  With  the  deformation  of  the 
asperities  and  the  exposure  of  fresh  new  surfaces,  the  adhesion  is  expected 


to  increase. 
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The  history  dependence  of  the  friction  force  is  then  the  consequence 
of  the  evolution  of  the  variables  that,  respectively,  affect  the  com¬ 
ponents  y, ,  y  and  y  of  the  friction  force:  the  roughness,  expressed  by 
u  p  a 

the  average  slope  of  the  contacting  asperities,  0;  the  ratio  of  the  width 
of  asperity  penetration  to  the  diameter  of  the  particle,  w/2R;  and  the 
interface  shear  strength  factor  T  /t  which  is  the  ratio  of  the  interface 
shear  strength  to  the  shear  flow  strength  of  the  softer  contacting 
material: 

y  =  y(0,  co/2R,  t  /t)  (4.11) 

4.3  Some  Remarks  on  Friction  and  Plasticity  Theory 

Although  specific  details  may  differ  (signif icatively . . . ) ,  most 
current  theories  of  the  origin  of  friction  single  out  plastic  deformation 
on  the  contact  neighborhood  as  a  major  phenomenon  associated  with  the 
frictional  sliding  of  metallic  bodies.  On  the  other  hand,  the  model  of 
a  particle  sliding  with  friction  on  a  surface  is  often  presented  as  a 
prototype  of  plastic  behavior.  A  question  which  then  arises  is  why 
classical  plasticity  theory  has  not  been  applied  successfully  to  fric¬ 
tion  problems  (viewed  as  boundary  plasticity  problems). 

One  reason  for  this  has  been  pointed  out  by  several  authors  (e.g. 
DRUCKER  [1954],  MICHALOWSKI  and  MROZ  [1978],  SANCHEZ-PALENCIA  and  SUQUET 
[1983],  KLARBRING  [1984]):  The  sliding  rule  dictated  by  the  Coulomb's 
law  of  friction  is  not  an  associated  flow  rule  as  defined  in  the  Clas¬ 
sical  Theory  of  Plasticity.  In  particular,  a  pointwise  version  of 
Coulomb  friction  is  encompassed  in  the  relations  (see  PIRES  and  ODEN 
[1983]) 
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ISTI  i“la„l 


a  <  0 
n  — 


A  -  0  If  la  I  <  jjla  I 
~T  ~  1 ~T 1  1  n 1 

"~T  =  X?T  if  =  UlaJ 


A  >  0 

Here  u  is  the  tangential  sliding  velocity  at  a  point  on  the  contact 

surface,  a  is  the  normal  stress  (the  "contact  pressure"),  A  is  a  real 

number,  and  a  is  the  tangential  stress  vector  (the  frictional  stresses). 

[Thus,  0  corresponds  to  N  and  0_  to  F  in  the  classical  rigid-body  case], 
n  ~l  ~ 

Thus,  one  can  define  a  slip  surface,  such  as  that  shown  in  Fig.  20,  in 
the  spirit  of  a  yield  surface  in  plasticity,  such  that  contact  stress 
states  on  the  surface  correspond  to  states  at  which  sliding  occurs. 
Since-u  =Aa  on  this  surface,  classical  friction  does  not  satisfy  the 
classical  normality  rule.  And,  therefore,  most  of  the  standard  arguments 
of  classical  plasticity  do  not  have  a  straightforward  extension  to 
frictional  phenomena  .[unless  the  normal  forces  on  the  contact  boundary 
are  known  a  priori,  in  which  case  the  frictional  forces  are  essentially 
derivable  from  a  convex  dissipation  functional  (see  DUVAUT  and  LIONS 
[1972])]. 

We  observe  that  the  theory  of  BOWDEN  and  TABOR  [1964]  suggests  a 
way  to  analyze  friction  contact  problems  in  terms  of  the  classical 
(associated)  plasticity  theory.  The  yield  criterion  (4.6)  can  be  re¬ 
written  in  terms  of  the  nominal  normal  and  tangential  stresses  in  the 
following  way 

f(p,  s,  6)  =  p2  +  as2  -  62H2  =  0 
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where  6  =  A^/A  6  [0,1]  is  the  bearing  area  ratio.  In  the  above  equa¬ 
tion  6  plays  the  role  of  a  hardening  (softening)  parameter,  which  can 
make  the  yield  surface  shrink  *  the  point  p  =  0,  s  =  0  if  there  is 
loss  of  contact  (6  =  0) . 

Similar  ideas  were  also  proposed  by  MICHALOWSKI  and  MROZ  [1978] 
who  included  an  additional  hardening  parameter  to  describe  the  contact 
hardening  due  to  tangential  sliding.  However,  these  authors  assume  that 
the  elastic  part  of  the  contact  deformation  is  linear  which,  as  seen 
earlier,  does  not  apply  to  rough  metallic  surfaces. 

Independently  of  the  difficulties  that  might  be  encountered  in 
trying  to  use  the  theories  of  BOWDEN  and  TABOR  [1964]  or  of  MICHALOWSKI 
and  MROZ  [1978]  in  continuum  models  for  quantitative  purposes  (and 
BOWDEN  and  TABOR  [1964]  advise  readers  not  to  try  to  push  their 
theory  too  far) ,  it  seems  that  their  essential  applicability  would  be  to 
quasistatic  friction  problems  (which  we  called  Type  I  in  the  Introduc¬ 
tion. 

The  recent  theories  mentioned  in  the  end  of  the  previous  section 
are  not  of  much  use  for  us  here  either,  since  it  is  not  our  final  goal 
to  analyze  the  effects  of  prolonged  sliding  and  wear  on  the  friction 
phenomena . 

We  still  have  to  look  for  further  experimental  evidence  on  the 
transition  between  static  and  kinetic  friction  before  we  can  try  to 
develop  a  model  for  computational  mechanics. 


64 


5.  STATIC  AND  KINETIC  FRICTION.  STICK-SLIP  MOTION 

5.1  The  Initiation  of  Sliding 

5.1.1  Preliminary  micro-displacements 

The  classic  laws  of  friction  described  in  Section  4.1  dictate  that, 
at  the  onset  of  sliding,  the  friction  force  is  equal  to  the  (static) 
coefficient  of  friction  times  the  normal  force.  However,  this  law  is 
valid  only  in  a  macroscopic  sense,  since,  before  macroscopic  sliding 
is  observed,  some  preliminary  micro-displacements  occur  at  lower  force 
levels.  Several  experimentors  have  studied  the  evolution  of  these 
microdisplacements  when  the  bodies  are  subjected  to  increasing  tangen¬ 
tial  forces.  References  for  these  works  can  be  found  in  BOWDEN  and 
TABOR  11964]. 

Typical  experimental  curves  of  the  tangential  force  coefficient 

$  =  F/N  versus  the  tangential  displacement  uT  ,  obtained  by  COURTNEY- 

PRATT  and  EISNER  [l957]  ,  are  shown  in  Fig.  21.  COURTNEY-PRATT  and 

EISNER  [1957]  also  measured  the  electrical  conductance  of  the  contacting 

surfaces  during  their  experiments.  Neglecting  the  conductance  of  surface 

films,  the  increase  in  conductance  c/c  is  related  with  the  increase  in 
’  o 

the  true  area  of  contact  A  /A  according  to: 

r  o 


Using  this  expression,  a  comparison  of  the  experimental  results  with 
the  predictions  of  (4.7)  is  given  in  Fig.  23. 

An  important  point  made  clear  by  the  experiments  of  these  authors 
was  the  nature  of  the  preliminary  displacements.  As  shown  in  Figs.  21, 
22,  these  displacements  are  essentially  irreversible  (plastic) .  The 


(COURTNEY -PRATT  and  EISNER  [1957.]) 
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Figure  22.  The  irreversible  part  of  the  tangential  displacement. 

Platinum  surfaces.  Normal  load  N  «  920  g. 

o  Observed  displacements  u,^  whilst  the  tangential  force 
is  applied. 

p 

•  Residual  displacements  u  obtained  by  subtracting  from 
the  observed  displacements  u^,  the  elastic  displacements 
(deduced  from  the  experimental  results  of  Fig.  21). 

—  Theoretical  elastic  displacements  computed  assuming  that 
they  result  from  elastic  distortion  of  the  slider  according 
to  the  equation  u_  =  [(2-V)/4Gr]F  due  to  MINDLIN  il949]  where 
G  =  modulus  of  rigidity  and  V  =  Poisson's  ratio. 

(BOWDEN  and  TABOR  11964]  after  experimental  results  of 
COURTNEY-PRATT  and  EISNER  [1957]). 


Figure  23.  Tangential  force  coefficient  versus  real  area  of  contact. 

Platinum  surfaces  :  o  clean;  •  lubricated.  The  full  line 
is  the  theoretical  curve  plotted  according  to  eq.  (4.7) 
with  a  =  12.  (COURTNEY-PRATT  and  EISNER  [1957]). 
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origin  of  the  elastic  component  was  not  clear  for  the  authors  who  were 
not  able  to  separate  it  into  the  elastic  deformation  of  the  specimen  and 
that  of  the  apparatus.  BOWDEN  and  TABOR  [l964] ,  analyzing  the  same  exper¬ 
imental  results  of  COURTNEY-PRATT  and  EISNER,  attribute  that  elastic 
component  to  bulk  elastic  distortion  of  the  slider.  On  the  other  hand, 
the  fact  that  the  plastic  component  of  the  displacement  (u,j,)  satisfies 

p 

du^/dcj)  2  0  at  $  *  0  strongly  suggests  that  when  the  tangential  force 
begins  to  be  applied  the  plastic  deformation  of  the  asperities  is  normal 
(rather  than  tangential)  exactly  as  predicted  by  the  junction  growth 
theory  of  BOWDEN  and  TABOR. 

5.12  Creep  Sliding 

In  all  the  experiments  mentioned  above,  whether  with  hard  or  with 
soft  metals,  COURTNEY-PRATT  and  EISNER  observed  that  the  micro-sliding 
of  the  specimens  continued  for  long  periods  of  time  after  each  increment 
of  tangential  force  was  applied.  For  example,  with  the  soft  metal  indium, 
at  the  end  of  half  an  hour  the  rate  of  movement  could  be  2  x  10  ^  cms 
For  harder  metals,  after  a  few  minutes,  the  movement  was  imperceptible. 

BURWELL  and  RABINOWICZ  [l 9 5 3 J  also  observed  the  creep  sliding  on 
low  melting  point  metals  at  low  velocities  (see  Figs.  24  ,  25  ).  In  Fig.  25 
each  point  on  the  initial  increasing  branch  of  the  y  -  u^,  curve  corres¬ 
ponds  to  the  final  steady  state  sliding  (constant  velocity  uT)  attained 
by  a  slider  acted  upon  by  a  constant  tangential  force,  which  (divided 
by  the  normal  force)  corresponds  to  a  friction  coefficient  y  .  The 
same  authors  point  out  that,  for  harder  metals,  such  as  aluminium, 
copper,  steel,  etc.,  it  is  probable  that  an  initial  increasing  branch 

of  the  y  -  u  curve  also  exists.  The  difficulty  to  provide  additional 

“9 

experimental  evidence  of  this,  at  velocities  of  the  order  10  cm/s,  is 
obviously  extraordinary.  In  any  case,  since  those  velocities  are  so 


Figure  24.  Creep  Slzding. 

Steel  on  indium,  clean;  normal  load  N  =  30 
tangential  force  F  =  200  g  (constant) 

(a)  displacement  as  a  function  of  time 

(b)  sliding  velocity  as  a  function  of  time 
(BURWELL  and  RABINOWICZ  [1953]). 
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t(days) 

(b) 


Figure  24.  Creep  Sliding. 

Steel  on  indium,  clean;  normal  load  N  =  300  g; 
tangential  force  F  =  200  g  (constant) 

(a)  displacement  as  a  function  of  time 

(b)  sliding  velocity  as  a  function  of  time 
(BURWELL  and  RABINOWICZ  [1953]). 


1 


ClT  (cms-1) 


Figure  25.  Coefficient  of  friction  as  a  function  of  final  sliding 
velocity . 

Clean  steel  on  indium;  normal  load  N  =  300  g.  Each 
point  on  the  increasing  branch  of  the  y  -  u^,  curve 
corresponds  to  the  final  steady  state  sliding- obtained 
by  imposing  a  constant  tangential  force.  For  tangential 
forces  corresponding  to  y  >  hmax  no  steady  state  sliding 
was  obtained.  The  points  on  the  decreasing  branch  of  the 
y  -  u^,  curve  were  obtained,  using  conventional  experimental 
techniques  by  imposing  a  much  higher  sliding  velocity. 
(BURWELL  and  RABINOWICZ  [1953]). 
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small,  it  can  be  concluded,  following  BOWDEN  and  TABOR  [l964],  that 
the  frictional  behavior  of  ordinary  engineering  metals  at  room  temper¬ 
ature  is  reasonably  well  explained  in  terms  of  their  plastic  properties 
without  introducing  the  part  played  by  creep. 

5.L3  Rate  Dependence  of  the  Static  Coefficient  of  Friction 

A  more  important  time  dependent  effect  on  the  static  friction  has 
been  observed  by  several  authors  who  have  studied  the  stick-slip  motion 
(see  section  5.2). 

Consider  a  slider  resting  on  a  surface  with  no  macroscopic  sliding 

motion  relative  to  the  surface  and  a  tangential  force  F  ,  applied  to 

the  slider,  being  increased  at  constant  rate  <J)  *  F/N  until  gross 

sliding  occurs.  Under  these  conditions,  it  can  be  observed  that  the 

value  u  of  (})  =  F/N  at  which  the  macroscopic  sliding  occurs  depends 
s 

on  the  rate  <J>  of  increase  of  the  tangential  force.  As  shown  in  Fig.  26, 

y  decreases  with  the  increase  of  the  rate  of  application  of  the  tangen- 
s 

tial  force. 

The  interpretation  given  by  early  researchers  for  this  phenomenon 

was  essentially  the  same  used  to  explain  the  difference  between  static 

and  kinetic  friction  (recall  section  4.2).  The  strength  of  the  contact 

junctions  would  increase  with  the  time  of  stationary  contact.  Consequently, 

the  smaller  the  time  of  stationary  contact,  t  ,  the  smaller  the 

s 

coefficient  of  static  friction  would  be  (see  Fig.  26 ) .  Expressions 
proposed  by  several  authors  for  this  time  dependence  of  the  coefficient 
of  static  friction  can  be  found  in  a  survey  paper  by  RICHARDSON  and 
NOLLE  [1976]  and  are  reproduced  in  TABLE  1.  Sketches  of  the  correspond¬ 
ing  behaviors  are  shown  in  Fig.  27  . 

However,  the  experimental  work  of  several  authors  has  made  clear 
that  this  interpretation  was  not  correct. 


Figure  26.  Rate  dependence  of  the  static  coefficient  of  friction 
and  S2  are  the  points  at  which  gross  sliding 
initiates;  and  $  are  the  rates  of  increase  of 

the  tangential  force  coefficients  4>^  and  (J^J  t  ^  and 
t  2  are  the  times  of  stationary  contact;  y^  and  y^ 
are  the  static  coefficients  of  friction. 
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Figure  27 .  Coefficient  of  static  friction  vs.  time  of  stationary  contact; 
equations  from  the  literature  (see  TABLE  l) . 

(RICHARDSON  and  NOLLE  11976]) 


SIMKINS  [1967]  carried  out  experiments  similar  to  those  of  COURTNEY- 
PRATT  and  EISNER  [1957]  to  observe  the  micro-displacements  before  gross- 
sliding.  He  found  that  higher  rates  of  loading  inevitably  led  to  macro¬ 
scopic  sliding  at  lower  force  levels.  However,  in  other  experiments 
designed  to  assess  the  influence  of  the  time  of  stationary  contact  on 
the  value  of  the  static  coefficient  of  friction,  he  could  not  find  any 
correlation  between  the  time  of  stationary  contact  and  the  value  obtained 
for  that  coefficient. 

JOHANNES,  GREEN  and  BROCKLEY  [1973]  and  RICHARDSON  and  NOLLE  [1976] 
carried  out  experiments  in  such  a  manner  that  they  could  vary  independently 
the  rate  of  application  of  the  tangential  force  and  the  time  of  stationary 
contact.  This  was  achieved  by  interrupting  the  increase  of  the  tangential 
force  or  delaying  its  application  (see  Fig.  28).  RICHARDSON  et  al  found 
that,  if  no  delay  was  introduced  (Figure  28  .a)  the  best  fit  to  the  exper¬ 
imental  data  (expressing  p  as  a  function  of  t  ,  which  in  this  case 

s  s 

is  not  independent  of  $  )  was  obtained  with  equation  4  of  TABLE  1. 

However,  if  the  time  of  stationary  contact  was  made  independent  of  the 
rate  of  application  of  the  tangential  load  (Fig.  28.b,c)  JOHANNES  et  al 
and  RICHARDSON  et  al  found  that  the  governing  variable  was  the  rate 

of  increase  of  the  tangential  force  and  not  the  time  of  stationary  contact. 

The  dependence  of  p  on  4>  obtained  by  RICHARDSON  and  NOLLE  £  1 9 7 6 J  is 
shown  in  Fig.  29  .  It  can  be  seen  that,  for  sufficiently  s.aall  load  rates, 
the  coefficient  of  static  friction  is  constant  and  equal  to  a  value  which 
is  the  usually  quoted  coefficient  of  static  friction.  For  large  loading  rates 
the  coefficient  of  static  friction  tends  to  be  constant  and  equal  to  a 
value  which  is  usually  interpreted  as  the  coefficient  of  kinetic  friction, 
although  comparative  measurements  have  rarely  been  made.  As  a  consequence 


I 


Figure  28.  Dependent  (a)  and  -independent  (b.  e).  rates  of  loading 
and  times  of  stationary  contact. 

Since  the  load  rate  <j>  is  the  same  for  the  three  exper¬ 
iments  y  would  be  the  same  regardless  of  the  differences 
s 

in  the  times  of  stationary  contact. 


S  ,  S.  .  S  are  the  points  at  which  gross  sliding  initiates 
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<j>  =  f/N 

Figure  29.  Static  friction  vs,  loading  rate.  Logarithmic  $  scale  and 
static  coefficient  of  friction  normalized  by  dividing  by  U 

ss 

the  coefficient  of  static  friction  at  low  loading  rates. 
(RICHARDSON  and  NOLLE  [1976]). 
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of  these  observations,  those  authors  conclude  that  empirical  expressions  of 
the  type  shown  in  TABLE  1  should  be  recast  with  loading  rate  as  the  indepen¬ 
dent  variable. 

5.2  Stick-Slip  Motion 
5.2J  Introduction 

When  two  surfaces  in  contact  slide  relative  to  each  other  (Fig.  30), 
intermittent  vibration  of  the  relaxation  type  can  often  be  observed, 
especially  at  low  driving  velocities  (Fig.  31).  This  phenomenon,  usually 
known  as  stick-slip  motion  (after  BOWDEN  and  LEBEN  [l939j  )  is  highly 
undesirable  in  tables  of  machine  tools,  the  precise  positioning  of  which 
is  fundamental  for  the  accuracy  of  the  work  performed.  In  other  applica¬ 
tions,  violent  oscillations  may  lead  to  surface  damage  and  failure  of 
machine  components. 

Many  authors  have  studied  the  stick-slip  phenomenon.  We  shall  des¬ 
cribe  here  some  of  the  more  fundamental  experimental  observations  on  stick- 
slip  motion  and  the  most  common  interpretations  for  the  phenomenon.  For 
more  complete  lists  of  references  see,  e.g.,  BELL  and  BURDEKIN  [l 969— 7oJ 
and  ANTONIOU,  CAMERON  and  GENTLE  [l976]  . 

BLOK  [l94o]  ,  analyzing  experimental  results  of  BOWDEN  and  LEBEN 

[l939]  (obtained  with  metallic  bodies)  and  of  PAPENHUYZEN  [l938j 
(obtained  with  a  rubber  specimen  on  glass  and  road  materials)  concluded 
that  the  essential  condition  for  the  occurrence  of  stick-slip  motion  is 
a  decrease  of  the  frictional  force  with  increasing  sliding  speeds.  In 
the  same  paper,  BLOK  provided  also  the  first  systematic  study  of  fric¬ 
tional  vibrations  and  established  a  quantitative  criterion  for  their 


appearance. 


Figure  30.  Models  of  two  Cequivalentl  sliding  systems  which 
may  have  stick-slip  oscillations. 

K  =  linear  stiffness;  C  =  linear  damping;  M  =  mass 

of  the  slider;  N  =  normal  load;  =  driving  velocity; 


u  =  friction  coefficient 


Typical  traces  of  stick-slip  ; 
(a)  and  ( b )  in  Figure  32. 

=  displacement  of  the  slid 

UT  =  driving  velocity. 


In  experiments  carried  out,  using  either  especially  designed  appara- 
ti  or  slightly  modified  machine  tool  tables  and  slideways,  it  has  been 
observed  that  the  amplitude  of  the  stick-slip  motion  decreases  when: 

( a )  the  driving  velocity  UT  increases  (see  Figs.  32-34) 

(b)  the  damping  coefficient  C  increases  (see  Fig.  32) 

(c)  the  spring  stiffness  K  increases  (see  Fig.  33) 

(d)  the  mass  M  of  the  slider  decreases  (see  Fig.  34). 

It  has  also  been  observed  that  the  frequency  of  the  stick-slip  motion 
increases  with  the  increase  of  the  driving  velocity  and  that  the  maximum 
value  of  this  frequency  approaches  the  undamped  natural  frequency  of  the 
system  (Figs.  35,  36,  37)  although  in  some  cases  the  oscillation  stops  at 
a  level  below  that  natural  frequency. 

The  frictional  properties  of  the  contacting  surfaces  and  the  dynamic 
characteristic  of  the  bodies  involved  (mass,  stiffness,  damping)  are  thus 
required  to  study  the  stick-slip  motion. 

The  equation  that  governs  the  tangential  movement  of  the  slider  in 
the  model  of  Fig.  30(b)  is  the  following  (note  that  appropriate  choice  of 
variables  leads  to  analogous  equation  for  the  model  of  Fig.  30(a)): 

M  iiT  +  C  uT  +  K  uT  =  F  (5.1) 

where  F  is  the  instantaneous  friction  force. 

The  difficulties  obviously  arise  from  the  friction  force  F  .  The 
essential  problems  have  not  been  mathematical,  they  have  been  mechanical 
and  experimental.  The  question  is:  how  does  the  friction  force  vary 
when  the  sliding  velocity  is  zero  or  in  the  neighborhood  of  zero? 
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Figure  32.  Relation  between  amplitude  of  the  stick-slip  motion  and 
velocity  for  various  damping  coefficients. 

THEORETICAL:  Model  (I)  &  (lie)  of  Fig.  38 

SYSTEM  (1)  C  =  0.01  lb.  sec.  in"1 

-1 

SYSTEM  (2)  C  =  0.32  lb.  sec.  in 
SYSTEM  (3)  C  =  0.59  lb.  sec.  in"1 


(BROCKLEY,  CAMERON  and  POTTER  [1967]) 
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Figure  33.  Relation  between  amplitude  of  the  stick-slip  motion 

and  driving  velocity  for  various  spring  stiffnesses  (K). 

Steel  on  steel  unlubricated 

Normal  load  N  =  1750  g 

(RABINOWICZ  [1965]). 
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Figure  3'*  .  Relation  between  amplitude  of  the  stick-slip  motion  and 
driving  velocity  for  various  masses  of  the  slider. 

K  =  86.9  Kg/nua 

-  Lubricant  A 

-  Lubricant  C 

•  M  =  53.2  Kg 

o  36.0  Kg 

®  11.3  Kg 


(KATO  and  MATSUBAYASHI  [1970J ) . 
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Figure  35.  Variation  of  stick-slip  frequency  with  driving  velocity 
for  various  normal  loads. 

Natural  frequencies  of  the  system 
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normal  mode  ~  200  Hz 
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(DOKOS  [1946]) 
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Figure  36.  Variation  of  stick-slip  frequency  with  driving  velocity 
for  various  lubrication  conditions. 

Natural  frequencies  of  the  systems  (tangential  motion) 
at  11.6  Hz 

Lubrication  conditions 
•  Unlubricated 
x  Lubricant  1 
o  Lubricant  2 

(BELL  and  BURDEKIN  [1969-79]) 


0T(in.sec"1) 

Figure  37.  Variation  of  stiek-slip  frequency  with  driving  velocity 
for  various  lubrication  conditions. 

Natural  frequencies  of  the  systems  (tangential  motion) 
at  25  Hz. 

Lubrication  conditions 
•  Unlubricated 
x  Lubricant  1 
o  Lubricant  2 

(BELL  and  BURDEKIN  [1969-70]) 
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5.2.2  The  Assumed  Friction  Characteristics 

Models  with  different  assumptions  on  the  dependence  of  the  friction 
on  time,  displacement,  velocity  and  acceleration  have  been  used  to 
characterize  the  occurrence  of  stick-slip  motion.  We  summarize  some  of 
the  most  representative  in  Fig.  38. 

These  models,  with  the  input  of  appropriate  data  will  simulate  stick- 
slip  motion.  The  accuracy  of  their  predictions  (e.g.,  the  amplitude  and 
frequency  of  the  oscillation  and  the  critical  driving  velocity  above 

which  no  oscillation  is  observed)  will  depend  on  how  close  the  assumed 
characteristic  for  the  friction  force  is  to  the  actual  characteristic 
existing  under  the  particular  conditions  involved.  Another  difficulty  is 
the  determination  of  unique  data  to  incorporate  in  the  above  models. 

5.2.3  The  Experimental  Friction  Characteristics 

The  limitations  of  the  above  models  can  be  better  appreciated  by 
analyzing  the  characteristics  obtained  experimentally  by  several  authors, 
as  shown  in  Figs.  39  -  4  3.  We  wish  to  point  out  the  following: 

(i)  In  obtaining  these  characteristics  some  authors  (e.g.  SAMPSON  et  al 
[l943j  )  have  smoothed  and  numerically  differentiated  the  data 
available:  traces  of  the  displacement  oscillations.  Some  irregular¬ 
ities  inherent  to  friction  processes  plus  effects  assumed  of  secondary 
importance  by  the  authors,  may  have  been  partially  obscured  by  using 
these  techniques. 

(ii)  Other  authors  (BELL  et  al  [l969-70}  ,  KATO  et  al  [l97oJ  )  have 
measured  simultaneously  displacements,  velocities  and  accelerations 
using  different  transducers.  As  pointed  out  by  ANTONIOU  et  al  [1976] 
since  each  transducer  has  its  own  response,  which  does  not  coincide 
absolutely  with  the  response  of  the  other  two  transducers,  the 
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Figure  38.  Models  of  the  variation  of  the  friction  force. 
(Ila)  BLOK  [1940] 

(I)  &  (Ila)  RABINOWICZ  [1957] 

(I)  &  (lib)  DERJAGUIN  et  al  [1957] 

(I)  &  (lie)  BROCKLEY  et  al  [1967] 

(lid)  BELL  et  al  [1969-70] 

(He)  BELL  et  al  [1969-70] 

(llf)  BANERJEE  [1968] 

(llg)  BO  and  PAVELESCU  [1982] 


0.002 


0.004 


Distance  slid  (mm) 

Figure  39.  Relation  between  coefficient  of  friction  and  distance  slid 
in  an  impact  experiment. 

Clean  copper  on  mild  steel.  The  sliding  of  a  body  on  an 
inclined  plane  is  initiated  with  the  impact  of  another  body. 
(&ABINOWICZ  [1951]). 


Figure  40.  Relation  between  coefficient  of  friction  and 
distance  slid  in  a  stick-slip  experiment. 

(KATO,  YAMAGUCHI  and  MATSUBAYASHI  11970]). 
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(b) 


0  8  16  24  32 


sliding  velocity  (cms-1) 

Figure  41.  Coefficient  of  friction  vs.  sliding  velocity . 

(a)  steel  on  Babbitt 

(b)  Wood's  metal  on  Babbitt  (1-first  slip;  2-second  slip) 

(c)  steel  on  steel 

(SAMPSON,  MORGAN,  REED  and  MUSKAT  [1943]). 


results  obtained  may  contain  the  effects  of  calibration  and  hysteresis 
errors . 

(iii)  Most  of  the  authors  have  studied  only  the  tangential  movement  of  the 
slider.  Consequently,  they  have  assumed  the  normal  load  on  the  con¬ 
tact  to  be  constant  and  equal  to  the  weight  of  the  slider  or  to  the 
force  which  produces  the  static  deflection  of  their  elastic  normal 
loading  system.  Explicitly  or  implicitly,  all  variations  observed 
on  the  friction  force  are  thus  assumed  to  be  the  conseauence  of  corres¬ 
ponding  variations  in  the  friction  coefficient.  We  note  that  BLOK, 
in  his  early  work,  rightly  emphasized  that,  in  general,  curves  showing 
a  decrease  of  frictional  force  with  velocity  do  not  implv  a  correspond¬ 
ing  decrease  of  the  coefficient  of  friction,  as  fluctuations  of  load 
might  also  occur. 

Some  comments  and  conclusions  suggested  by  Figs.  38  -  42  are  presented 

in  the  following: 

(a)  The  change  of  friction  force  from  the  maximum  static  value  to  the 

minimum  kinetic  value  occurs  in  a  very  short  distance  —  of  the 
_2 

order  of  10  mm  in  the  impact  experiments  of  RABINOWICZ,  and  of 
the  order  of  10  '  mm  in  the  stick-slip  experiments  of  KATO  et  al. 

(Figs.  39,  40). 

(b)  The  existence  of  a  second  slip  in  the  experiments  of  SAMPSON  et  al. 

(Fig.  41b)  shows  again  the  rate  dependence  of  the  maximum  static 
friction  force:  in  the  second  slip  (large  load  rate,  short  time  of 
stationary  contact)  the  friction  remains  at  a  level  close  to  the 
kinetic  friction. 

(c)  The  friction  force  is  not  a  reversible  function  of  the  sliding 
velocity.  The  accelerating  and  decelerating  branches  of  the  curves 
friction-velocity  are  usually  distinct  during  stick-slip  motion 
(Figs.  41  -  43). 
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(d)  Different  pairs  of  metals  in  contact  and  different  lubrication  con- 

I 

ditions  may  produce  characteristics  with  shapes  clearly  distinct 
(Fig.  41);  according  to  SAMPSON  et  al.  interchange  of  the  metal  on 
the  slider  and  on  the  driving  surface  often  makes  a  great  difference. 

(e)  Even  for  the  same  combination  of  metals  and  lubrication ,  the  shape 
and  slope  of  the  characteristic  and  the  separation  between  the  two 
branches  will  depend  on  the  dynamic  properties  of  the  system  (mass, 
stiffness )  and  on  the  driving  velocity  (see  Fig.  42).  This  means 
that  the  experimental  friction-sliding  velocity  curve  is  not  defined 
uniquely  by  the  nature  of  the  surfaces  in  contact  —  it  is  a  conse¬ 
quence  of  all  the  dynamic  variables  involved. 

From  the  above  observations,  mainly  from  (e),  we  conclude  that  an 
acceptable  theory  of  stick-slip  motion  will  have  to  explain  the  com¬ 
plex  relationship  between  the  friction  force  and  the  sliding  velocity, 
rather  than  assume  a  simplified  relationship  from  which  the  experimental 
evidence  will  deviate  often. 

5.3  Memory-dependent  Friction 

Several  authors  have  attempted  to  interpret  the  stick-slip  oscilla¬ 
tions  on  the  basis  of  friction  laws  which  allow  for  a  dependence  of  the 
friction  force  on  the  previous  sliding  history.  Early  studies  with 
metals  suggesting  friction  laws  of  this  type  are  due  to  RABINOWICZ 
[1958,  1959];  more  recent  studies  with  rocks  are  due  to  RUINA  [1980, 

1983],  RICE  and  RUINA  [1983],  and  GU  et  al.  [1983]. 

RUINA  [1983]  assumes  that  the  frictional  surface  has  at  any 
instant  in  time,  a  state  which  can  be  characterized  by  a  collection 
of  variables  ( internal  variables).  These  variables  represent  the 


surface  memory  of  previous  sliding  —  in  general  they  are  weighted 
averages  of  some  function  of  the  recent  sliding  velocity.  If  only 
constant  normal  stress  histories  are  considered  and  the  friction  stress 
is  assumed  to  be  proportional  to  the  normal  stress,  then: 

i)  the  instantaneous  coefficient  of  frictron  is  assumed  to  be  a 
function  of  the  sliding  velocity  and  the  state; 
ii)  the  instantaneous  rate  of  change  of  state  is  assumed  to  depend 
also  on  the  sliding  velocity  and  on  the  state. 

These  assumptions  can  be  expressed  in  the  form: 


t  =  p  f(u„  , a ) 
<*i  -  8i^x»  a  ) 


(  CL  ~  C  Ot  j  ,  C*2  ,  •  .  •  ]  ) 


For  each  constant  sliding  velocity  the  functions  g^  are  such 

s  s  * 

that  the  state  variables  evolve  towards  steady  state  values  (uT) 

•  g  g 

satisfying  g^(uT  ,  a  )  =  0  i  =  1,  2,  ... 

Several  examples  of  state  variables  friction  laws  were  provided  by 
the  same  authors.  The  simplest  one,  involving  only  one  state  variable, 
is  of  the  form: 


x  =  pLy  +  a  +  A  X,n(uT/uT)J 

a  »  -  rp  [a  +  B  £n(uT/uT)] 

c 


(52) 


•  is  0  is 

where  u  is  an  arbitrary  positive  constant  and  y  ,  A  ,  B  and  d£ 

are  positive  empirical  constants .  According  to  this  law,  the  frictional 

stress  has  a  positive  instantaneous  rate  dependence  expressed  by  A  £n 


(u^/u^)  and  the  state  variable  a  evolves  towards  the  steady  state 


ss 


*  -B  &n(uT/uT)  with  a  characteristic  slip  distance  d£ 


value  a 


A  dynamic  stability  analysis  of  the  sliding  with  constant  driving 


speed  ,  leads  to  the  following  conclusion  (for  single  state  variable 

friction  laws):  for  values  of  the  spring  stiffness  K  (Fig.  30. a) 

greater  than  a  critical  value  K  the  steady  motion  is  stable;  for 

values  of  K  smaller  than  K  the  steady  motion  is  unstable.  The 

cr 

value  of  K  is  found  to  be: 
cr 


•  q  C  •  * 

uTdx  (uT)/duT 


dc3T^T*  a)/3ux 


evaluated  at  u^,  =  .  A  necessary  condition  for  instability  is 


then: 


dTSS(uT) 


UT  =  UT 


For  the  simple  law  (5.2),  K  reduces  to 

cr 


(B-A)] 

d 


d  p  A 
cr 


This  expression  predicts  the  following  experimental  observation:  an 
increase  of  the  stiffness  or  a  decrease  of  the  mass  increases  the 
stability  of  the  steady  motion. 

However,  this  expression  seems  to  contradict  a  known  experimental 
fact  (see  the  results  of  RABINOWICZ  in  Fig.  33):  the  larger  the 
driving  velocity,  the  smaller  the  minimum  stiffness  required  for  a 


smooth  sliding. 


Finally,  we  note  that  a  detailed  study  on  the  applicability  of 
state  variable  friction  laws  to  metal  surfaces  has  not  yet  been  done. 
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5. 4  The  Importance  of  the  Normal  Degree  of  Freedom  in  Sliding 
Friction 

Here  we  follow  essentially  the  fundamental  paper  of  TOLSTOI  [l967] . 
Other  relevant  contributions  will  also  be  mentioned. 


It  was  seen  in  Sections  5.2  that  the  stick-slip  motion  is  a 
consequence  of  the  decrease  of  the  friction  force  with  sliding  velocity. 


i.e. 


J 


(5.2) 


As  in  section  4.2,  let  F  be  equal  to  the  product  of  the  real 
area  of  contact  by  the  average  shear  strength  of  the  material  being 
sheared, 


F  -  A  t 
r 

As  in  section  3.1.,  let  A^  be  a  (differentiable)  function  of  the  separa 
tion  between  the  contacting  bodies, 


A  =  A  (h) 
r  r 

Finally,  assume  that  both  x  and  h  are  differentiable  functions  of 
the  tangential  velocity 

T  =  x(uT) 
h  «  h(uT) 


Then, 

dF  _  ^Ar  dh  .  dx 
du^  "  T  dh  du^  r  du^ 


(5.3) 


We  shall  analyze  this  relation  in  hopes  of  establishing  conditions  for 
the  occurence  of  stick-slip  motion. 
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The  second  tern  in  the  second  member  of  (5.3)  is 

always  positive  (or  zero)  at  a  constant  temperature,  because  it  represents 
the  rate  of  increase  of  tha  rheological  shear  strength  with  the  rate  of 
shear.  It  is  also  clear  that  dF/dh  =  TdA^/dh  ,  in  the  first  term  of 
the  second  member,  is  always  negative  since,  as  h  is  increased,  the 
true  area  of  contact  decreases.  Furthermore,  TOLSTOI  confirmed  experi¬ 
mentally  tnat  dF/dh  <0  in  static  experiments ,  the  results  of  which  are 
shown  in  Fig.  44*  In  this  figure  it  can  be  seen  that,  for  the  particular 
surfaces  used  in  the  experiment,  friction  is  halved  as  the  slider  rises 


by  only  about  0.2  ym  .  The  friction  force  is  thus  very  sensitive  to  the 
normal  separation. 

In  the  same  paper,  TOLSTOI  presents,  among  others,  the  following 
argument  to  show  that  dh/duT  >  0  :  an  increase  in  the  speed  of  the 
slider  increases  the  upward  vertical  components  of  the  impulses  exerted 
on  the  slider  asperities  as  they  collide  with  those  of  the  underlying 
surfaces;  this  increases  the  amplitude  of  the  normal  natural  vibrations 
of  the  slider,  the  frequency  of  which  depends  on  the  contact  stiffness 
and  mass  of  the  slider;  due  to  the  nonlinearity  of  the  normal  force- 
separation  relationship  (recall  section  3.1  and  Fig.  44),  the  normal 
vibrations  of  the  slider  are  highly  asymmetric;  consequently,  an  increase 
of  the  amplitude  of  those  normal  vibrations  raises  the  mean  level  on  which 
the  slider  moves.  We  note  that,  for  steady  state  sliding  of  lubricated 
surfaces,  BELL  and  BURDEKIN  [l 969— 7oJ  obtained  experimental  results 
which  confirm  that  dh/du^>  0  under  those  conditons. 

Having  shown  that  dF/dh  =  t  dAr/dh  <0,dh/du,j,  >  0  and  A^  dx/du^  0  , 
it  follows  that  condition  (5.2)  for  occurrence  of  stick-slip  motion  may 


44.  Normal  load  (N,  curve  1)  and  static  friction  force  (F,  curve  2) 
versus  the  separation  h  (arbitrarily  fi  is  taken  to  te  zero  for 
the  maximum  normal  load  used  in  the  experiments) . 

Dry  steel  surfaces.  (TOLSTOI  [1967]). 
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not  be  always  fulfilled.  For  example,  whenever  h  is  kept  constant, 
dF/du.j,  is  positive  and  the  sliding  must  be  smooth.  To  illustrate  this 
point,  TOLSTOI  shows  two  ways  to  keep  h  constant  and  obtaining  a  smooth 
sliding: 

(1)  By  using  an  extremely  small  driving  velocity 

(2)  By  using  an  external  damping  of  the  normal  oscillations. 

The  author  presents  experimental  results  corresponding  to  both 
the  above  situations  (.see  Fig.  45). 

In  case  1,  a  lower  bound  for  the  extremely  small  velocity  required 
is  estimated  by  the  author  in  terms  of  the  asperity  height  and  spacing, 
the  coefficient  of  friction,  the  normal  yield  stress  and  the  creep  vis- 
cisity.  In  fact,  the  small  velocity  obtained  is  that  required  for  the 
creep  deformation  of  the  asperities  without  normal  displacements  (recall 
the  information  on  Section  5.1.2  on  creep  sliding). 

In  case  2  ,  the  effect  of  the  normal  damping  is  such  that  no  reduc¬ 
tion  on  the  friction  force  is  observed  when  the  sliding  velocity  is 
increased.  There  exists  no  quantitative  difference  between  kinetic  and 
static  friction.  Furthermore,  if  the  normal  damping  is  introduced,  with 
no  change  on  the  driving  velocity  during  a  run  which  shows  stick-slip 
oscillations,  it  is  observed  that  these  oscillations  disappear  and  that 
the  value  of  the  coefficient  of  friction  for  the  subsequent  smooth  sliding 
is  even  greater  than  the  maximum  values  obtained  at  the  end  of  the  stick 
periods  of  the  stick-slip  motion. 

The  responsibility  of  the  freedom  of  normal  displacement  for  both 
the  falling  F-u^  relation  and  the  stick-slip  motion  was  corroborated 
in  two  ways : 

(a)  Sufficiently  heavy  external  damping  of  the  tangential  vibrations  alone 
could  suppress  these  vibrations  but  failed  to  affect  the  negative 
slope  of  the  F-u.^  curve. 


I 


0.4 


o — O  —  o 


Sliding  velocity  (cms  ) 

Figure  45.  Coefficient  of  friction  as  a  function  of  s tiding 

velocity  under  external  normal  damping  (curves  1-8) 
or  without  normal  damping  (curves  9-10). 

1,  2,  5,  6,  7  -  unlubricated  steel 

3,  4,  9,  10  -  lubricated  steel 

-  smooth  sliding 

-  stick-slip  oscillations  (the  points  plotted 

correspond  to  mean  values  of  the  friction 
coefficient  during  the  oscillations). 

Note  that  only  when  no  external  normal  damping  is 
imposed,  the  stick-slip  oscillations  do  appear  (curve  lO'1 

(TOLSTOI  [1967]). 
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(b.)  Direct  observations  of  the  normal  displacements  of  the  slider  carried 
out  at  low  driving  speeds  showed  that  the  forward  movements  of  the 
undamped  slider  invariably  ocurred  in  strict  synchronism  with  its 
upward  jumps  (see  Fig.  46). 

We  note  that  several  authors  have  also  done  observations  analogous 
to  (b).  In  their  early  study  BOWDEN  and  TABOR  [l939_]  measured  the  elec¬ 
trical  conductance  during  stick-slip  motion  and  observed  a  marked  fall  of 
electrical  conductance  during  the.  slip  phase  of  the  motion.  This  sug¬ 
gests  the  occurrence  of  normal  jumps  of  the  slider  during  the  slip  phase. 
Oscillograms  of  the  normal  contact  oscillations  which  occur  in  stages  of 
tangential  jumps  of  the  slider  during  stick-slip  motion  were  obtained  by 
TOLSTOI,  BORISOVA  and  GRIGOROVA  [1971].  Study  of  those  oscillations  re¬ 
vealed  a  fundamental  frequency  consistent  with  the  normal  interface 
stiffness  properties  (see  BUDANOV,  KUDINOV  and  TOLSTOI  [1981]).  Direct 
measurements  of  the  separation  of  unlubricated  and  lubricated  surfaces 
during  stick-slip  motion  were  also  done  by  BO  and  PAVELESCU  [1982]  and 
by  TUDOR  and  BO  [1982],  respectively.  Other  experimental  evidence  of 
the  influence  of  the  normal  oscillations  on  the  stick-slip  motion  was 
presented  by  ANT0NI0U,  CAMERON  and  GENTLE  [1976]. 

In  summary,  TOLSTOI  [l 96 7j  emphasized  the  fundamental  role  played 
by  the  normal  stiffness  of  the  contacting  asperities  in  sliding  friction. 
In  particular,  he  concludes  that  the  decreasing  friction-velocity  charac¬ 
teristic  at  small  velocities  is  not  an  intrinsic  property  of  sliding 
friction:  it  is  solely  the  consequence  of  the  microvibrations  in  the 
direction  normal  to  the  sliding. 
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Figure  46.  Simultaneous  tangential  (uT  vs.  t)  and  normal  fh  vs.  t) 

jumps  of  the  slider  during  stick-slip  motion.  Unlubricated. 
(TOLSTOI  [1967J ) . 
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As  might  be  expected,  the  influence  of  the  normal  contact  vibration 
on  the  friction  force  is  not  exclusive  of  the  stick-slip  motion. 

TOLSTOI  [l967]  and  GODFREY  [  1 9 6 7 j  report,  for  both  lubricated  and 
unlubricated  conditions,  a  resonance  fall  of  the  friction  force  when 
a  normal  forced  oscillation  of  the  slider  (or  of  the  surface  on  which 
it  slides)  is  produced  by  an  external  vibrating  source.  TOLSTOI  [  1 9 6 7J 
also  reports  a  fall  of  friction  force  due  to  the  normal  oscillations 
produced  by  the  normal  impact  of  a  small  sphere  on  the  slider. 

Large  normal  and  frictional  force  oscillations  during  smooth  unlubri¬ 
cated  sliding  of  two  steel  surfaces  are  reported  by  SOOM  and  KIM  [1983a,  b] . 
Again,  syncronism  between  those  two  oscillations  is  observed  experimentally 
and,  except  at  very  low  speeds,  those  oscillations  are  dominated  by  fre¬ 
quencies  in  the  vicinity  of  the  normal  contact  resonance. 

ARONOV,  D’ SOUZA,  KALPAKJIAN  and  SHAREEF  [1983,  1984]  observed,  with 
a  pin-on-disk  experimental  apparatus,  that  an  increase  in  normal  load 
leads  to  the  transition  from  a  steady-sliding  with  small  oscillations  to 
high  frequency  self-excited  oscillations  (intermediate  states  were  also 
observed).  Strong  coupling  among  the  different  degrees  of  freedom  (tan¬ 
gential,  normal,  rotational)  was  observed  in  the  course  of  those  self- 
excited  oscillations.  Modelling  the  sliding  system  as  a  rigid  body  with 
plane  motion  and  using  experimentally  determined  modal  frequencies  and 
damping  ratios,  those  authors  provided  an  approximate  solution  for  the 
oscillatory  motion  and  showed  its  dynamic  instability. 

The  experimental  evidence  collected  in  this  section  leads  us  to  the 
conclusion  that  an  appropriate  model  for  sliding  friction  must  incorporate 
physically  reasonable  normal  contact  conditions .  In  addition3  all  the  major 
sources  of  coupling  between  normal  and  tangential  degrees  of  freedom  should 
be  taken  in  account  in  such  a  model. 
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PART  II  -  CONTINUUM  MODELS  AND  VARIATIONAL  PRINCIPLES 

FOR  DYNAMIC  FRICTION 


In  this  part  of  our  study,  we  formulate  continuum  models  of  dynamic 
friction  .effects  in  elastodynamics  which  reflect  some  of  the 
experimentally-observed  properties  discussed  in  Part  I. 

6.  A  SIMPLE  MODEL  OF  NORMAL  INTERFACE  RESPONSE 

6.1  Preliminary  Considerations 

Before  describing  our  formulation  it  is  important  to  clarify  what  our 
major  physical  assumptions  will  be.  In  other  words,  we  need  to  explain  why, 
among  so  many  experimental  observations  and  theoretical  interpretations,  we 
will  pick  up  only  a  few  key  features  which  we  hope  will  provide  the  essential 
ingredients  for  modeling  a  large  class  of  dynamic  friction  phenomena. 

The  first  key  ingredient  of  our  model  is  a  normal  contact  constitutive  law 
which  takes  in  account  the  normal  deformability  of  the  interface  asperities. 

We  will  not  introduce  in  our  model  a  statistical  description  of  the  interface 
geometry  and  of  the  normal  contact,  which  would  be  a  possible,  defensible  and 
interesting  approach.  Instead,  we  will  introduce  a  mechanical  law  for  normal 
contact  which  does  not  contradict  the  results  of  thjse  statistical  theories. 

As  shown  in  Section  3.1,  both  theory  and  experiment  seem  to  agree  in  an  exponen¬ 
tial  type  growth  of  the  normal  contact  stresses  with  the  penetrating  approach 
for  sufficiently  high  loads.  For  sufficiently  small  loads,  which  we  will  be 
working  with,  the  vanishing  normal  stiffness  for  small  penetrations  predicted 
by  the  theory  and  the  power  laws  obtained  by  experim "mtors  are  also  not  con¬ 
tradictory.  We  will  thus  assume  that  on  the  contact  boundary  each  macroscopic 
point  "contains"  a  number  of  microscopic  interface  asperities  sufficiently  laxge 
that  a  power  law  of  the  type  (3.5,  7)  holds  locally. 
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We  incorporate  this  nonlinear  boundary  effect  in  our  model  instead  of  a 
classical  non-penetration  unilateral  contact  condition  for  the  following  reasons: 

(i)  For  small  linear  elastic  bodies  with  rough  surfaces  subjected  to  small 
loads  the  deformation  of  the  contact  interface  may  be  of  an  order  of 
magnitude  comparable  with  the  bulk  linear  elastic  deformation  of  the 
contacting  bodies. 

(ii)  The  experimental  results  of  ANDREW,  COCKBURN  and  WARING  [1967-68]  imply 
that  such  nonlinear  contact  beaavior  cannot  be  disregarded  in  dynamic 
situations  even  if  no  relative  sliding  occurs. 

(iii)  The  high  sensitivity  of  the  friction  force  to  the  normal  separation  of 
the  surface,  the  resonance  dropt--  i\  friction  force  due  to  normal  defor¬ 
mation  of  the  contacts  and  the  striking  observations  of  TOLSTOI  [1967] 
and  other  authors  described  in  Section  5.4  strongly  suggest  that  a  phys¬ 
ically  reasonable  normal  interface  model  has  to  be  used. 

Obviously,  for  highly  polished  surfaces  it  is  physically  reasonable  to 
expect  that  the  results  obtained  by  including  the  effects  of  the  normal  inter¬ 
face  deformability  should  approach  those  given  by  assuming  a  classical  unilater¬ 
al  contact  condition.  Analytical  solutions  for  frictionless  static  Hertz  con¬ 
tact  problems  suggest  that  this  is  the  case,  although  these  questions 
have  not  been  studied  from  a  mathematical  point  of  view. 

The  use  of  the  power  law  (3.5)  for  the  normal  contact  response  corresponds 
to  an  assumption  of  an  exclusively  elastic  normal  response.  Although  this  is 
the  expected  behavior  after  sufficient  rubbing  and  smoothing  of  the  interface 
it  is  questionable,  and  some  experimental  evidence  contradicts  that  this  is 
the  case  if  the  sliding  is  not  smooth  and  severe  oscillations  (accompanied  prob¬ 
ably  by  small  normal  jumps)  do  occur.  Possible  physically  consistent  ways  to 
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Incorporate  in  our  model  the  dissipative  effect  of  normal  plastic  defor¬ 
mation  on  the  interface  will  be  discussed  in  Part  III  of  this  work. 

The  other  key  ingredient  of  our  model  is  obviously  the  friction  law. 

We  will  not  provide  a  model  for  the  preliminary  (plastic)  micro¬ 
displacements  which  occur  before  gross-sliding.  Essentially  we  will  as¬ 
sume  that  these  effects  can  be  disregarded  when  compared  with  the  much 
larger  rigid-body  type  displacements  which  occur  during  sliding.  The 
classic  law  of  friction  (i)  in  Section  4.1  will  be  thus  assumed  to  hold 
locally  on  the  contact  boundary.  Only  a  small  generalization  will  be 
introduced  in  order  to  accommodate  a  possible  dependence  of  the  coeffi¬ 
cient  of  friction  on  normal  contact  forces. 

Following  the  ideas  of  TOLSTOI  [1967],  we  will  not  consider  any 
distinction  between  static  and  kinetic  coefficients  of  friction.  We 
will  try  to  show,  by  analysis  and  by  numerical  results,  that  apparent 
decreases  of  measurable  coefficients  of  friction  are  the  result  of  dy¬ 
namic  instabilities  which  are  a  consequence  of  the  inherent  non-symmetry 
of  friction  problems.  No  other  variations  of  the  coefficient  of  friction  with 

velocity  will  be  considered  either,  although  they  could  easily  be  incor¬ 
porated.  We  avoid  here  the  potentially  stabilizing  or  destabilizing 
effects  that  those  variations  would  introduce. 

Finally,  we  observe  that  creep  sliding  will  not  be  modeled  but  the  possi¬ 
bility  of  the  experimental  rate  dependence  of  the  static  coefficient  of  friction 
being  an  inertia  related  phenomenon  will  be  investigated  in  a  forthcoming  work. 
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6.2  Interface  Models 

The  Interface  between  contacting  bodies  is  a  hypothetical  medium  of 
vanishing  thickness,  the  mechanical  response  of  which  depends  upon  the 
various  geometrical  and  physical  properties  of  the  surfaces  in  contact  as 
discussed  in  Part  I.  For  the  class  of  problems  addressed  here,  we  wish 
to  characterize  the  response  of  such  an  interface  to  normal  deformations 
in  a  way  that  is  consistent  with  the  experimental,  and  statistical 
theoretic  results  summarized  in  Section  3.1. 

Consider  a  continuous  material  body  JL  ,  in  contact  with  another 
material  body  Ji-j  over  a  contact  surface  rcC3&  .  We  establish  a  fixed 
Cartesian  coordinate  system  x^,  i=l,2,3,  and  consider  motions  of  Si 
relative  to  this  fixed  spatial  frame  of  reference.  We  denote  by  u  the 

in/ 

displacement  of  particles  X  €  SL  relative  to  a  fixed  reference 

configuration,  the  configuration  occupied  by  the  "undeformed"  body, 

and  by  u.  the  Cartesian  coordinates  of  u. 

The  notion  of  a  contact  surface  F  is  worthy  of  some  elaboration.  This 

L* 

surface  should  represent  the  boundary  of  the  parent  bulk  material  of  which  the 
body  is  composed.  One  can  regard  it  as  parallel  to  a  surface  marking  average 
surface  heights.  We  suppose  that  T  has  a  well-defined  exterior  normal  vector 
n  and  that  the  actual  interface  (asperities,  oxide  film,  gas,  work-hardened 
material,  etc.)  is  initially  of  thickness  tg,  as  shown.  The  initial  gap  be¬ 
tween  F  and  the  other  body  that  may  come  in  contact  with  Jl  is  defined  as  the 

w 

distance,  measured  along  a  line  normal  to  T  ,  between  the  highest  as- 
perities  of  the  bodies  in  the  reference  configuration.  The  interface 
thickness  after  deformation  is  denoted  t  in  Fig.  47  and  the  actual  displacement 
of  F„  in  the  direction  of  n  is  u  =  u»n.  Thus,  the  approach  of  the  material 

V  ~  T1  ~ 

contact  surfaces  is 


a 


t0  -  t  -  (un  -  g>+ 


(6.1) 


This  definition  is  valid  whenever  the  body  is  rigid  and  bounded  by 

a  surface  opposing  T  which  is  ideally  flat,  as  indicated  in  the  figure. 

I* 

This  model  is  easily  generalized  to  cases  in  which  is  deformable  and 
rough  (see  KIKUCHI  and  ODEN  [1985]). 

Let  a  denote  the  Cartesian  components  of  the  Cauchy  stress  tensor  at 

a  particle  X  in  li.  The  normal  stress  developed  on  Tr  is  then 

a  =  O, .n  n  (6.2) 

n  ij  i  ] 

where  the  standard  summation  is  used. 

The  constitutive  properties  of  the  interface  are  characterized  by  a  re¬ 
lationship  between  and  a.  In  Fig.  48  we  have  plotted  the  variation  in  the 
normal  force  N  on  a  unit  block  of  the  material  with  the  penetrating  approach. 
Except  for  very  large  penetrations,  where  the  response  is  essentially  expo¬ 
nential,  a  large  body  of  evidence  points  to  a  power-law  relation  between 
(=N/unit  area)  and  a: 

m 

-a  =  c  (u  -  g)  n  (6.4) 

n  n  n  + 

Here  c  and  m  are  material  parameters  independent  of  a. 
n  n 

Much  of  the  experimental  evidence  described  earlier  indicated  that  the 
material  parameters  cn  and  m^  can  be  easily  determined  from  interface  compression 
experiments.  Other  effects  could  be  added  to  this  interface  law,  but  it  is 
sufficiently  general  to  model  the  principal  mechanisms  of  interest  here. 

Remarks.  1.  The  form  of  the  interface  constitutive  equation  (6.4)  is  consis¬ 
tent  with  the  experimental  observations  outlined  earlier  (recall  Section  3.1): 
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N 


N  oc  exp(a) 
(exponential 
law) 


a 


Figure  48.  Normal  force  versus  approach  relationship  suggested  by 
experimental  and  statistical  theoretic  results . 

Portion  A,  "light"  normal  loads  characteristic  of  sliding 
interfaces. 

Portion  B,  "heavy"  normal  loads  characteristic  of  fixed 


interfaces. 
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da 

a)  —a  |  .  0  («-(t.n-  g)+) 

a=0 


m 

b)  -a  aa  n  with  2  <  m  <3.33 
n  —  n  — 


2.  Finite  element  computations  of  static  contact  problems  using  such  a 
normal  interface  constitutive  equation  were  first  done  by  BACK,  BURDEKIN  and 
COWLEY  [1974], 


7.  A  PROBLEM  IN  ELASTODYNAMICS 

We  shall  now  formulate  a  class  of  initial-value  problems  in 

elas to dynamics  which  include  sliding  friction  effects.  We  begin  by 

considering  a  linearly  elastic  body,  the  interior  of  which  is  an  open 

bounded  domain  ft  in  ]R  (N  =  2  or  3)  with  boundary  T  consisting  of 

portions  1^,  1^,  and  1^,  ,  as  indicated  in  Fig.  49  .  Displacements 

u  are  prescribed  on  ,  tractions  t  are  prescribed  on  Tp  ,  the  body 

may  be  elastically  supported  by  springs  of  modulus  k..  on  ,  and  T- 

ij  t  C 

is  the  candidate  contact  surface.  We  suppose  that  the  body  may  come  in 
contact  with  a  rough  sliding  belt  (a  moving  foundation)  which  slides  by 
the  material  contact  surface  with  a  velocity  U  .  This  foundation 

is  initially  (at  t  =  0)  at  a  gap  distance  g  from  While  strains  in  the 

body  are  small,  certain  large  rigid  motions  relative  to  the  reference  config¬ 
urations  are  permitted  in  our  theory.  For  simplicity,  the  body  is  assumed 
to  be  composed  of  a  Hookean  material  characterized  by  the  classical  consti¬ 


tutive  equation, 
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Figure  49.  Geometry  and,  notation  used  in  elastody  1 amics  'problem 
with  sliding  friction. 
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aij(~)  "  Eijk£uk,£ 


(7.1) 


where  u^  £  =  ,  1  £  i,j,k,£  £  N  and  the  elasticities  satisfy 


Eijld  Ejik£  "  Eij£k  =  Ek£ij 


EijW,£  L  ^ 


Eijk£  Ak£  Aij  ’  -  aAijAij 


(a.e.  in  Q,  for  every  symmetric  tensor  , 

a  being  a  positive  constant) 


(7.2) 


Wll  ±M 


For  a  time  interval  [ 0 ,T] ,  the  equations  governing  this 

elastodynamics  problem  are  grouped  as  follows: 

Linear  Momentum  (Constitutive  Equations,  Angular  Momentum,  and  Mass 
Conservation) 


aij(u),j  +  bj^  =  Pi^  in  ft  x  [0,  T] 


(7.3) 


where 
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bj  -  the  components  of  body  force  per  unit  volume,  assumed  to 
be  sufficiently  smooth  functions  of  x  =  (x, ,  x_,...,x  ) 
p  =  the  mass  density,  ^  =  0  , 

ot 

OO 

P  6  L  (ft)  ,  p  >  p_  >  0 

°  9?u- 

Uj  =  particle  acceleration  = 

3t2 


(a«  ■  V 


Boundary  Conditions  on  Tp,  rp,  rp 
.D 


"!  =  "i 


on  rD  x  10, T] 


aij^)nj  =  ti  on  Ff  x 


■  'Vuj  -  D 


on  r£  x  [  0 , T] 


(7.4) 


D  E 

where  U  ^ ,  t^,  and  U  ^  are  given  sufficiently  smooth  functions. 


Initial  Conditions 


u  =  u. 


u  =  u. 


>  in  at  t  =  0 


(7.5) 


with  Uq,  Uj  given  smooth  functions  of  x  . 

Normal  Interface  Response.  Recalling  (6.4),  we  have 


m 


a  (u)  =  -c  (u  -  g) 
n  ~  n  n  + 


n 


on  rc  x  [0,T] 
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Friction  Conditions.  Let  a  be  the  tangential  (frictional)  stresses 
developed  on  the  contact  surface  (a,^  =  a^jnj  ~  °nni^  *  ^en  we  employ 

the  following  friction  law: 


u  <  g 
n  s 


a  (u)  -  0 

<vl 


\>s 


|aT(u) 1 

< 

ct(V8)+T 

and 

lgT(u)| 

< 

> 

(7.7) 

and 

l2T(u>l 

= 

m 

CT<'Un"s)+T  ^a^-0’ 

4  -  D„  =-Aa 
~T  ~T  ~T 


on  rc  x  [0,T] 


Here  c^  and  mT  are  material  constants  depending  on  interface 
properties  (indeed,  characterizing  further  constitutive  properties  of  the 

interface),  is  the  tangential  velocity  of  material  particles  on  Tc  , 

•  C 

and  UT  is  the  prescribed  tangential  velocity  of  the  foundation 

(adjacent  body)  with  which  comes  in  contact. 

The  friction  law  (7.7)  is  a  generalization  of  the  Coulomb's  friction 

law,  which  is  recovered  if  m  =  m„  .  In  such  a  case,  u  =  cm/c  is  the 

n  T  T  n 

usual  coefficient  of  friction.  The  law  (7.7)  allows  for  a  dependence  of 

the  friction  coefficient  on  the  normal  contact  pressure  (recall  Sections  4.1,2). 
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We  remark  that  other  possible  choices  would  also  have  been  easy  to  Intro¬ 
duce  at  this  point.  The  one  presented  above  can  be  interpreted  as  allowing  for 
deviations  of  Amonton's  first  law  due  to  a  variation  of  the  true  area  of  contact 
with  a  power  of  the  normal  stress  level  different  from  unity 

y  =  C|anja 

“t 

with  a  - - 1  and  C  = 

m 

u 

n 


8 .  VARIATIONAL  FORMULATIONS 


8.1  A  Nonlinear  Hyperbolic  Variational  Inequality.  We  now  consider 
the  formulation  of  a  variational  principle  corresponding  to  the 
elasto-dynaraics  problem  outlined  above.  Following  steps  similar  to  those  of 
DUVAUT  and  LIONS  [1976],  the  nonlinear  elastodynamics  problem  can  be  shown 
to  be  formally  equivalent  to  the  following  variational  problem: 

Problem  1. 

Find  the  function  t  ->  u(t)  of  [0,T]  V  such  that 

N 


<u(t)  ,  v  -  u(t£>  +  a(u(  t) ,  v  -  u(t)) 

^  /v  »v  IV  IV 

+  <J>(u(t)),  V  -  u(t^> 

«>N/  <V  (V 


(8.1) 


+  j(u(t),  v)  -  j(u(t)  ,  u(t)) 
><^(t)  ,  V  -  u(t)>  V  v£  v 


with  the  initial  conditions: 


u(0)  =  u 
'i(O)  =  Uj 


\ 


(8.2) 
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We  have  assumed  here  for  simplicity,  that  p  =  1  and  that  U*T 
The  following  notations  and  definitions  were  also  used: 

V  =  the  space  of  admissible  displacements  (velocities) 

=  {v  =  (Vj,  v2,...,vN)|viC  H1^),  y(vj[)  =  0  a.e.  on  1^} 


=  0 


a  :  V  x  V  -*■  U 

a(u,  v)  =  a_(u,  v)  +  a.(yu,  yv) 

V  (M  M  -*•  (V  fV 


(8.3) 


where: 


a_(u,  v)  =  virtual  work  (power)  produced  by  the  action  of  the 

U  r*  ~ 

stresses  cr,(u)  on  the  strains  (strain  rates)  e . .  (v)  = 
ij  ~  iJ  ~ 

(vi,3  +v3,l)/2> 


=  f  ai,(u)  e^v)  dx 
J  ~ 

-Ja  Ei«ukAvi,j dx  “’“G  v 


(8.4) 


where  u  is  a  displacement  and  v  is  a  virtual  displacement  (velocity). 

a^(£.  ri)  =  virtual  work  (power)  produced  by  the  deformation  of  the 
linear  springs  on  r_. 


=  /r  Kij?j  ds  5*ne(L2(rE))N 


(8.5) 


P  :  V  -*  V’ 

<P(u),  =  virtual  work  (power)  produced  by  the  normal  contact 

pressure  on  the  displacement  (velocity)  v, 

r  m 

=  Jr,  c  (u  -  g)  n  v  ds 
n  n  +  n 


(8.6) 
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j  :  V  x  V  ->  m 

j(u,  v)  *  virtual  power  produced  by  the  frictional  force  on  the  velocity  v  , 


f(t)€ 


* L  ct( vg)  tLvH?i  ds 

X 

V' 


(8.7) 


i 

i 

i 


<f(t),  v>  *  virtual  work  (power)  produced  by  the  external  forces  (body 
forces,  prescribed  tractions,  initial  deflection  of  the 
linear  springs)  on  the  displacement  (velocity)  v. 


fB  biVx  +  /  v<vi>ds  + 
+  L  vEjY  (vds 


(8.8) 


Here  <^»  denotes  duality  pairing  on  V'  x  V  where  V'  is  the 

1  N 

topological  dual  of  V  ;  y  is  the  trace  operator  mapping  (H  (ft))  onto 

L  fj 

(H  (ft))  which  may  be  decomposed  into  normal  components  y  n(v)  and 
tangential  components  y^,(v)  .  For  simplicity  of  notation,  the  latter  are 
denoted  as  v^  and  v^  ,  respectively.  We  also  observe  that  the  boundary 
integrals  on  Tc  are  well  defined  for  1  £  m^  ,  m^,  <_  3  if  N  =  3  and 
for  1  <_  m^  ,  ntp  if  N  =  2  ,  because,  for  e  H*(ft),  y(v^)  £  Lq(T)  ,  with 
1  S  q  —  4  ,  for  N  =  3  ,  and  with  1  <  q  ,  for  N  =  2  .  We  can  thus  con¬ 
clude  that  for  the  case  N=2  all  the  experimentally  observed  values  of  m^ 

(m^  6  [2,  3.33],  recall  Section  3.1)  can  be  used  in  a  mathematically  consis¬ 
tent  way  in  our  formulation.  For  the  case  N=3,  the  mathematical  restriction 

m  <  3  will  be  harmless  in  most  of  the  cases:  in  BACK,  BURDEKIN  and  COWLEY 
n  — 

[1973]  only  one  experimental  value  of  m^  is  shown  (3.125)  which  slightly  ex¬ 
ceeds  3,  the  most  common  value  for  m  being  2. 

*  n 


i 


124 


Of  course,  for  a  complete  definition  of  Problem  1  we  would  need  to 
make,  precise  the.  assumed  regularity  of  the  displacements,  velocity  and 
acceleration  relatively  to  the  time  variable  t  ,  i.e.,  we  would  need  to 
specify  the  spaces  of  abstract  functions  to  which  the  mappings  :  u(t)  , 
u(t)  :  [ 0 , T]  V  and  ii(t)  :  [0,T]  V*  belong.  Intimately  related  to 
this  question  is  the  need  for  the  initial  conditions  (8.2)  to  make  sense 
and  the  need  to  be  precise  in  what  sense  the  variational  inequality  (8.1) 
is  satisfied  along  the  time  interval  (0,T)  . 

These  questions  will  not  be  addressed  here.  Instead,  we  will  proceed 
formally,  and  seek  a  computational  algorithm  to  obtain  approximate 

solutions  for  the  problem  above. 

8.2  A  Regularization  of  the  Friction  Functional.  Our  first  goal  is 
to  approximate  Problem  1  by  a  family  of  regularized  problems  which  lead 
to  the  solution  of  a  variational  equation  instead  of  a  variational 
inequality  (as  (8.1)). 

Toward  this  end,  we  approximate  the  friction  functional  j  :  V  x  V 
H  which  is  nondif f erentiable  in  the  second  argument  (velocity)  by  a 
family  of  functionals  j  convex  and  differentiable  on  the  second 
argument : 


j  :  V  x  V  ->  m 
Je 


C  *r 

j£(u,  v)  =J  cT(un  -  g)+  \|>e(vT  -  UT)  ds 


(8.9) 


N  J 

where  the  function  \};  :  (L^T.,))  ->Lq(r_)  is  an  approximation  of  the 

C  C  L* 

function  |«|  :  (Lq(rc))N  -*■  1^(1^)  and  is  defined  for  e>0  ,  £e(Lq(rc)) 
and  a.e.  x  6  ^  (  according  to  : 


N 


e|^|  (1  ’if)  if  I? (x) |  £  e 
e(ljl  "  j)  if  I  ^*>1  >  e 


(8.10) 


6 


Remark:  Since  for  a.e.  x  e  rc  ,  Ke  (£  (x))  |  <  |£(x)|  and  since  |£  | 
Lq(T  )  then  ip  (E)  e  Lq(rr)  as  stated  above  and  no  additional  restrictions 
on  are  required  for  the  integral  on  (8.9)  to  be  well  defined.  Q 


The  partial  derivative  of  j  relative  to  the  second  argument,  at 
(u,  w)  in  the  direction  of  v,  is  then  given  by: 


<5_(u>w),  \r>  =<$?j  (u,w) ,  v> 


=  /rcCT(Vg)+  IV”T"HT)(^T)]dS 
Vu,  v,  w  €  V 


(8.11) 


Q  N 

where,  for  e>0  ,  £  ,  ^  6  (LM(yc)  )  and  a.e.  x  d  Tq  , 


MCXn)  =  ip;<0(n) 

£  #s»  #s/  £  ^  ~ 


7“ 

\  (2  -  |f|)(£-n)  if  |C(x)l<  e 
\t\  (C*n)  if  U(x)  |  >  e 

I  c,  I  ~  ~  ~  ~ 


(8.12) 


is  the  directional  derivative  of  at  £  in  the  direction  of  n  . 


Remark:  From  the  definition  it  is  easy  to  see  that  for  a.e.  x  £  r.,|  [<J>  ( ^ ) 

.  ~  t  M 

(n)](x)j  £  2|n(x)  |  and  consequently  [4>e(0(n)]  ^  Lq(rc)  .  The  integral 
on  (8.11)  is  then  well  defined  without  any  further  restrictions  on  mT  .  Q 
We  now  define  the  regularized  form  of  Problem  1: 

Problem  1_: 

1  C. 


Find  the  function  u  (t)  of  [0,T]  ->•  V  s.t. 

~£ 


<Cu£(t),  V>  +  a(u£(t),  v) 


+  <!P(ue(  t) ) ,  £>  +  <?e^£(t),  u£(t)),  v> 


=  <f  (  t)  ,  v>  \/  v  C  V 


(8.13) 


with  the  initial  conditions 


E 


E 


We  observe  that  now  we  have  a  variational  equation  instead  of  a 
variational  inequality.  We  also  observe  that  this  problem  is  formally 
equivalent  to  the  problem  defined  by  a  system  of  equations  analogous  to 
(7.3)  -  (7.6)  with  u£  replacing  u  and  again  p  =  1  ,  E  0  . 
However,  the  regularized  friction  conditions  on  T  are  now  of  the  form: 


Regularized  Friction  Conditions 


.  *r  •  *C 

u  -tC  u  -If 

(2  - 

if  <  c 


,(u,  fl)  =  -^(u^g)^ 


(8.15; 


*  ,'rc 

u  -  UT 

A* 

ii^r 


lif  '  '  >  e 

In  order  to  visualize  the  regularization  procedure  used  above,  we 
consider  now  its  particularization  for  the  case  of  a  two  dimensional  (N=2) 


domain  ft  with  a  boundary  T  sufficiently  smooth  that  we  can  define  a 
unit  vector  i  tangent  to  T  .  In  this  case  each  vector  £  tangent  to 

~1.  V>  <"w 

rc  is  determined  by  the  real  number  £  such  that  £  =  £  i^.  .  The 

functions  and  <)>e  =  are  then,  essentially,  real  valued  functions 

of  a  real  variable,  defined  by: 


- 


e|f|2(l-j  |f|)  if  \K\±e 


[  e( if  I  "  T>  if  1^1  >  e 

f  (2-|||)  |  if  |C|  <  e 


(8.16) 


f(£) 

sgn(C) 


€ 


c 


The  graphs  of  these  functions  are  depicted  in  Fig.  50  together  with  the 
graphs  of  the  functions  which  they  approximate:  the  function  |»|  and  its 
derivative  in  1R  -  {0},  sgn(»)  . 

Finally  we  observe  that  a  physical  interpretation  might  be  attributed  to  the 
above  regularization  procedure.  The  increasing  branch  of  the  <j)£(C)  curve  is 
essentially  of  the  same  nature  as  those  experimentally  observed  during  creep 
sliding  (recall  Section  5.1.2).  However,  we  prefer  to  regard  it  exclu¬ 

sively  as  a  mathematical  procedure  to  obtain  an  approximate  problem  which  cor¬ 
responds  to  a  variation  equation.  We  do  so  because  the  values  of  £  which  will 
be  used  in  actual  computations  will  not  be  so  small  as  those  required  for  the 
above  physical  interpretation  to  be  valid  in  the  case  of  hard  metals  with 
which  we  will  be  considering.  The  choice  of  £  will  be  dictated  only  by  the 
desired  proximity  of  the  solutions  of  problems  1  and  1^  and  the  corresponding 
computational  costs  associated. 

Many  examples  of  similar  regularization  procedures  in  dynamic  friction 
problems  can  be  found  in  the  literature.  Among  them  we  refer  to  DUVAUT  and 
LIONS  [1972.1  for  mathematical  aspects;  to  THRELFALL  [1978]  and  ROONEY  and 
DERAVI  [1982]  for  computational  applications  to  Theory  of  Mechanisms  and  to 
MARTINS  and  ODEN  [1983]  for  a  finite  element  analysis  of  a  simplified  friction 
problem. 

We  are  aware  of  some  of  the  difficulties  associated  with  these  regulari¬ 
zation  procedures:  a  perfect  zero  relative  velocity  between  the  contacting 
surfaces  is  not  possible  —  when  "stuck"  the  contact  surfaces  will  creep  with 
a  relative  velocity  smaller  than  £;  trying  to  impose  a  very  small  £  will  im¬ 
ply  the  need  to  reduce  substantially  the  steps  of  the  time  integration  scheme, 
particularly  in  unloading  situations  —  transitions  from  "slip"  to  "stick". 
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PART  III  -  FINITE  ELEMENT  MODELS  AND  NUMERICAL  ANALYSIS 


9.  FINITE  ELEMENT  APPROXIMATIONS  OF  THE  REGULARIZED  PROBLEM 


Using  standard  finite  element  procedures,  approximate  versions  of  Prob¬ 


lem  1  can  be  constructed  in  finite  dimensional  subspaces  V.  (CVCV1)  . 
e  n 


For  a  certain  mesh  (h)  the  approximate  displacements,  velocities  and 


accelerations  at  each  time  t  are  elements  of  V,  : 

n 


vh(t),  vh(t),  vh(t)  £  V 


Within  each  element  (esl ,2, . . . ,E^)  the  components  of  the 


displacements,  velocities  and  accelerations  are  expressed  in  the  form: 


e  I 


v.(x,t)  =  £  v.(t)N  (x) 

J  ~  J=  ]  J  ^  ~ 

h  Ne  T 

v.(x,t)  =  £  v.(t)  N  (x) 

J  ~  1=1  J 


v^(x,  t)  =  £  v?(t)  N  (x) 

J  ~  1=1  J 


where  j=l,2,...,N;  N  =  number  of  nodes  of  the  element  ;  v*(t),  v*(t), 

J  J 

••1/  .  \  _  ■  1  _  3  _  -t  i  r  •*  .  •  .  .  .  .  •  >  i 


Vj(t)  are  the  nodal  values  of  the  displacements,  etc.,  at  time  t  and 


is  the  element  shape  function  associated  with  the  node  I  . 


The  finite  element  version  of  Problem  ]g.  is  then: 


Problem  1„  : 


Find  ’■he  function  t -►  u  (t)  of  lO,T]  s.t. 

~e 


(i£(t),  vh)  +  a(u^(t),  vh)  +  <CP(uh(t)),  vh^> 


rAflstsntnsm-HJ’  r  n  >lh.m  i  k-n  «  «  t  mi  x  i  \  i>  m  w  •»  t*  . tv  ^  ^ 
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T 

«' 

\ 


+  <\^(u£(t),  u^(t)),  yh>=  </f(t),  vh^> 


Vvh€  vh 


(9.1) 


with  the  initial  conditions 


(u^(0),  vh)  =  (uQ,  vh) 


(u^(0),  vh)  =  (u.^  vh) 


V  Vhevh 


(9.2) 


If  is  the  number  of  nodes  of  the  finite  element  mesh,  then  this 


problem  is  equivalent  to  the  following: 


Find  the  function  t  -*•  r(t)  of  [ 0 ,T]  -*■  1R  ,  s.t. 


M  r  (t)  +  K  r(t)  +  P(r(t))  + 


+  J(r(t) ,  r(t))  =  F(t) 


(9.3) 


with  the  initial  conditions 


r(0)  =  p, 


5<0>  =  Ei 


(9.4) 


Here  we  have  introduced  the  following  matrix  notations: 


r(t),  r(t),  r(t)  =  the  column  vectors  of  nodal  displacements,  velocities 


and  accelerations,  respectively 


=  standard  mass  matrix 


=  standard  stiffness  matrix 


=  consistent  nodal  force  vector 


P  (r(t))  =  vector  of  consistent  nodal  normal  forces  on  F 


M-.xnm  i*  mm  w  '~j>  r  rfJV*/'  V)  VI  W*\*»  UV  l/VVV  Lf¥  UKU  XV*  W> 
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J(r(t),  r(t)) 
*0^1 } 


vector  of  consistent  nodal  friction  forces  on 
initial  nodal  displacement  (velocity) 


C  6  ) 

The  components  of  the  element  vector  P  are  of  the  form: 

(e). 


’p  -  -/(e),,  a  n  ds 
j  c 


(9.5) 


<e) 


and  the  components  of  the  element  vector  'J  are  of  the  form: 

*>w 

(e) 


'  '/(e> r.  Ni 


ds 


(9.6) 


10.  ALGORITHMS  FOR  NONLINEAR  DYNAMICAL  SYSTEMS 


The  algorithms  that  we  shall  use  for  solving  the  discrete  dynamical  system 
involve  variants  of  standard  schemes  in  use  in  nonlinear  structural  dynamics 
calculations.  An  outline  of  these  computational  procedures  follows. 


10.1  A  Newmark-Type  Algorithm. 

We  begin  with  the  presentation  of  an  implicit  Newmark-type  algorithm  that 

has  proved  to  be  effective  for  many  problems  in  the  general  class  under  study. 

Let  P  denote  a  partition  of  the  time  domain  [0,  T]  into  M  intervals  of 

length  At  such  that  0  =  t^  ,  t^,  ...,  t^  ,  ...»  t^  =  T  with  -  tR  =  At. 

We  shall  approximate  the  velocities  and  accelerations  at  time  t^  by  expressing 

them  as  functions  of  the  displacements  velocities  and  accelerations  at  time 

t  1  and  of  the  displacements  at  time  t  ,  by  the  following  relations: 

K.— _L  K 


’  =  }  +  (1  "  fe-l  +  At(1  ‘ 


U~K  3At 


0At‘ 


^"ilK-1^  "  RAt  i-1  (2B  "  1>1?K-1 


(10.1) 


-  h. 


where  y  and  B  are  the  so-called  Newmark  parameters  and  u^,  -  u£(tj,),  etc*  are 
used  in  the  foregoing,  if  no  confusion  is  likely  to  arise. 
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Introducing  the  above  relations  into  the  variational  equation  (9.1  ),  we 
obtain  the  following  variation  equation  at  time  t  : 


BAt 


[uR,  v)  +  a(uR»  v)  +  <P(uR),  £> 


+  £>  =  <&>  £>  VY£vh 


(10.2) 


Here 


<?<Sk)’  V’ 


<?&>  sfe  +  «  -P“K-1  +  At(l  -X)  UK_X)  ,  v> 


<fK-  !>'<£k’  v>  + 


+  -Vi +  8sr  Vi  +  -  ‘Ivr  ;> 


The  above  variational  equation  can  be  put  in  the  operator  form: 


V-V  -  “  1:1  v'h 


(10.3) 


with  being  the  nonlinear  operator  defined  by 


<V!!k)>  l>  V.  XV’  -  -S’  (V  v>  +  a(V  v) 

h  h  BAt  ~  ~ 

+  <P(uR),  V^+C^J^),  v>> 

■<-£k>  j>  Vuk>  j  e  vh 


(10.4) 
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Let  K^u^)  =  DRj^u^)  £  j£(Vh,  Vh^  136  the  derivative  of  the  map 


at  uK  £  .  Then  the  Newton-Raphson  iteration  technique  for  solving 

/in  n\  _ rr-11 _  _  i-i ..  —  J ~ 


(10.  3)  is  as  follows:  given  a  starting  value  u  6  V  ,  obtain  sue- 

•vlX 


cessive  approximations  of  the  solution  u^  by  using  the  recurrence  formula, 


Hk<1+1)=  ~R(i)  '  IVHK<1)>rl  rk<“k1>) 


(10.  5) 


where  i=0,l,...  is  the  iteration  counter.  Thus,  at  each  iteration  i  , 


the  following  variational  equation  has  to  be  solved: 


)(A~K  ^’^V  xV'  -  <[Rk(uk^  h,  £>v  xV  » 
h  n  n  n 


V:"i  uo.6) 


a  (1)  (i+1)  (i) 

with  Au  =  u  -  u 

'vKk  *vl\  lx 


In  the  present  problem  we  have,  for  u  ,  w  ,  v£  v. 


xV'  " 
n  n 


=  a(w,  v)  +  - j  (w,  v)  +  <D  P(u  >(w),  v>  +  <J>  J(u„)(w),  v> 

~  ~  8At^  ~  ~  ~  ~  *  ~K  ~ 


where 
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=  J  CTmT(unK_g)+T  l<{>e(~TK_HT)(^T)lwndS 
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where,  for  e  >  0  ,  | ,  n  ,  £  e  (Lq(rc))N  and  a.e.  x£  p 


*  (5)(n,0  =  <i>"U)(n,e) 

£  ~  ~  ~  t  —  ~  ~ 
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f  -7^j-(«  •  D)(S-p  +  7<2  -ifptn  •  c> 
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(£  *  n)(£  •  O  *  i*^j-  (n  •  O 


if  |£(x)|  >  e 


(10.7) 


is  the  second  directional  derivative  of  at  E  on  the  directions  of 


0  and  C . 

Remark:  If  N  =  2  the  ahove  definition  reduces  to 


f  O  -  £|Xn-?>  if  l§C5)l  <  e 


tt  (0(n,c)  -  ^ 

C,  *v  <v 


0  if  |£(x)|  >  e 


for  e  >  0  and  E  ,ri  ,£  £  (kq(rr))  .  □ 


We  next  recast  equation  (10.6)  in  matrix  form: 


L(1)  AfK(i>  =-~RK(i> 


where 


hi 


(i)  . 


5  +  ^  +  ^1>+^1)+*iT<1) 
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pAt 

The  element  components  of  the  matrices  Kn  ,  K^n  ,  and  for  the  case 

N  =  2  are  of  the  form: 


<e)K.n 


m  -1 


M.N,  =  f  c  m  (u  -g),  n  n.n.NMNv,  ds 
i  j  J(e)p  n  n  n  6  +  i  j  M  N 


(e)r 


where  q  is  the  part  of  rc  belonging  to  element  e: 
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M  N. 

i  1 
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where  ^e^^(<)  is  the  part  of  where 


uT-UTj  £  e 


Physically,  the  matrix  K  produces  a  variation  of  the  normal  nodal 


forces  due  to  a  variation  of  the  nodal  normal  displacements;  the  matrix 
Tn 

K  produces  a  variation  of  the  nodal  friction  forces  due  to  a  variation 

T 

of  the  nodal  normal  displacements;  finally,  the  matrix  C  produces  a 

variation  of  the  nodal  friction  forces  due  to  a  variation  of  the  nodal 

sliding  velocities. 

n  T 

Both  K  and  C  are  symmetric  matrices  because  both  are  associated 
with  terms  of  <^Kj,(u^)(w),  which  can  be  considered  as  second  derivatives 
of  functionals.  In  fact, 


<DP(u)  (w),  £>  =  <J>2p(u)(w),  £> 

for  u,v,w  6  V  ,  with  p  :  V  -*■  R  being  the  deformation  energy  associated 
with  the  nonlinear  springs  on  T  : 


/ 


We  also  have. 


<^32Je^u*  ^ 


<C32  -Mu>  u> 


£  ~  ~ 


is  the  second  partial  derivative  of  the  friction  functional  j  relative 


to  the  second  variable  (velocity).  We  also  observe  that,  in  the  two-dimen¬ 


,T  , 


sional  case  (N=2)  C  is  the  incremental  damping  matrix  associated  with 


the  region  of  T  which  is  "stuck"  at  some  instant  (the  modulus  of  the 

v 


sliding  velocity  is  smaller  than  e  )  .  Finally,  we  observe  that  the  matrix 


Tn 

K  is  not  symmetric.  This  results  from  the  fact  that,  according  to  the  nor 


mal  contact  behavior,  the  normal  stresses  on  depend  only  on  the  pene¬ 


tration  (un~g)+  .  A  change  on  the  sliding  velocity  does  not  produce  any 


change  on  the  normal  contact  forces 
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10.2  Semi-Explicit,  Central  Difference  Scheme 


It  is  well  known  that  the  central-difference  technique  is  equivalent  to 
the  Newmark  method  for  3  =  0,  y  =  1/2.  This  special  case  of  the  algorithm  described 


above  is  now  outlined.  The  accelerations  and  the  displacements  at  time  t  are 

K 


related  with  the  displacements,  velocities  and  accelerations  at  time  t  and 

K— 1 


with  the  velocities  at  time  t  ,  according  to: 

K 


~K-1  +  At  ^~K  ~K-l'; 


At 


(10.8) 


“K  =  Hk-1  +  At  "k-1  +  —  ~K-1 


Introducing  these  relations  into  (9.1),  the  following  variational  equa¬ 


tion  is  obtained  at  time  t  : 

K 


At  (^K*  +  <5  (“K}  •  Z>  '  <£k* 


V  v  6  Vh  (10.9) 


where,  now: 


(~K>  *  Z>  E  *  ~‘> 


At 


"  <¥"k-i  +  At!iK-i  +  T  Hk-i*  “*>•  Z> 


(10.10) 


<£k,  v>  =  <C£k,  y> 
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-  a(“K-l  +  At  "K-l  +  2  Hk-1’  Z> 


At 


-  <P(“K-1  +  At  UK-1  +  2  ~K-1 


V) ) .  v> 


(10.11) 


It  is  important  to  observe  that  the  second  equation  in  (10.8)  determines 


uniquely  the  displacement  at  time  t  from  the  values  of  the  velocity  and  ac- 

K 


celeration  at  the  previous  time  t  . .  The  nonlinearity  involved  in  equation 

K— ± 

(10.9)  results  only  from  the  nonlinear  dependence  of  the  frictional  forces 


on  the  tangential  velocities.  In  the  iterative  process  required  to  solve 
(10.9),  the  displacements  (and  consequently  the  normal  stresses)  will  thus 
remain  constant. 

Using  again  the  Newton-Raphson  method  to  solve  equation  (10.9)  for  the 
unknown  velocities  u  ,  the  following  linear  equation  is  obtained  at  each 
iteration: 

<KK(u£l))(Au£i)),  y>  -  -  <Rk(u£1)),  v>  V  y  6  vh  (10.12) 


In  this  case, 


~K 


y> 


<RK(yK) »  = 


(w,  v)  +  <^DJ  (uR)  (w)  ,  v> 

Vv  r e  vh 

,  y)  +  <j(uK),  v>-<?K,  v > 


V  6  V 
~  n 


<DJ(u  )(w),  v>  =  <3  J.(u  ,  u  )(w),  v> 

~1\  ~  Z  «vl\  <vl\  «v 


(10.13) 


The  discrete  variational  equation  leads  to  the  following  system  of  al¬ 
gebraic  linear  equations  to  be  solved  at  each  iteration: 


~(i)  A.(i)  _  _D(i) 


K 


(10.14) 


cT(i) 


where 
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~K 
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+  At  r  ,  + 
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+  At  r„  ,  + 

At 
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Some  computational  advantages  of  this  algorithm  over  the  Newmark  scheme 

described  earlier  are  the  following: 

T(i)  ~(i) 

(i)  Both  M  and  (hence  )  are  symmetric  matrices. 

(ii)  If  a  diagonalized  mass  matrix  M  is  used,  the  only  nondiagonal  entries 

-  —  - 

of  result  from  the  points  on  which  are  "stuck".  The  structure 

~(i) 

of  for  an  appropriate  numbering  of  the  nodes  of  the  mesh  on  Tc 

suggests  that  the  use  of  a  symmetric  skyline  solver  for  the  solution  of 
(10.14)  may  prove  to  be  very  efficient. 

(iii)  If,  in  addition  to  the  diagonalization  of  the  mass  matrix: 

(iii-1)  a  rotation  of  the  degrees  of  freedom  on  is  performed 
in  such  a  way  that  the  final  degrees  of  freedom  have  the 
direction  of  the  normal  and  the  tangent  to  F,-.  at  each  con¬ 
tact  node; 

T 

Ciii-2)  the  element  contributions  to  C  are  computed  with  a  quad¬ 
rature  rule  which  uses  the  nodes  as  integration  points 
(e.g.,  trapezoidal  rule  for  linear  elements  and  Simpson's 

rule  for  quadratic  elements) ; 

T(i) 

then  the  matrix  C  is  also  diagonal  and  the  solution  of  (10.14) 


becomes  very  simple 


140 


Remark:  The  discontinuity  of  the  Coulomb's  friction  law  at  zero  sliding 
velocity  is  a  major  source  of  computational  difficulties  in  friction  problems. 
Even  though,  in  the  algorithms  described  in  this  and  the  previous  section,  a 
regularized  form  of  that  law  is  used,  those  difficulties  cannot  be  completely 
avoided.  The  situation  which  may  arise  when  using  the  methods  described  here 
with  a  constant  time  step  is  the  following:  in  unloading  situations  (passage 
from  sliding  to  adhesion)  the  Newton-Raphson  iterative  techniques  may  fail  to 
converge  if  e  is  very  small  and  the  time  step  too  large.  This  situation  does 
not  seem  to  be  dependent  on  the  specific  iterative  technique  used  or  on  the 
type  of  regularization  of  Coulomb's  law  adopted.  In  a  previous  work  (MARTINS 
and  ODEN  [1983]),  we  found  similar  problems  though  we  used  a  different  iter¬ 
ative  scheme  (successive  approximations)  and  a  different  regularization  (piece- 
wise  linear  <j>^).  Our  numerical  results  suggest  that  for  small  values  of  £ 
the  radius  of  convergence  of  the  iterative  scheme  used  is  very  small  due  to 
the  steep  changes  in  (j)£  on  the  interval  [-£,  e) . 

The  critical  situations  arise  in  transitions  from  sliding  to  adhesion 
because  it  is  then  that  the  most  important  changes  in  the  solution  occur.  We 
observe  that,  for  the  nonregular ized  problem,  a  discontinuity  of  the  tangen¬ 
tial  acceleration  arises  when  a  point  on  the  contact  surface  ceases  to  be 
sliding  and  a  period  of  adhesion  is  initiated. 

One  simple  remedy  for  these  difficulties  is  to  decrease  the  time  step 
until  two  successive  solutions  are  not  too  far  apart.  Obviously  this  increases 
the  computational  costs  and  other  more  sophisticated  and  less  expensive  rem¬ 
edies  are  desirable.  We  feel,  however,  that  radical  improvements  in  effi¬ 
ciency  should  not  be  expected. 

In  the  computer  code  that  we  have  developed,  reduction  of  the  time  step 
At  ^  is  only  performed  if  the  Newton-Raphson  iteration  fails  to  converge  in 


I 


S 


1 
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a  prescribed  number  MAXITE  of  iterations.  In  that  situation,  successively 
smaller  time  steps  At  are  tried  until  a  convergent  solution  is  obtained.  This 
smaller  time  step  is  then  kept  for  a  prescribed  number  of  steps  KSTEPR  during 
which  the  most  drastic  changes  in  the  solution  are  expected  to  occur.  After 
this  time  interval,  the  time  step  is  gradually  increased  again  to  At 

max 

It  was  found  that  for  the  smaller  time  steps,  a  reduction  below  O.lAt 

max 

was  rarely  needed  if  a  ^max  was  used  which  provided  simultaneously  for  accur¬ 
ate  computation  of  dynamic  response  and  for  the  stability  of  the  central  dif¬ 
ference  technique  and  if  the  values  of  £  used  were  not  too  small.  The  para¬ 
meter  MAXITE  has  been  set  usually  to  5  and  KSTEPR  to  a  number  in  the  range 
10—20.  An  indication  of  the  extra  computational  work  resulting  from  the 
reduced  time  steps  is  given  by  the  quotient  T’/T,  where  T  is  the  total  phys¬ 
ical  time  spanned  by  our  analysis  (with  the  time  step  reductions)  with  some 
total  number  of  time  steps  and  T*  is  the  total  physical  time  that  would 
be  possible  to  span  with  the  same  total  number  of  time  steps  if  no  reduction 
of  time  steps  were  needed.  In  most  of  the  computations,  values  of  this 
ratio  of  order  1.2  -  1.5  were  experienced.  If  more  restrictive  values  of 
£  are  used,  this  value  will  increase  signif icatively:  for  example,  for 
£  =  0.02  cm  s  ^  in  the  example  of  the  next  section  (11.1),  T'/T  =  4.2. 

11.  NUMERICAL  STUDIES 


In  this  section,  we  discuss  several  applications  of  the  models  and  me¬ 
thods  discussed  earlier  to  representative  problems. 

11.1  An  Elastic  Slab  Subjected  to  Periodic  Loading. 

Our  first  example  is  designed  more  to  test  the  performance  and  robustness 
of  the  algorithms  described  in  the  preceding  section  than  as  a  model  of  complex 


dynamic  friction  phenomena.  We  begin  with  the  analyses  of  motion  of  the  elas¬ 
tic  slab  shown  in  Fig.  51.  The  dimensions  of  the  slab  are  16  cm  x  2  cm,  as  in¬ 
dicated,  and  it  is  assumed  to  be  of  unit  thickness  and  to  be  constructed  of  a 

linearly  elastic  material  with  properties  close  to  those  of  cast  iron.  Young's 

6  3  *"1  -2 

modulus  was  assumed  to  be  E  =  1.4  x  10  (10  Kg  •  cm  •  s  ),  Poisson's  ratio 

—6  3  —  3 

V  =  0.25  and  the  mass  density  p  =  7  x  10  (10  Kg  cm  ).  The  slab  is  assumed 

to  be  in  a  state  of  plane  strain;  it  is  simply  supported  on  the  portion  F^ 

of  its  boundary  and  is  compressed  along  a  frictional  interface  by  a  flat 

"rigid"  surface,  the  vertical  downward  displacement  of  which  is  prescribed. 

This  corresponds  to  prescribing  an  initial  uniformly  distributed  gap  g  = 

-4 

-5  x  10  cm.  The  prescribed  tangential  velocity  of  the  "rigid"  surface  is 

•C 

zero  (U_  =  0) .  The  normal  contact  properties  of  the  interface  were  taken  from 

Table  1  of  BACK,  BURDEK1N  and  COWLEY  [1973],  assuming  that  both  surfaces  in 

contact  are  of  cast  iron  hand-scraped  with  a  surface  finish  (peak  to  valley 

distance) in  the  range  6-8  ym  .  The  coefficient  m^  is  then  equal  to  2  and, 

8  3  -.3  _2 

after  a  change  of  units,  c^  =  10  (10  Kg  •  cm  '  s  ).  The  friction  coeffi¬ 
cient  along  T  was  arbitrarily  assumed  to  be  y  =  0.3  and  independent  of 

g  3  _3  _2 

the  normal  load.  Consequently,  nVj,  =  2  and  c^  =  0.3  x  10  (10  Kg  cm  s  ). 

On  one  of  its  ends  (r_),  the  slab  is  subjected  to  a  time  dependent  uni- 

r 

formly  distributed  force 

t  =  t  sin  wt 
x 


where  t  =  30(10^Kgcm  ^  s 


and  a)  =  3  x  10  rad. 
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The  prescribed  initial  conditions  are  as  follows:  the  initial  velocities 
in  all  the  slab  are  zero  (r(0)  =  0)  and  the  initial  displacements  are  the  sta¬ 
tic  equilibrium  displacements  of  the  slab  due  to  the  normal  compression  exerted 
by  the  flat  surface  on  alone  (no  friction  on  and  no  applied  tractions 

on  rF) . 

We  observe  that,  due  to  the  normal  deformation  of  the  interface,  the  equi¬ 
librium  normal  displacements  and  normal  pressure  on  are  not  known  a  priori. 

The  initial  equilibrium  displacements  solution  is  obtained  by  solving  the  fol¬ 
lowing  system  of  nonlinear  algebraic  equations: 

K  r(0)  +  P(r (0)  =  0  (11.1) 


«NM 
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The  numerical  solution  of  this  frictionless  unilateral  contact  problem  is 
obtained  by  using  a  standard  Newton-Raphson  algorithm  analogous  to  that  dis¬ 
cussed  in  Section  10. 1. 

The  finite  elements  mesh  used  in  this  analysis  consists  of  16  nine-node 
isoparametric  quadratic  elements,  as  illustrated  in  Fig.  51.  The  regulariza¬ 
tion  parameter  e  for  the  Coulomb's  friction  law  was  taken,  successively,  to 
be  equal  to  1,  0.1  and  0.02  cm  s  The  dynamical  equations  of  this  discrete 
model  were  integrated  using  Newmark's  method,  as  discussed  earlier,  with  para¬ 
meters  8  =  0.25,  Y  =  0.5  and  a  maximum  time  step  of  At  =  10  ^  s. 

max 

The  distributions  of  normal  stresses  at  several  time  instants  obtained 
with  E  =  0.1  cm  s  \  are  shown  in  Fig.  52. 


i 
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The  distributions  of  friction  stresses  on  T  at  several  time  instants 


are  shown  in  Fig  53.  The  travelling  wave  type  evolution  of  these  stresses  can 
be  clearly  observed  on  that  figure  as  can  the  sharp  transition  between  the 


i 

m 


tV*  njK  AV  AiK'A/f 


time 

000 

500 


147 


sliding  and  adhesion  regions  on 

The  effect  of  the  regularization  parameter  £  on  the  evolution  of  the 
displacements,  velocities,  and  friction  stresses  at  the  contact  nodes  55  and 
58  are  shown  in  Figs.  54-59.  As  might  be  expected,  smaller  values  of  e 
lead  to  sharper  transitions  between  sliding  and  adhesion  states;  large  values 
of  £  smooth  out  the  fronts  in  these  solutions.  It  can  be  also  observed  that 
the  major  discrepancies  between  results  obtained  for  different  values  of  £ 
occur  in  the  tangential  displacements  and  velocities,  not  on  the  friction 

stresses.  This  results  from  the  fact  that  when  sliding  occurs,  the  value  of 
the  friction  stress  is  necessarily  y  times  the  normal  pressure.  Consequent¬ 
ly,  the  only  possible  discrepancy  is  in  the  relative  duration  of  stick  and 
slip  states  and  in  values  of  the  friction  stress  during  the  stick  state. 

Since,  in  the  present  case,  the  alternating  stick  and  slip  states  are  closely 
conditioned  by  the  variation  of  the  periodic  forcing  function,  it  is  not  sur¬ 
prising  that  the  discrepancies  are  small.  It  is  likely  that  for  other  situa¬ 
tions  the  use  of  large  values  for  £  could  produce  friction  stress  solutions 
much  less  accurate  than  in  the  present  example. 

An  examination  of  the  computed  variation  in  the  tangential  displacements 
and  velocities  leads  us  to  the  conclusion  that  the  essential  effect  of  the  re¬ 
duction  of  £  is  to  produce  more  accurate  adhesion  states.  High  values  of  £ 
(e.g.  £  =  1)  lead  to  solutions  which  only  vaguely  resemble  those  obtained  with 
smaller  values  of  £,  no  meaningful  conclusion  relative  to  stick  or  slip  being 
possible.  From  a  practical  point  of  view  it  appears  reasonable  to  choose  values 
of  £  ,  that  are  sufficiently  small  relative  to  the  order  of  magnitude  of  the 
tangential  velocities  that  occur  during  the  sliding  states.  In  this  manner, 
sliding  and  adhesion  will  be  essentially  relative  and  not  absolute  concepts. 

We  observe  that  for  tangential  velocities  on  the  order  of  1  cm  s  \  the  value 
of  £  =  1  cm  s  ^  is  obviously  inadequate  while  results  obtained  with  £  =  0.1, 
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Figure  55.  Tangential  displacements  at  Bode  853  using  different  values  of  e. 


t  values  of  e. 
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Friction  stresses  at  Bode  85,  using  different  values  of  e. 
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0.02  cm  s  seem  to  be  physically  reasonable. 

11.2  Steady  Sliding  and  Dynamic  Stability  of  a  Rigid  Block 


As  a  prelude  to  examples  dealing  with  nonlinear  dynamics  of  a  model  of  a 
continuous  body,  we  shall  now  consider  applications  of  the  algorithms  developed 
in  Section  10  to  some  very  simple  dynamical  systems.  Our  mission  here  is  to 
construct  very  simple  models  of  stick-slip  motion,  consistent  with  the  conti¬ 
nuum  models  discussed  previously,  that  have  sufficient  degrees  of  freedom  to 
capture  qualitatively  dynamic  friction  phenomena  of  the  type  observed  experi¬ 
mentally  by  TOLSTOI  t 1967 ]  and  others:  in  stick-slip  motions  of  sliding  bodies, 
tangential  jumps  occur  in  strict  syncronization  with  jumps  normal  to  the  sliding 
surfaces. 

For  this  purpose,  a  minimum  of  three  degrees-of-freedom  seems  to  be  neces¬ 
sary;  two  (normal  and  tangential)  is  inadequate  since  this  would  not  provide 
the  coupling  between  tangential  and  normal  motions  observed  in  Part  I.  A 
simple  mechanical  system,  not  unlike  many  seen  in  friction  experiments,  that 
can  capture  this  coupling  is  a  simple  rigid  body  in  plane  motion,  the  rotational 
degree-of-freedom  providing  the  coupling  mechanism.  Such  a  system  is  shown  in 
Fig.  60  and  61  where  we  have  indicated  as  degrees-of-freedom  the  sliding  and  pene¬ 
trating  (normal)  displacements  u^  and  u  ^  and  the  rotation  Uq. 

Assuming  small  rotations  and  applying  the  same  contact-interface  laws  used 
in  our  continuum  interface  model,  we  arrive  at  a  (regularized)  system  of  equa¬ 
tions  for  this  model  which  is  again  of  the  form, 


M  r (t)  +  C  r (t)  +  K  r(t) 


+  P(r (t) )  +  J(r(t) ,r (t) )  =  F 


(11.2) 


where 
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Figure  60.  A  Block  Sliding  with  Friction  on  a  Moving  Belt, 
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Here  M  Is  the  mass  of  the  body,  W  Its  weight,  I  Its  moment  of  inertia  relative 

to  the  mass  center,  Cx>  C^,  Cg  are  damDine  coefficients  introduced  to  model  dissipation 
effects  if  needed,  Kx  the  (tangential)  spring  stiffness,  Bq  is  the  trans¬ 
verse  width  of  the  contact  surfaces  L  is  the  length  of  the  sliding  surface 
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in  the  direction  of  sliding,  H/2  the  height  of  the  center  of  mass  above  F  , 

Kj 

and  the  remaining  terms  are  recognizable  from  our  discussions  in  Sections  6-10. 

The  penetrating  approach  at  a  point  s  on  the  contact  surface  is  a(s)  = 

(u  suQ) , 

yG  0  + 


We  now  use  this  to  study  the  dynamic  stability  of  steady  sliding.  By 
ste?  /  sliding  we  mean,  with  reference  to  Fig.  60  ,  the  static  equilibrium 
position  of  the  body  in  frictonal  contact  with  the  moving  belt.  The  equili¬ 
brium  position  is  characterized  by  the  equilibrium  equation  (set  r  =  r  =  0  in 
equation  (11.2)): 


K  r  +  P(rn)  +  J (r  ,  0)  =  F  (11.4) 

During  stick-slip  oscillations,  the  body  does  not  stay  in  this  equilibrium 
position,  although  at  the  conclusion  of  the  "stick  part"  of  a  cycle,  one  would 
expect  the  body  to  be  in  a  state  very  close  to  such  an  equilibrium  position. 

The  questions  that  then  arise  are  the  following:  does  a  solution  for  equa¬ 
tion  (3)  exist?;  is  it  unique?;  how  can  it  be  computed?;  is  it  dynamically 
stable?  We  will  address  here  only  the  last  two  questions.  While  we  hope  to 
consider  some  questions  of  existence  and  uniqueness  of  the  general  problem 
in  a  later  paper,  our  numerical  results  are  physically  reasonable  and  were 
obtained  without  computational  difficulty. 

In  order  to  compute  solutions  of  (11.4)  for  a  certain  range  [0,  c  ^,1 
of  values  of  the  parameter  c;„,  we  subdivide  the  interval  [0,  c„]  into  a 
prescribed  number  NINCT  of  increments  Ac^  =  and,  for  each  increment, 

K  =  1,  NINCT,  we  again  solve  the  nonlinear  system  of  equations  (11.4)  by  the 

Newton-Raphson  method.  For  the  Newton  process,  starting  values  for  the 

til 

equilibrium  solutions  at  the  end  of  the  K  increment  are  taken  to  be  the 
(converged)  solution  for  the  (K-l)t^1  increment.  The  initial  solution  at  c^  =  0 
obviously  corresponds  to  a  frictionless  unilateral  contact  problem  of  the  type 


in  (11.1). 
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Having  computed  the  steady  sliding  equilibrium  position  for  each  value 
of  c^,,  we  now  explore  the  sensitivity  of  these  equilibrium  states  to  small 
perturbations.  If  in  equation  (11.2)  we  let 
r(t)  =  rQ  +  w(t) 


and  subtract  the  equilibrium  equation  (11.4)  from  the  resulting  equation, 
we  obtain: 


Mw  +  Cw  +  Kw  +  [P(rQ  +  w)  -  P(rQ)]  +  [J(rQ  +  w,  w)  -  J(rQ,  0)]  = 

(11.5) 

Assuming  that  the  perturbation  w(t)  is  small,  we  consider  only  the  linear  con¬ 
tributions  for  the  terms  in  brackets  in  the  above  equation: 


P(r0  +  w>  -  Hr0)  =  || 


rro 


w  Kn(rQ)w 


J(r  +  w,  w)  -  J(r  ,  0)  = 

*v  U  "V  ~  vJ 


u 


,  12 
r  +af 


~  ~0 
r=0 


w 


r=r 

t  to 

r=0 


KTn(r0,  +  CT(rQ,  0)w 


(11.6) 


__  Trt  m 

The  matrices  KT,  K  amd  C  are  precisely  the  three-degrees-of-freedom 

analogs  of  the  finite  element  matrices  derived  in  Section  10.  It  also  turns 

T 

out  that  the  matrix  C  (r^,  0)  is  null  because,  in  the  equilibrium  positions 
of  the  body,  the  velocity  of  the  body  is  zero  but  the  prescribed  velocity 
of  the  moving  belt  is  different  from  zero  | |  >  e >  0.  This  means  that  the 
equilibrium  points  correspond  to  velocities  on  the  horizontal  branches  of  the 
friction  law  depicted  in  Fig.  50  and  consequently  the  derivative  of  the  friction 


o  > 


stresses  with  respect  to  the  sliding  velocity  is  zero.  This  also  means 
that  the  linearization  of  the  equations  of  motion  about  the  equilibrium 
position  is  independent  of  the  regularization  approximation  used:  we  are 
outside  neighborhoods  of  the  discontinuity  of  Coulomb's  friction  law. 

The  final  system  of  linearized  equations  for  the  small  oscillations 
about  the  equilibrium  position  is  then 

Mw  +  Cw  +  [K  +  Kn  +  KTn]w  =  0  (11.7) 

The  behavior  of  this  linear  system  is  determined  by  the  solutions  of 
the  following  quadratic  eigenvalue  problem 

[A2M  +  AC  +  IK  +  Kn  +  KTn]]W  =  0  (11.8) 

This  system  differs  from  standard  structural  dynamics  eigenvalue  prob¬ 
lems  because  of  the  absence  of  symmetry  due  to  friction  effects  entering 
the  non-symmetric  structure  of  the  matrix  K  .  To  illustrate  numerically 
what  can  happen  due  to  this  lack  of  symmetry,  we  compute  the  equilibrium 
positions  (11.4)  of  a  block  sliding  on  a  moving  surface,  as  shown  in  Fig. 60, 
for  different  values  of  c^  and  for  each  of  these  positions  we  compute  the 
corresponding  eigenvalues  for  (11.7). 

The  data  used  is  close  to  the  experimental  setting  devised  by  BELL  and 

BURDEKIN  [1969-70].  This  does  not  mean  that  we  will  be  trying  to  reproduce 

any  experimental  results:  several  geometrical  features  are  distinct  and 

other  data  (n^,  cT»  y,  I)  is  arbitrarily  assumed  or  approximated  from  tables 

provided  by  other  authors  (e.g.  c  ,  m  ).  For  the  time  being,  we  are  pri- 

n  n 

marily  concerned  with  qualitative  behaviors  of  this  simple  system.  The 
data  used  for  our  numerical  experiments  is  the  following: 


K  =  2388 


450 


10  Kg  s 


Kg 


M  9  J 

~  (IT  +  H  )  =  124.2 


0.0 


C  =  0.0;  57.16;  571.6;  5715 


3  2 

10  Kg  cm 


3  -1 

10  Kg  s 


ce  = 
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n 
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H  = 
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0.0;  1.388  x  10A;  1.388  x  105;  1.388  x  106 
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10 


in,  -  2.5 
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n 
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in3  v  -3.5  -2 

10  Kg  cm  s 


y  e  [o,  1.6) 


3  -2 

10  Kg  cm  s 


3  2-1 

10  Kg  cm  s 


>  (H.9) 


We  observe  that  the  geometry  of  the  problem  and  the  conditions  of 

support  of  the  body  imply  that  a  necessary  condition  for  equilibrium  is  that 

L  ,  , 
y  <  —  =  1.6. 

H 

The  values  of  the  equilibrium  displacements  u^,  u^,  and  u^  for  the 
different  values  of  y  are  plotted  in  Fig.  62.  It  can  be  seen  that  u^  in¬ 
creases  linearly  with  y  and  that  u^.  and  u^  increase  indefinitely  as  y 
approaches  1.6. 

Figure  63  is  a  plot  of  the  computed  evolution  in  the  complex  plane  of 
the  eigenvalues  associated  with  the  normal  and  rotational  degrees  of  freedom 
for  the  successive  equilibrium  positions  obtained  with  increasing  values  of  y. 
In  the  absence  of  normal  and  rotational  damping,  these  eigenvalues  are  found 
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Figure  63.  Evolution  of  the  Eigenvalues  on  the  Complex  Plane 
for  theSuooessive  Equilibrium  Positions  Obtained 
with  Increasing  Coefficient  of  Eviction  (No  damping) 
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to  be  pure  imaginary  and  distinct  for  y  *  0.  As  p  is  increased,  at  some 
value  of  y  they  coalesce  and  start  having  real  components  symmetric  with 
respect  to  the  imaginary  axis.  At  a  larger  value  of  y  they  become  pure  im¬ 
aginary  again  and  a  conjugate  pair  of  the  eigenvalues  approaches  the  point 
(0,  0)  along  the  imaginary  axis  as  y  1.6.  This  behavior  as  y  -*■  1.6 
is  exactly  what  should  be  expected  since  for  y  =  1.6  the  equili¬ 

brium  of  the  body  cannot  be  sustained  and  the  body  acquires  a  free  rigid  body 
mode  (tumbling).  We  also  observe  that  for  the  geometric  and  stiffness  prop¬ 
erties  of  this  problem  the  eigenvalues  associated  with  the  tangential  degree 
of  freedom  are  uncoupled  with  the  normal  and  rotational  freedoms  and  re¬ 
main  constant  for  all  the  values  of  y. 

However,  the  essential  fact  to  be  gained  from  this  analysis  is  that  the 
non-symmetry  of  the  frictional  contact  laws  can  produce,  for  sufficiently 
large  values  of  the  coefficient  of  friction,  eigenvalues  with  positive  real 
part  which  determine  the  dynamic  instability  of  the  steady  sliding.  It  can 
also  be  shown  that,  due  to  the  uncoupling  of  the  tangential  mode  the  intro¬ 


duction  of  external  tangential  damping  (Cx)  would  not  stabilize  the  system ; 


for  viscous  stabilization  an  external  normal  (C  )  and  rotational  (C_)  damping 

y  0 


are  required.  This  behavior  is  made  clear  in  the  complex  plane  diagram  of  Fig.  64. 
The  essential  effect  of  such  viscous  damping  is  to  shift  the  eigenvalues 
towards  the  negative  real  axis.  For  sufficiently  high  damping,  all  the  eigen¬ 
values  associated  with  the  normal  and  rotation  freedoms,  for  all  admissable 
values  of  y,  will  have  negative  real  parts. 


A  short  note  on  the  damping  factors  £  and  £Q  in  Fig.  64  is  needed:  These 

y  y 


are  merely  indicators  of  the  amount  of  damping  introduced.  They  are  defined 


using  the  natural  frequencies  and  of  the  linearized  system  for  the 


equilibrium  position  at  y  =  0:  £  =  C  / 2Mw  £  =  CQ/2lconrt  . 
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Figure  64.  Influence  of  the  normal  an  '  ^tational  damping  on  the  eigenvalues 
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We  will  now  show  the  effects  of  the  eigenvalues  structure  of  (11.8)  on 


the  motion  of  the  rigid  body.  For  each  set  of  data,  we  perform  a  transient 


analysis  of  the  motion  of  the  body  assuming  that  the  initial  conditions  are 


the  following:  the  displacements  corresponding  to  the  steady  sliding  equi¬ 


librium  and  a  small  perturbation  on  the  velocities  —  an  upwards  translational 


velocity,  i.e.. 


r(0)  =  {0.0,  -0.01  cm  s-1,  0.0}T 


(11.10) 


We  first  consider  a  case  in  which  y  is  smaller  than  the  value  at  which 


some  eigenvalues  start  to  have  real  components.  The  data  used  is: 


M,  I,  c  ,  m  ,  m  ,  L,  H,  B  ,  W  as  in  (11.9) 
n  n  i  u 


K  =  11100 
x 


3  -2 

10  Kg  s 


(11.11) 


C  T  =  yCn 


y  =  0.15 


UC  =  0.080  cm  s-1 
x 


At  =  1  x  10  s 
max 


e  =  0.002  cm  s 


(11.12) 


In  Figs.  65,  66  and  67  we  show  phase  plane  plots  for  the  tangential 


(u  „) ,  normal  (u  _)  and  rotational  (uQ)  displacements.  It  is  clear  that  the 
xG  yG  U 


resulting  motion  is  a  small-amplitude  stable  oscillation.  Obviously,  if  we 


had  considered  some  damping  this  oscillation  would  be  damped  out  and  the 


steady  sliding  would  be  essentially  attained  after  a  period  of  time. 


A  more  interesting  situation  arises  when  y  is  such  that  some  eigenvalues 


have  positive  real  parts.  Starting  from  a  basic  data  set  and  modifying,  one  at 


a  time,  several  physical  parameters,  it  was  possible  to  assess  the  influence  of 


these  parameters  on  the  dynamic  responses  obtained  in  this  case.  The  data  used 


is  the  following  (when  more  than  one  value  is  given  for  some  parameter,  the 


first  belongs  to  the  basic  data  set): 


X.W  iwkkr*  ?av.v.v.v;«,v  , 


Figure  66.  Phase  plant  plot  for  the  normal  motion  of  the  center  of  mass 


Figure  67.  Phase  plane  plot  fov  the  votation  of  the  oentev  of  mass  (y 
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L,  H,  cn>  mT  as  in  (11.9) 

Kx  =  2388,  11100,  44410  103  Kg  s"2 

Kg 

3  2 

10  Kg  cm 

Sx  "  CJ2  M(0x  "  °*005»  0.05,  0.5,  1.6,  5 

iy  -  y2M“>y0  =  °*°4 

^0  “  C0^2Ia)0O  =  0,04 

Bq  =  14  cm 

■t  ’  “  cn  V  -  1 

Uc  -  0.127,  0.060,  0.220,  0.287  cm  s”1 

At  =  1  x  10-5  s 
max 

E/Ux  6  [0.07,  0.25] 

The  following  remarks  summarize  our  interpretation  of  the  numerical 
results: 

(i)  Due  to  the  instability  of  the  normal  and  rotational  modes,  the  normal 
and  rotational  oscillations  grow  (see  Fig.  68), 

(ii)  The  variation  of  the  normal  force  on  the  contact  produces  changes  in 
the  sliding  friction  force  which  in  turn  produce  a  tangential  os¬ 


cillation; 


)RM 
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(lii)  The  tangential  oscillation  may  then  become  sufficiently  large  that, 

for  small  values  of  the  velocity  of  the  belt  the  points  of  the  body 

X 

•  C 

on  the  contact  surface  attain  the  velocity  and  the  body  sticks  for 
short  intervals  of  time  (see  Fig.  69); 

(iv)  With  the  growing  of  the  normal  oscillation  actual  normal  jumps  of  the 
body  may  occur  (see  Fig.  68  ); 

(v)  The  repeated  periods  of  adhesion  have  the  result  of  decreasing  the 
average  value  of  the  friction  force  on  the  contact  and,  due  to  the 
absence  of  equilibrium  with  the  restoring  force  on  the  tangential 
spring,  the  tangential  displacement  of  the  center  of  mass  decreases 
(see  Fig.  70-73). 

(vi)  Then,  one  of  the  two  following  situations  may  occur: 

(a)  the  normal,  rotational  and  tangential  oscillations  evolve  to  what 
appears  to  be  a  steady  oscillation  (see  Fig.  74)  with  successive  periods 
of  adhesion  and  sliding,  the  average  values  of  the  friction  force 

and  of  the  spring  elongation  being  smaller  than  those  corresponding 
to  the  steady  sliding  equilibrium  position  (see  Fig.  70,  U^=0.287  cm  s  "S ; 

(b)  at  a  sufficiently  small  value  of  the  spring  elongation  the  normal 
and  rotational  damping  are  able  to  damp  out  the  corresponding  os¬ 
cillations  and  the  body  sticks  (see  e.g.  Fig.  72)  since  the  re¬ 
storing  force  of  the  spring  is  then  smaller  than  the  maximum 
available  friction  force. 

Thus  monitoring  the  spring  elongations  as  is  usually  done  in  friotion 
experiments ,  case  (a)  would  be  perceived  as  an  apparently  smooth  sliding 
with  a  coefficient  of  kinetic  friction  smaller  than  the  coefficient  of  static 
friction  and  case  (b)  would  be  perceived  as  stick-slip  motion. 
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Figure  70.  Influence  of  the  velocity  of  the  support  on  the  tangential  motion  of 
the  center  of  mass  (y  =  1), 


Figure  71.  Influence  of  the  tangential  damping  on  the  tangential  motion  of  the 
centev  of  mass  (  y  =  !)• 
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Figure  74.  Phase  -plane  plot  of  the  steady  normal  oscillation  obtained 

during  apparently  smooth  sliding  with  an  apparent  coefficient 
of  kinetic  friction  tower  than  the  coefficient  of  static  friction 

(y  =  1,  if  =  0.287  cm  s"1) 


i 


B 

I 

B 

1 


i 


33 

■ 


179 


Figure  70  shows  that,  as  ohseryed  experimentally,  an  increase 

•  Q 

on  the  velocity  produces  an  increase  of  the  frequency  of  the  stick- 
slip  motion  and  a  decrease  of  the  tangential  displacement  amplitude 
(for  our  data  this  decrease  is  small}.  It  also  shows  that,  as 
expected,  the  transition  from  stick-slip  motion  to  an  ap¬ 
parently  smooth  sliding  occurs  abruptly  at  some  critical  velocity  0.240 

cm  s  1) .  However,  for  the  data  that  we  have  used,  that  transition  occurs 
when  the  stick-slip  frequency  (  ~  1.5  Hz)  is  excessively  small  compared  with 
the  tangential  natural  frequency  of  the  body  (  ~  11.6  Hz). 

Figure  71  shows  that  an  increase  in  tangential  damping  produces  a 
(small)  decrease  of  the  displacement  amplitude  and  that  sufficiently  large 
tangential  damping  produces  an  apparently  smooth  sliding  with  an  apparent 
coefficient  of  kinetic  friction  smaller  than  the  coefficient  of  static  fric¬ 
tion  (in  agreement  with  observations  of  TOLSTOI  [1967]  mentioned  earlier). 

Figures  72  and  73  show  that,  also  as  expected,  decreases  of  the  tan¬ 
gential  displacements  amplitude  can  be  produced  by  increasing  the  tangential 
stiffness  or  by  decreasing  the  mass. 


11.3  A  Model  for  the  Normal  Contact  Dissipation 


In  Section  (6.1),  we  remarked  that  if  the  sliding  is  not  smooth  and 

-4 

small  ( ~ 10  cm)  normal  jumps  occur,  it  is  possible  that  some  plastic  defor¬ 
mation  of  the  asperities  will  take  place.  In  the  numerical  experiments 
presented  in  the  previous  section,  it  appears  that  damping  effects  asso¬ 
ciated  with  normal  and  rotational  degrees  of  freedom  may  play  a  role  in  the 
deceleration  phase  of  the  sliding  portion  of  the  stick-slip  cycles.  In  the 
previous  section,  we  assumed  arbitrarily  that  all  the  damping  effects  dif¬ 
ferent  drom  dry  friction  could  be  simulated  by  linear  viscosity.  A  conse¬ 
quence  of  that  assumption  is  the  fact,  visible  e.g.  in  Fig.  68,  that  normal 
and  rotational  damping  act  not  only  during  contact  but  also  while  the  body 
jumps.  On  the  other  hand,  trying  to  formulate  a  detailed  model  of  the  normal 
elasto-plastic  behavior  of  the  asperities  to  apply  in  dynamic  situations  ex¬ 
pected  to  have  not  too  severe  dissipative  effects  does  not  seem  to  be  worth¬ 
while. 

We  can  borrow  one  idea  for  a  simple  model  of  such  effects  from  an  anal¬ 
ysis  of  vibroimpact  phenomena  and  impact  damping  by  HUNT  and  CROSSLEY  [1975]. 
The  key  point  is  to  model  hysteresis  loops  of  the  form  presented  in  Fig. 75(a), 
which  may  occur  if  the  surfaces  in  contact  are  allowed  to  unload  completely 
(THORNLEY  et  al.  [1965]),  by  approximate  loops  of  the  form  presented  in 
Fig.  75(b). 

If  the  normal  elastic  pressure-penetration  curve  is  of  the  power  law 
m 

type  (p  =  c^  a  n)  possible  additional  dissipative  terms  that  can  produce 
loops  of  the  type  indicated  in  Fig. 75(b)  can  be  taken  of  the  form 


Figure  75.  Hysteresis  loops  for  the  normal  deformation  of  the  interface 
(schematic). 

(a)  Experimentally  observed  loop ,  under  quasi-static  loading 
conditions 3  if  total  separation  of  the  surfaces  is  allowed 
during  unloadings. 

(b)  Hysteresis  loop  modelled  by  t'he  constitutive  equation  (11.13) 
under  dynamic  loading  conditions. 


where  a  is  the  penetrating  approach  and  is  a  nonlinear  damping  coefficient 
to  be  determined  experimentally. 

HUNT  and  CROSSLEY  [1975]  do  propose  for  macroscopic  Hertzian  contacts 
slight  generalizations  of  the  term  adopted  here  and  they  also  show  that  for 
the  case  2,_^  =  m^  and  small  energy  losses  the  correlation  between  the  energy 
loss  per  cycle  of  contact  and  the  damping  coefficient  b^  is  very  simple. 

With  this  new  term,  the  normal  contact  law  (6.4)  for  our  continuum 


model  becomes 


m  Si 

-a  =  c  (u  -g)  n  +  b  (u  -g)  "  u 
n  nn  +  nn  +  n 


(11.13) 


Accordingly  a  new  term  must  be  included  on  the  left  hand  side  of  the 


variational  statements  (8.1,  8.13): 


.  x, 

<Q(u,  u),  v>  =  J  b  (u  -g)n  u  v  ds 
~  ~  ~  »nn  +  nn 


(11.14) 


Obviously  for  the  case  of  the  three-dimensional  domain,  restrictions 
must  be  placed  upon  2,^  in  order  for  the  integral  in  the  above  expression 
to  be  well  defined. 

The  above  term  does  not  introduce  any  other  complication  or  restric¬ 
tion  on  either  the  finite  elements  discretization  nor  on  the  algorithms 
described  earlier.  In  Figs. 76, 77  we  show  the  hysteresis  loops  obtained  for 
the  case  of  the  rigid  body  model  with  the  geometric,  mass  and  contact  prop¬ 
erties  (11.9)  used  in  the  previous  section,  for  typical  values  of  the  im¬ 
pact  velocities  u^,  two  values  for  the  constant  b  and  2,  =  m  -  2.5. 

^  y’  n  n  n 

we  also  have  checked  that  complex  plane  diagrams  of  the  type  of  (Fig. 

64)  could  also  be  obtained  with  this  type  of  damping  and  that  the  typical 
stick-slip  motion  curves  of  the  type  of  Figs.  70-73  can  also  be  obtained  with 


this  more  realistic  non-linear  damping  term. 
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APPROACH  (10~*cm) 
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11.  4.  The  Steady  Sliding  of  a  Linearly  Elastic  Body  and  its  Dynamic 
Stability 


All  the  computations  performed  in  Section  11.2  with  a  rigid  body  can  also 
be  done  with  finite  element  models  of  linearly  elastic  bodies.  The  essential 
effects  observed  there  are  also  observed  here. 

In  order  to  compare  the  results  of  both  the  computations  we  consider 
the  block  of  Fig.  60  as  a  linearly  elastic  body  and  assume  a  plane  strain 
state.  The  geometry,  total  mass,  total  tangential  stiffness  and  contact 
properties  are  those  given  by  the  data  (11. 9).  The  Young’s  modulus  is  E  = 

O  1  _o 

1.4  x  10  (10  Kg  cm  s  )  and  the  Poisson's  ratio  is  V  =  0.25. 

The  finite  elements  mesh  used  consists  of  4  x  3  nine  mode  isopara¬ 
metric  elements.  In  Fig.  78  we  show  the  deformed  mesh  configurations  for 
the  steady  sliding  equilibrium  positions  at  several  values  of  the  coefficient 
of  friction  y.  As  expected,  for  the  level  of  forces  in  presence,  the  block 
behaves  much  like  a  rigid  body. 

As  in  section  11.2  we  have  also  computed  the  eingenvalues  of  (11.8)  for  the 
successive  equilibrium  positions  at  values  of  y  6  [0,  1.6),  under  the  assump¬ 
tion  of  no  damping.  In  Fig.  79  we  compare  the  rigid  body  model  eigenvalues 
associated  with  the  normal  and  rotation  displacements  with  those  from  the 
finite  elements  model  which  are  associated  with  similar  modes.  The  results 
are  very  close,  showing  that  the  def ormability  of  the  body  reduces  slightly  the 
imaginary  components.  All  the  eigenvalues  in  the  first  quadrant  of  the  com¬ 


plex  plane  are  plotted  in  Fig.  80  .  It  can  be  observed  that  in  addition  to 
those  corresponding  to  the  rigid  body  normal  and  rotational  displacements, 
there  are  other  unstable  eigenvalues  with  larger  imaginary  components  and 
that  one  of  them  does  appear  for  values  of  y  below  the  value  at  which  the  rigid 
body  instability  initiates. 
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Figure  79.  Comparison  of  the  eigenvalues  associated  with 
"normal  and  rotational  displacements"  obtained 
with  the  rigid  body  and  the  finite  elements 
model . 
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Similar  computations  were  also  performed  for  a  case  where  no  rigid 

body  motions  of  the  elastic  body  could  occur.  The  body  is  a  slab  similar 

to  the  one  presented  in  Fig.  51  .  The  material  and  normal  contact  properties 

and  the  boundary  conditions  on  are  the  same  as  those  in  section  11.1. 

The  dimensions  are  now  16  cm  x  1  cm  and  the  finite  elements  mesh  consists  of 

9x2  nine  node  isoparametric  elements.  The  initial  negative  normal  gap  on 

-4  -3 

rc  was  taken  successively  to  be  5  x  10  cm  and  1  x  10  cm.  The  "rigid" 
flat  surface  that  compresses  the  body  along  is  assumed  to  have  a  velocity 
towards  the  right.  Deformed  mesh  configurations  and  distributions  of  normal 
stresses  on  Tc  are  shown  in  Figs.  81  and  82  for  g  *  -1x10  cm.  The  effect  of 
the  increase  of  the  compression  on  the  real  parts  of  the  eigenvalues  is  shown 
in  Fig  83  :  larger  compression  produces  increased  instability. 

Our  final  numerical  examples  consist  of  obtaining  the  dynamic  response 
of  the  block  of  Fig.  60  assuming  that  it  is  linearly  elastic  and  that  it  starts 
from  the  perturbed  equilibrium  position  considered  earlier  in  section  11.2. 

In  these  examples  we  will  use  the  central  difference  method  with  a  dia¬ 
gonalized  mass  matrix.  The  maximum  time  step  used  is  At  x  =  3  x  10  ^  s. 

The  geometry,  normal  contact  properties  and  total  mass  are  again  as  in 

(11.  9).  We  will  now  assume  the  total  tangential  stiffness  Kx  eaual  to  11100 
3  -2 

(10  Kg  s  ),  the  coefficient  of  friction  y  equal  to  0.15  or  0.60  and  the 

•  P  —  1 

velocity  successively  equal  to  0.01,  0.08,  0.80  cm  s  .  No  damping  ef¬ 
fects  will  be  considered  when  modelling  the  interior  of  the  linearly  elastic 
body. 

We  observe  that  for  y  =  0.15  all  the  eigenvalues  are  still  imaginary 
(in  the  absence  of  any  damping) .  Perturbing  the  equilibrium  position  a  small 
stable  oscillation  is  obtained  exactly  as  in  the  rigid  body  case.  For  this 
reason  we  don't  reproduce  here  the  corresponding  plots. 
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Figure  81.  Undeformed(X) 3  and  deformed  mesh  configurations  of  a  compressed  stab 
in  steady  sliding  equilibrium  for  y=0(*),  y=0.5(a),  and  y=1.0(o)* 
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Length  along  Fc(cm) 

g  82.  Distribution  of  normal  stresses  on  V  for  the  steady  sliding  equilibrium  at 
various  values  of  y  .  L 
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For  the  case  y  =  0.6  (some  eigenvalues  with  positive  real  parts)  we 
have  assumed  normal  dissipation  effects  on  the  contact  boundary  of  the  form 

TO  9  /  C  1 

described  in  section  11.3  with  b  *  0.381  x  10  10  Kg  cm  *  s  and 

n 

»  2.5.  The  description  of  the  motion  of  the  body  is  essentially  the  same 
given  for  the  rigid  body  case. 

The  evolution  of  the  elongation  of  the  spring  is  shown  in  Fig.  84  for 
•C 

the  three  velocities  U  considered.  The  decrease  of  the  amplitude  of  the 

•C 

stick-slip  motion  with  the  increase  of  U  is  again  observed.  In  Fig.  85 

x 

we  show  a  phase  plane  plot  of  the  normal  oscillations  of  the  contact  node  29 
and  in  Figs.  86  and  87  we  show  the  evolution  of  the  normal  and  friction 
stresses  on  the  same  node.  In  order  to  assure  the  reader  that  the  numerous 
spikes  in  Figs.  86  and  87  are  not  erroneous  we  show  in  Figs.  88 
and  89  what  one  of  those  spikes  may  look  like.  In  Fig.  89  the  occurrence 
at  each  cycle  of  contact  of  periods  of  adhesion  and  sliding  is  indicated. 

It  is  also  important  to  observe  that  having  used  values  for  the  normal  con¬ 
tact  properties  (c  ,  m  )  taken  from  Table  1  of  BACK,  BURDEKIN  and  COWLEY 

n  n 

[1973]  and  having  considered  a  block  with  dimensions  and  weight  close  to  those 

used  by  BELL  and  BURDEKIN  [1969-70],  it  turns  out  that  the  frequency  of  the 

normal  contact  oscillations  observed  in  Figs.  85  ,  86  and  88  is  of  the  order 

3 

of  magnitude  indicated  by  TOLSTOI  [1967]  as  typical:  10  Hz. 

•  C  —1 

For  the  case  of  the  larger  velocity  (U  =  0.8  cm  s  )  that  velocity  is 

X 

sufficiently  large  that  the  tangential  oscillation  of  the  body  is  not  suffi¬ 
cient  to  produce  any  stick  state.  Consequently  the  average  coefficient  of 
friction  during  sliding  is  equal  to  the  static  coefficient  of  friction.  We 
note  also  that  the  instability  of  the  equilibrium  position  makes  it  impossi¬ 
ble  for  the  contact  damping  to  damp  out  the  normal  oscillation.  A  steady  self- 
excited  oscillation  is  then  attained.  That  can  be  observed  in  Fig.  89  (hor¬ 
izontal  oscillation  of  the  node  connected  to  the  spring)  and  in  Fig.  90 


Figure  84.  Evolution  of  the  spring  elongation  for  different 


Figure  85(b).  Phase  plane  plots  of  the  normal  oscalllation  of  the  contact 
node  29.  Pinal  decreasing  oscillation. 
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Figure  88.  Evolution  of  the  normal  contact  stress  at  node  29 
during  the  contact  portion  of  a  cycle  of  normal 
oscillation  (\f  =  0.08  cm  s-1^. 
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Figure  89.  Evolution  of  the  friction  stress  at  node  29  during 

the  contact  portion  of  a  cycle  of  normal  oscillation. 
($  =  0.08  cm  s-1) . 
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(normal  oscillation  of  the  contact  node  29).  In  our  computations  we 
have  observed  these  self -excited  oscillations  (also  in  the  rigid  body  case) 
whenever  the  velocity  is  sufficiently  large  compared  with  the  critical 

X 

velocity  at  which  the  transition  from  (global)  stick-slip  motion  to  (appar¬ 
ently)  smooth  sliding  occurs  and  the  coefficient  of  friction  is  large  enough 
that  the  steady  sliding  equilibrium  position  is  not  stable. 

12.  SUMMARY,  CONCLUSIONS  AND  FINAL  COMMENTS 

In  this  paper  a  large  body  of  literature  on  the  frictional  contact  of 
metallic  surfaces  was  reviewed.  Particular  emphasis  was  given  to  the  normal 
def ormability  of  rough  interfaces,  the  transition  between  static  and  kinetic 
friction,  the  occurrence  of  stick-slip  motion  and  the  fundamental  role 
played  by  normal  contact  oscillations  on  the  frictional  behavior. 

To  provide  some  focus  for  this  investigation,  we  have  proposed  that  dry 
metallic  friction  phenomenon  be  divided  into  (at  least)  three  basic  categories 
which  perhaps  require  different  theories  for  their  characterization:  type  I  - 
quasi-static  plastic  deformation  of  interfaces,  type  II  -  dynamic  friction 
phenomena  of  materially  stable  interfaces,  and  type  III  -  cases  in  which 
plowing  and  wear  are  significant.  The  present  investigation  focuses  on  type 
II  -  friction. 

Simple  continuum  mechanics  models  for  normal  and  tangential  interface 
response  which  reflect  essential  experimental  observations  were  formulated. 
Variational  principles  which  govern  a  large  class  of  frictional  contact  prob¬ 
lems  in  elastodynamics  were  also  presented. 

Numerical  algorithms  for  the  analysis  of  finite  element  approximations 
of  the  continuum  problems  formulated  earlier  were  also  proposed.  Several 
numerical  experiments  show  the  feasibility  of  the  techniques  proposed. 
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The  dynamic  stability  of  steady  sliding  and  its  consequences  on  the  be¬ 
havior  of  sliding  bodies  ("rigid"  or  linearly  elastic)  was  numerically  studied 
Occurrence  of  apparently  smooth  sliding  with  an  apparent  coefficient  of  kin¬ 
etic  friction  smaller  than  the  static  one  and  occurrence  of  stick-slip  oscilla 
tions  were  numerically  demonstrated.  These  were  shown  to  be  a  consequence  of 
dynamic  instabilities  which  occur  for  sufficiently  large  values  of  the  coef¬ 
ficient  of  (static)  friction  and  appropriate  values  of  the  other  physical 
parameters  involved.  Those  instabilities  were  shown  to  be  the  consequence  of 
the  inherent  non-symmetry  of  frictional  contact  problems  and  they  may  occur 
even  if  the  coefficient  of  kinetic  friction  is  assumed  to  be  equal  to  the 
coefficient  of  static  friction. 

Preliminary  numerical  results  on  the  effect  of  the  variation  of  several 
physical  parameters  on  the  qualitative  behavior  of  the  sliding  at  small  velo¬ 
cities  give  evidence  that  our  models  are  very  promising  simulators  of  these 
classes  of  dynamic  response. 

The  experimental  facts  summarized  in  this  work  and  the  numerical  results 
obtained  here  do  have  implications  for  both  the  experimental  and  numerical 
study  of  frictional  dynamic  problems. 

Experimental  results  on  the  normal  def ormability  of  metallic  interfaces 
at  low  contact  pressures  for  all  the  combinations  of  materials  and  surface 
finishes  of  interest  in  practical  applications  are  desperately  needed.  A 
thorough  test  of  the  models  proposed  here  will  only  be  possible  with  close  co¬ 
operation  between  experimentalists  and  numerical  analysts. 

The  dispersion  of  the  values  reported  in  the  literature  for  the 
coefficient  of  friction  is  not  too  surprising  in  view  of  our  results. 

Our  results  and  the  experiments  of  several  authors  (a.g.  TOLSTOI 
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[1967])  strongly  suggest  that  the  (apparently)  smooth  sliding 
motion  and  the  corresponding  apparent  value  of  the  coefficient  of  kinetic 
friction  are  highly  dependent  on  the  geometric  and  dynamic  properties  of 
the  apparatus  used  (stiffness,  mass,  damping).  The  following  observations 
of  MADAKSON  [1983]  are  easily  understandable  in  the  above  context:  "It  has 
been  demonstrated  that  the  friction  of  a  given  material  depends  also  on  the 
test  system.  Samples  of  an  identical  material  were  distributed  to  different 
laboratories  to  measure  the  friction  at  given  conditions.  Using  different 
measuring  systems  each  laboratory  reported  a  different  value  of  the  friction." 

The  normal  def ormability  of  the  interface  is  an  essential  feature  of 
dynamic  contact  problems  involving  metallic  bodies.  Not  taking  into  account 
this  def ormability  in  finite  elements  models  of  these  phenomena  (even  if 
frictional  effects  are  negligible)  leads  to  some  serious  physical  inconsis¬ 
tencies.  For  example,  the  absence  of  normal  def ormability  leads  to  models 
which  can  provide  oscillations,  depending  on  the  mesh  used,  with  frequencies 

as  high  as  10^  Hz  but  these  models  may  be  incapable  of  delivering  experimen- 

3 

tally  observed  contact  oscillations  of  frequencies  of  the  order  of  10  Hz. 

We  believe  that  some  difficulties  encountered  in  practical  finite  element 
calculations  of  dynamic  friction  problems  (rapidly  alternating  contact-re¬ 
lease  and  stick-slip  situations)  have  the  same  origin  as  the  oscillations 
observed  in  our  numerical  computations  —  inherent  non-symmetry  of  contact 
friction  contributions  to  the  governing  dynamical  equations  and  consequent 
dynamic  instabilities. 

These  effects  are  certainly  unrealistically  exaggerated  by  the  use  of 
an  infinitely  rigid  normal  contact  interface  response  and  a  programmed 
(artificial)  discontinuous  drop  of  the  friction  coefficient  from  a  static 
value  to  a  lower  kinetic  value  upon  sliding. 
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